The salmonid MHC class I: more ancient loci uncovered

Salamanders reveal novel trajectories of amphibian MHC evolution

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

Intra-specific copy number variation of MHC class II genes in the Siamese fighting fish

Duplicates of genes for major histocompatibility complex (MHC) molecules can be subjected to selection independently and vary markedly in their evolutionary rates, sequence polymorphism, and functional roles. Therefore, without a thorough understanding of their copy number variation (CNV) in the genome, the MHC-dependent fitness consequences within a species could be misinterpreted. Studying the intra-specific CNV of this highly polymorphic gene, however, has long been hindered by the difficulties in assigning alleles to loci and the lack of high-quality genomic data. Here, using the high-quality genome of the Siamese fighting fish (Betta splendens), a model for mate choice studies, and the whole-genome sequencing (WGS) data of 17 Betta species, we achieved locus-specific amplification of their three classical MHC class II genes - DAB1, DAB2, and DAB3. By performing quantitative PCR and depth-of-coverage analysis using the WGS data, we revealed intra-specific CNV at the DAB3 locus. We identified individuals that had two allelic copies (i.e., heterozygous or homozygous) or one allele (i.e., hemizygous) and individuals without this gene. The CNV was due to the deletion of a 20-kb-long genomic region harboring both the DAA3 and DAB3 genes. We further showed that the three DAB genes were under different modes of selection, which also applies to their corresponding DAA genes that share similar pattern of polymorphism. Our study demonstrates a combined approach to study CNV within a species, which is crucial for the understanding of multigene family evolution and the fitness consequences of CNV.

MHC class I evolution; from Northern pike to salmonids

Background

Salmonids are of major importance both as farmed and wild animals. With the changing environment comes changes in pathogenic pressures so understanding the immune system of all salmonid species is of essence. Major histocompatibility complex (MHC) genes are key players in the adaptive immune system signalling infection to responding T-cells populations. Classical MHC class I (MHCI) genes, defined by high polymorphism, broad expression patterns and peptide binding ability, have a key role in inducing immunity. In salmonids, the fourth whole genome duplication that occurred 94 million years ago has provided salmonids with duplicate MHCI regions, while Northern Pike, a basal sister clade to salmonids, represent a species which has not experienced this whole genome duplication.

Results

Comparing the gene organization and evolution of MHC class I gene sequences in Northern pike versus salmonids displays a complex picture of how many of these genes evolved. Regional salmonid Ia and Ib Z lineage gene duplicates are not orthologs to the Northern pike Z lineage sequences. Instead, salmonids have experienced unique gene duplications in both duplicate regions as well as in the Salmo and Oncorhynchus branch. Species-specific gene duplications are even more pronounced for some L lineage genes.

Conclusions

Although both Northern pike as well as salmonids have expanded their U and Z lineage genes, these gene duplications occurred separately in pike and in salmonids. However, the similarity between these duplications suggest the transposable machinery was present in a common ancestor. The salmonid MHCIa and MHCIb regions were formed during the 94 MYA since the split from pike and before the Oncorhynchus and Salmo branch separated. As seen in tetrapods, the non-classical U lineage genes are diversified duplicates of their classical counterpart. One MHCI lineage, the L lineage, experienced massive species-specific gene duplications after Oncorhynchus and Salmo split approximately 25 MYA. Based on what we currently know about L lineage genes, this large variation in number of L lineage genes also signals a large functional diversity in salmonids.

Teleost cytotoxic T cells

Cell-mediated cytotoxicity is one of the major mechanisms by which vertebrates control intracellular pathogens. Two cell types are the main players in this immune response, natural killer (NK) cells and cytotoxic T lymphocytes (CTL). While NK cells recognize altered target cells in a relatively unspecific manner CTLs use their T cell receptor to identify pathogen-specific peptides that are presented by major histocompatibility (MHC) class I molecules on the surface of infected cells. However, several other signals are needed to regulate cell-mediated cytotoxicity involving a complex network of cytokine- and ligand-receptor interactions. Since the first description of MHC class I molecules in teleosts during the early 90s of the last century a remarkable amount of information on teleost immune responses has been published. The corresponding studies describe teleost cells and molecules that are involved in CTL responses of higher vertebrates. These studies are backed by functional investigations on the killing activity of CTLs in a few teleost species. The present knowledge on teleost CTLs still leaves considerable room for further investigations on the mechanisms by which CTLs act. Nevertheless the information on teleost CTLs and their regulation might already be useful for the control of fish diseases by designing efficient vaccines against such diseases where CTL responses are known to be decisive for the elimination of the corresponding pathogen. This review summarizes the present knowledge on CTL regulation and functions in teleosts. In a special chapter, the role of CTLs in vaccination is discussed.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Major histocompatibility complex (MHC) fragment numbers alone - in Atlantic cod and in general - do not represent functional variability

This correspondence concerns a publication by Malmstrøm et al. in Nature Genetics in October 2016. Malmstrøm et al. made an important contribution to fish phylogeny research by using low-coverage genome sequencing for comparison of 66 teleost (modern bony) fish species, with 64 of those 66 belonging to the species-rich clade Neoteleostei, and with 27 of those 64 belonging to the order Gadiformes. For these 66 species, Malmstrøm et al. estimated numbers of genes belonging to the major histocompatibility complex (MHC) class I lineages U and Z and concluded that in teleost fish these combined numbers are positively associated with, and a driving factor of, the rates of establishment of new fish species (speciation rates). They also claimed that functional genes for the MHC class II system molecules MHC IIA, MHC IIB, CD4 and CD74 were lost in early Gadiformes. Our main criticisms are (1) that the authors did not provide sufficient evidence for presence or absence of intact functional MHC class I or MHC class II system genes, (2) that they did not discuss that an MHC subpopulation gene number alone is a very incomplete measure of MHC variance, and (3) that the MHC system is more likely to reduce speciation rates than to enhance them. Furthermore, their use of the Ornstein-Uhlenbeck model is a typical example of overly naïve use of that model system. In short, we conclude that their new model of MHC class I evolution, reflected in their title "Evolution of the immune system influences speciation rates in teleost fish", is unsubstantiated, and that their "pinpointing" of the functional loss of the MHC class II system and all the important MHC class II system genes to the onset of Gadiformes is preliminary, because they did not sufficiently investigate the species at the clade border.

Major histocompatibility complex (MHC) fragment numbers alone - in Atlantic cod and in general - do not represent functional variability

This correspondence concerns a publication by Malmstrøm et al. in Nature Genetics in October 2016. Malmstrøm et al. made an important contribution to fish phylogeny research by using low-coverage genome sequencing for comparison of 66 teleost (modern bony) fish species, with 64 of those 66 belonging to the species-rich clade Neoteleostei, and with 27 of those 64 belonging to the order Gadiformes. For these 66 species, Malmstrøm et al. estimated numbers of genes belonging to the major histocompatibility complex (MHC) class I lineages U and Z and concluded that in teleost fish these combined numbers are positively associated with, and a driving factor of, the rates of establishment of new fish species (speciation rates). They also claimed that functional genes for the MHC class II system molecules MHC IIA, MHC IIB, CD4 and CD74 were lost in early Gadiformes. Our main criticisms are (1) that the authors did not provide sufficient evidence for presence or absence of intact functional MHC class I or MHC class II system genes, (2) that they did not discuss that an MHC subpopulation gene number alone is a very incomplete measure of MHC variance, and (3) that the MHC system is more likely to reduce speciation rates than to enhance them. Furthermore, their use of the Ornstein-Uhlenbeck model is a typical example of overly naïve use of that model system. In short, we conclude that their new model of MHC class I evolution, reflected in their title “Evolution of the immune system influences speciation rates in teleost fish”, is unsubstantiated, and that their “pinpointing” of the functional loss of the MHC class II system and all the important MHC class II system genes to the onset of Gadiformes is preliminary, because they did not sufficiently investigate the species at the clade border.

Small-scale intraspecific patterns of adaptive immunogenetic polymorphisms and neutral variation in Lake Superior lake trout

Many fishes express high levels of intraspecific variability, often linked to resource partitioning. Several studies show that a species' evolutionary trajectory of adaptive divergence can undergo reversals caused by changes in its environment. Such a reversal in neutral genetic and morphological variation among lake trout Salvelinus namaycush ecomorphs appears to be underway in Lake Superior. However, a water depth gradient in neutral genetic divergence was found to be associated with intraspecific diversity in the lake. To investigate patterns of adaptive immunogenetic variation among lake trout ecomorphs, we used Illumina high-throughput sequencing. The population's genetic structure of the major histocompatibility complex (MHC Class IIβ exon 2) and 18 microsatellite loci were compared to disentangle neutral and selective processes at a small geographic scale. Both MHC and microsatellite variation were partitioned more by water depth stratum than by ecomorph. Several metrics showed strong clustering by water depth in MHC alleles, but not microsatellites. We report a 75% increase in the number of MHC alleles shared between the predominant shallow and deep water ecomorphs since a previous lake trout MHC study at the same locale (c. 1990s data). This result is consistent with the reverse speciation hypothesis, although adaptive MHC polymorphisms persist along an ecological gradient. Finally, results suggested that the lake trout have multiple copies of the MHC II locus consistent with a historic genomic duplication event. Our findings indicated that conservation approaches for this species could focus on managing various ecological habitats by depth, in addition to regulating the fisheries specific to ecomorphs.

EGAR, A Food Protein-Derived Tetrapeptide, Reduces Seizure Activity in Pentylenetetrazole-Induced Epilepsy Models Through α-Amino-3-Hydroxy-5-Methyl-4-Isoxazole Propionate Receptors

A primary pathogeny of epilepsy is excessive activation of α-amino-3-hydroxy-5-methyl-4-isoxazole propionate receptors (AMPARs). To find potential molecules to inhibit AMPARs, high-throughput screening was performed in a library of tetrapeptides in silico. Computational results suggest that some tetrapeptides bind stably to the AMPAR. We aligned these sequences of tetrapeptide candidates with those from in vitro digestion of the trout skin protein. Among salmon-derived products, Glu-Gly-Ala-Arg (EGAR) showed a high biological affinity toward AMPAR when tested in silico. Accordingly, natural EGAR was hypothesized to have anticonvulsant activity, and in vitro experiments showed that EGAR selectively inhibited AMPAR-mediated synaptic transmission without affecting the electrophysiological properties of hippocampal pyramidal neurons. In addition, EGAR reduced neuronal spiking in an in vitro seizure model. Moreover, the ability of EGAR to reduce seizures was evaluated in a rodent epilepsy model. Briefer and less severe seizures versus controls were shown after mice were treated with EGAR. In conclusion, the promising experimental results suggest that EGAR inhibitor against AMPARs may be a target for antiepilepsy pharmaceuticals. Epilepsy is a common brain disorder characterized by the occurrence of recurring, unprovoked seizures. Twenty to 30 % of persons with epilepsy do not achieve adequate seizure control with any drug. Here we provide a possibility in which a natural and edible tetrapeptide, EGAR, can act as an antiepileptic agent. We have combined computation with in vitro experiments to show how EGAR modulates epilepsy. We also used an animal model of epilepsy to prove that EGAR can inhibit seizures in vivo. This study suggests EGAR as a potential pharmaceutical for the treatment of epilepsy.

Evolution of innate-like T cells and their selection by MHC class I-like molecules

Until recently, major histocompatibility complex (MHC) class I-like-restricted innate-like αβT (iT) cells expressing an invariant or semi-invariant T cell receptor (TCR) repertoire were thought to be a recent evolutionary acquisition restricted to mammals. However, molecular and functional studies in Xenopus laevis have demonstrated that iT cells, defined as MHC class I-like-restricted innate-like αβT cells with a semi-invariant TCR, are evolutionarily conserved and prominent from early development in amphibians. As these iT cells lack the specificity conferred by conventional αβ TCRs, it is generally considered that they are specialized to recognize conserved antigens equivalent to pathogen-associated molecular patterns. Thus, one advantage offered by the MHC class I-like iT cell-based recognition system is that it can be adapted to a common pathogen and function on the basis of a relatively small number of T cells. Although iT cells have only been functionally described in mammals and amphibians, the identification of non-classical MHC/MHC class I-like genes in other groups of endothermic and ectothermic vertebrates suggests that iT cells have a broader phylogenetic distribution than previously envisioned. In this review, we discuss the possible role of iT cells during the emergence of the jawed vertebrate adaptive immune system.

MHC and Evolution in Teleosts

Major histocompatibility complex (MHC) molecules are key players in initiating immune responses towards invading pathogens. Both MHC class I and class II genes are present in teleosts, and, using phylogenetic clustering, sequences from both classes have been classified into various lineages. The polymorphic and classical MHC class I and class II gene sequences belong to the U and A lineages, respectively. The remaining class I and class II lineages contain nonclassical gene sequences that, despite their non-orthologous nature, may still hold functions similar to their mammalian nonclassical counterparts. However, the fact that several of these nonclassical lineages are only present in some teleost species is puzzling and questions their functional importance. The number of genes within each lineage greatly varies between teleost species. At least some gene expansions seem reasonable, such as the huge MHC class I expansion in Atlantic cod that most likely compensates for the lack of MHC class II and CD4. The evolutionary trigger for similar MHC class I expansions in tilapia, for example, which has a functional MHC class II, is not so apparent. Future studies will provide us with a more detailed understanding in particular of nonclassical MHC gene functions.

A comprehensive analysis of teleost MHC class I sequences

BACKGROUND

MHC class I (MHCI) molecules are the key presenters of peptides generated through the intracellular pathway to CD8-positive T-cells. In fish, MHCI genes were first identified in the early 1990's, but we still know little about their functional relevance. The expansion and presumed sub-functionalization of cod MHCI and access to many published fish genome sequences provide us with the incentive to undertake a comprehensive study of deduced teleost fish MHCI molecules.

RESULTS

We expand the known MHCI lineages in teleosts to five with identification of a new lineage defined as P. The two lineages U and Z, which both include presumed peptide binding classical/typical molecules besides more derived molecules, are present in all teleosts analyzed. The U lineage displays two modes of evolution, most pronouncedly observed in classical-type alpha 1 domains; cod and stickleback have expanded on one of at least eight ancient alpha 1 domain lineages as opposed to many other teleosts that preserved a number of these ancient lineages. The Z lineage comes in a typical format present in all analyzed ray-finned fish species as well as lungfish. The typical Z format displays an unprecedented conservation of almost all 37 residues predicted to make up the peptide binding groove. However, also co-existing atypical Z sub-lineage molecules, which lost the presumed peptide binding motif, are found in some fish like carps and cavefish. The remaining three lineages, L, S and P, are not predicted to bind peptides and are lost in some species.

CONCLUSIONS

Much like tetrapods, teleosts have polymorphic classical peptide binding MHCI molecules, a number of classical-similar non-classical MHCI molecules, and some members of more diverged MHCI lineages. Different from tetrapods, however, is that in some teleosts the classical MHCI polymorphism incorporates multiple ancient MHCI domain lineages. Also different from tetrapods is that teleosts have typical Z molecules, in which the residues that presumably form the peptide binding groove have been almost completely conserved for over 400 million years. The reasons for the uniquely teleost evolution modes of peptide binding MHCI molecules remain an enigma.

Multiple divergent haplotypes express completely distinct sets of class I MHC genes in zebrafish

The zebrafish is an important animal model for stem cell biology, cancer, and immunology research. Histocompatibility represents a key intersection of these disciplines; however, histocompatibility in zebrafish remains poorly understood. We examined a set of diverse zebrafish class I major histocompatibility complex (MHC) genes that segregate with specific haplotypes at chromosome 19, and for which donor-recipient matching has been shown to improve engraftment after hematopoietic transplantation. Using flanking gene polymorphisms, we identified six distinct chromosome 19 haplotypes. We describe several novel class I U lineage genes and characterize their sequence properties, expression, and haplotype distribution. Altogether, ten full-length zebrafish class I genes were analyzed, mhc1uba through mhc1uka. Expression data and sequence properties indicate that most are candidate classical genes. Several substitutions in putative peptide anchor residues, often shared with deduced MHC molecules from additional teleost species, suggest flexibility in antigen binding. All ten zebrafish class I genes were uniquely assigned among the six haplotypes, with dominant or codominant expression of one to three genes per haplotype. Interestingly, while the divergent MHC haplotypes display variable gene copy number and content, the different genes appear to have ancient origin, with extremely high levels of sequence diversity. Furthermore, haplotype variability extends beyond the MHC genes to include divergent forms of psmb8. The many disparate haplotypes at this locus therefore represent a remarkable form of genomic region configuration polymorphism. Defining the functional MHC genes within these divergent class I haplotypes in zebrafish will provide an important foundation for future studies in immunology and transplantation.

Evolution by selection, recombination, and gene duplication in MHC class I genes of two Rhacophoridae species

BACKGROUND

Comparison of major histocompatibility complex (MHC) genes across vertebrate species can reveal molecular mechanisms underlying the evolution of adaptive immunity-related proteins. As the first terrestrial tetrapods, amphibians deserve special attention because of their exposure to probably increased spectrum of microorganisms compared with ancestral aquatic fishes. Knowledge regarding the evolutionary patterns and mechanisms associated with amphibian MHC genes remains limited. The goal of the present study was to isolate MHC class I genes from two Rhacophoridae species (Rhacophorus omeimontis and Polypedates megacephalus) and examine their evolution.

RESULTS

We identified 27 MHC class I alleles spanning the region from exon 2 to 4 in 38 tree frogs. The available evidence suggests that these 27 sequences all belong to classical MHC class I (MHC Ia) genes. Although several anuran species only display one MHC class Ia locus, at least two or three loci were observed in P. megacephalus and R. omeimontis, indicating that the number of MHC class Ia loci varies among anuran species. Recombination events, which mainly involve the entire exons, played an important role in shaping the genetic diversity of the 27 MHC class Ia alleles. In addition, signals of positive selection were found in Rhacophoridae MHC class Ia genes. Amino acid sites strongly suggested by program to be under positive selection basically accorded with the putative antigen binding sites deduced from crystal structure of human HLA. Phylogenetic relationships among MHC class I alleles revealed the presence of trans-species polymorphisms.

CONCLUSIONS

In the two Rhacophoridae species (1) there are two or three MHC class Ia loci; (2) recombination mainly occurs between the entire exons of MHC class Ia genes; (3) balancing selection, gene duplication and recombination all contribute to the diversity of MHC class Ia genes. These findings broaden our knowledge on the evolution of amphibian MHC systems.

Selection and phylogenetics of salmonid MHC class I: wild brown trout (Salmo trutta) differ from a non-native introduced strain

We tested how variation at a gene of adaptive importance, MHC class I (UBA), in a wild, endemic Salmo trutta population compared to that in both a previously studied non-native S. trutta population and a co-habiting Salmo salar population (a sister species). High allelic diversity is observed and allelic divergence is much higher than that noted previously for co-habiting S. salar. Recombination was found to be important to population-level divergence. The α1 and α2 domains of UBA demonstrate ancient lineages but novel lineages are also identified at both domains in this work. We also find examples of recombination between UBA and the non-classical locus, ULA. Evidence for strong diversifying selection was found at a discrete suite of S. trutta UBA amino acid sites. The pattern was found to contrast with that found in re-analysed UBA data from an artificially stocked S. trutta population.

Molecular cloning and characterization of sea bass (Dicentrarchus labrax, L.) MHC class I heavy chain and β2-microglobulin

In this work, the gene and cDNA of sea bass (Dicentrarchus labrax) β2-microglobulin (Dila-β2m) and several cDNAs of MHC class I heavy chain (Dila-UA) were characterized. While Dila-β2m is single-copy, numerous Dila-UA transcripts were identified per individual with variability at the peptide-binding domain (PBD), but also with unexpected diversity from the connective peptide (CP) through the 3' untranslated region (UTR). Phylogenetic analysis segregates Dila-β2m and Dila-UA into each subfamily cluster, placing them in the fish class and branching Dila-MHC-I with lineage U. The α1 domains resemble those of the recently proposed L1 trans-species lineage. Although no Dila-specific α1, α2 or α3 sub-lineages could be observed, two highly distinct sub-lineages were identified at the CP/TM/CYT regions. The three-dimensional homology model of sea bass MHC-I complex is consistent with other characterized vertebrate structures. Furthermore, basal tissue-specific expression profiles were determined for both molecules, and expression of β2m was evaluated after poly I:C stimulus. Results suggest these molecules are orthologues of other β2m and teleost classical MHC-I and their basic structure is evolutionarily conserved, providing relevant information for further studies on antigen presentation in this fish species.

Signatures of natural selection among lineages and habitats in Oncorhynchus mykiss

Recent advances in molecular interrogation techniques now allow unprecedented genomic inference about the role of adaptive genetic divergence in wild populations. We used high-throughput genotyping to screen a genome-wide panel of 276 single nucleotide polymorphisms (SNPs) for the economically and culturally important salmonid Oncorhynchus mykiss. Samples included 805 individuals from 11 anadromous and resident populations from the northwestern United States and British Columbia, and represented two major lineages including paired populations of each life history within single drainages of each lineage. Overall patterns of variation affirmed clear distinctions between lineages and in most instances, isolation by distance within them. Evidence for divergent selection at eight candidate loci included significant landscape correlations, particularly with temperature. High diversity of two nonsynonymous mutations within the peptide-binding region of the major histocompatibility complex (MHC) class II (DAB) gene provided signatures of balancing selection. Weak signals for potential selection between sympatric resident and anadromous populations were revealed from genome scans and allele frequency comparisons. Our results suggest an important adaptive role for immune-related functions and present a large genomic resource for future studies

Sequence analysis of MHC class I α2 from sockeye salmon (Oncorhynchus nerka)

Most studies assessing adaptive MHC diversity in salmon populations have focused on the classical class II DAB or DAA loci, as these have been most amenable to single PCR amplifications due to their relatively low level of sequence divergence. Herein, we report the characterization of the classical class I UBA α2 locus based on collections taken throughout the species range of sockeye salmon (Oncorhynchus nerka). Through use of multiple lineage-specific primer sets, denaturing gradient gel electrophoresis and sequencing, we identified thirty-four alleles from three highly divergent lineages. Sequence identity between lineages ranged from 30.0% to 56.8% but was relatively high within lineages. Allelic identity within the antigen recognition site (ARS) was greater than for the longer sequence. Global positive selection on UBA was seen at the sequence level (dN:dS = 1.012) with four codons under positive selection and 12 codons under negative selection.

Molecular characterization of MHC class I and beta-2 microglobulin in a clonal strain of ginbuna crucian carp, Carassius auratus langsdorfii

Clonal ginbuna crucian carp is, a naturally gynogenetic fish, and is a useful model animal for studying T-cell-mediated immunity. To gain molecular information on MHC class I molecules from this species, we have identified four types of MHC class I (caauUA-S3n, caauUF-S3n, caauZE-S3n, and caauZB-S3n) and five beta 2-microglobulin (β(2)m) (caauβ2m-1a, caauβ2m-1b, caauβ2m-2, caauβ2m-3a and caauβ2m-3b) by an expressed sequence tag (EST) analysis and using homology cloning with degenerated primers. Like UA class I genes in other cyprinid fish, the caauUA-S3n shows features of classical MHC class I, such as conservation of all key amino acids interacting with antigenic peptides, and ubiquitous tissue expression. A phylogenetic analysis shows that the β(2)m-1 and β(2)m-2 isoforms are clustered with those of other cyprinid fishes, while β(2)m-3 isoforms make a cluster that is separated from a common ancestor of salmonid and cyprinid fishes. This finding suggests that the β(2)m isoforms of ginbuna cruician carp comprise two lineages and may possess different functions. The MHC class I and β(2)m sequences from one clonal strain will facilitate our understanding of the interaction of MHC class I with β(2)m in teleosts.

Outlier SNPs show more genetic structure between two Bay of Fundy metapopulations of Atlantic salmon than do neutral SNPs

Atlantic salmon of Eastern Canada were once of considerable importance to aboriginal, recreational, and commercial fisheries, yet many populations are now in decline, particularly those of the inner Bay of Fundy (iBoF), which were recently listed as endangered. We investigated whether nonneutral SNPs could be used to assign individual Atlantic salmon accurately to either the iBoF or the outer Bay of Fundy (oBoF) metapopulations because this has been difficult with existing neutral markers. We first searched for markers under diversifying selection by genotyping eight captively bred Bay of Fundy (BoF) populations for 320 SNP loci with the Sequenom MassARRAY™ system and then analysed the data set with four different F(ST) outlier detection programs. Three outlier loci were identified by both BayesFST and BayeScan whereas seven outlier loci, including the three previously mentioned, were identified by both Fdist and Arlequin. A subset of 14 nonneutral SNPs was more accurate (85% accuracy) than a subset of 67 neutral SNPs (75% accuracy) at assigning individual salmon back to their metapopulation. We then chose a subset of nine outlier SNP markers and used them to inexpensively genotype archived DNA samples from seven wild BoF populations using Invader™ chemistry. Hierarchical AMOVA of these independent wild samples corroborated our previous findings of significant genetic differentiation between iBoF and oBoF salmon metapopulations. Our research shows that identifying and using outlier loci is an important step towards achieving the goal of consistently and accurately distinguishing iBoF from oBoF Atlantic salmon, which will aid in their conservation.

Retained orthologous relationships of the MHC Class I genes during euteleost evolution

Major histocompatibility complex (MHC) class I molecules play a pivotal role in immune defense system, presenting the antigen peptides to cytotoxic CD8+ T lymphocytes. Most vertebrates possess multiple MHC class I loci, but the analysis of their evolutionary relationships between distantly related species has difficulties because genetic events such as gene duplication, deletion, recombination, and/or conversion have occurred frequently in these genes. Human MHC class I genes have been conserved only within the primates for up to 46-66 My. Here, we performed comprehensive analysis of the MHC class I genes of the medaka fish, Oryzias latipes, and found that they could be classified into four groups of ancient origin. In phylogenetic analysis using these genes and the classical and nonclassical class I genes of other teleost fishes, three extracellular domains of the class I genes showed quite different evolutionary histories. The α1 domains generated four deeply diverged lineages corresponding to four medaka class I groups with high bootstrap values. These lineages were shared with salmonid and/or other acanthopterygian class I genes, unveiling the orthologous relationships between the classical MHC class I genes of medaka and salmonids, which diverged approximately 260 Ma. This suggested that the lineages must have diverged in the early days of the euteleost evolution and have been maintained for a long time in their genome. In contrast, the α3 domains clustered by species or fish groups, regardless of classical or nonclassical gene types, suggesting that this domain was homogenized in each species during prolonged evolution, possibly retaining the potential for CD8 binding even in the nonclassical genes. On the other hand, the α2 domains formed no apparent clusters with the α1 lineages or with species, suggesting that they were diversified partly by interlocus gene conversion, and that the α1 and α2 domains evolved separately. Such evolutionary mode is characteristic to the teleost MHC class I genes and might have contributed to the long-term conservation of the α1 domain.

MHC evolution in three salmonid species: a comparison between class II alpha and beta genes

The genes of the major histocompatibility complex (MHC) are amongst the most variable in vertebrates and represent some of the best candidates to study processes of adaptive evolution. However, despite the number of studies available, most of the information on the structure and function of these genes come from studies in mammals and birds in which the MHC class I and II genes are tightly linked and class II alpha exhibits low variability in many cases. Teleost fishes are among the most primitive vertebrates with MHC and represent good organisms for the study of MHC evolution because their class I and class II loci are not physically linked, allowing for independent evolution of both classes of genes. We have compared the diversity and molecular mechanisms of evolution of classical MH class II alpha and class II beta loci in farm populations of three salmonid species: Oncorhynchus kisutch, Oncorhynchus mykiss and Salmo salar. We found single classical class II loci and high polymorphism at both class II alpha and beta genes in the three species. Mechanisms of evolution were common for both class II genes, with recombination and point mutation involved in generating diversity and positive selection acting on the peptide-binding residues. These results suggest that the maintenance of variability at the class IIalpha gene could be a mechanism to increase diversity in the MHC class II in salmonids in order to compensate for the expression of one single classical locus and to respond to a wider array of parasites.

Comprehensive analysis of MHC class I genes from the U-, S-, and Z-lineages in Atlantic salmon

BACKGROUND

We have previously sequenced more than 500 kb of the duplicated MHC class I regions in Atlantic salmon. In the IA region we identified the loci for the MHC class I gene Sasa-UBA in addition to a soluble MHC class I molecule, Sasa-ULA. A pseudolocus for Sasa-UCA was identified in the nonclassical IB region. Both regions contained genes for antigen presentation, as wells as orthologues to other genes residing in the human MHC region.

RESULTS

The genomic localisation of two MHC class I lineages (Z and S) has been resolved. 7 BACs were sequenced using a combination of standard Sanger and 454 sequencing. The new sequence data extended the IA region with 150 kb identifying the location of one Z-lineage locus, ZAA. The IB region was extended with 350 kb including three new Z-lineage loci, ZBA, ZCA and ZDA in addition to a UGA locus. An allelic version of the IB region contained a functional UDA locus in addition to the UCA pseudolocus. Additionally a BAC harbouring two MHC class I genes (UHA) was placed on linkage group 14, while a BAC containing the S-lineage locus SAA (previously known as UAA) was placed on LG10. Gene expression studies showed limited expression range for all class I genes with exception of UBA being dominantly expressed in gut, spleen and gills, and ZAA with high expression in blood.

CONCLUSION

Here we describe the genomic organization of MHC class I loci from the U-, Z-, and S-lineages in Atlantic salmon. Nine of the described class I genes are located in the extension of the duplicated IA and IB regions, while three class I genes are found on two separate linkage groups. The gene organization of the two regions indicates that the IB region is evolving at a different pace than the IA region. Expression profiling, polymorphic content, peptide binding properties and phylogenetic relationship show that Atlantic salmon has only one MHC class Ia gene (UBA), in addition to a multitude of nonclassical MHC class I genes from the U-, S- and Z-lineages.

Evolutionary analysis of two classical MHC class I loci of the medaka fish, Oryzias latipes: haplotype-specific genomic diversity, locus-specific polymorphisms, and interlocus homogenization

The major histocompatibility complex (MHC) region of the teleost medaka (Oryzias latipes) contains two classical class I loci, UAA and UBA, whereas most lower vertebrates possess or express a single locus. To elucidate the allelic diversification and evolutionary relationships of these loci, we compared the BAC-based complete genomic sequences of the MHC class I region of three medaka strains and the PCR-based cDNA sequences of two more strains and two wild individuals, representing nine haplotypes. These were derived from two geographically distinct medaka populations isolated for four to five million years. Comparison of the genomic sequences showed a marked diversity in the region encompassing UAA and UBA even between the strains derived from the same population, and also showed an ancient divergence of these loci. cDNA analysis indicated that the peptide-binding domains of both UAA and UBA are highly polymorphic and that most of the polymorphisms were established in a locus-specific manner before the divergence of the two populations. Interallelic recombination between exons 2 and 3 encoding these domains was observed. The second intron of the UAA genes contains a highly conserved region with a palindromic sequence, suggesting that this region contributed to the recombination events. In contrast, the alpha3 domain is extremely homogenized not only within each locus but also between UAA and UBA regardless of populations. Two lineages of the transmembrane and cytoplasmic regions are also shared by UAA and UBA, suggesting that these two loci evolved with intimate genetic interaction through gene conversion or unequal crossing over.

MH class IIalpha polymorphism in local and global adaptation of Arctic charr (Salvelinus alpinus L.)

Arctic charr, a highly plastic salmonid that inhabits the circumpolar region, colonized its current environment after the last glaciation. Recent colonization limits the capacity of many techniques to define and characterize constituent populations. As a novel approach, we used the major histocompatibility (MH) class IIalpha gene polymorphism as a marker that would characterize the genetic divergence of global Arctic charr populations caused by drift and by local adaptation to pathogens. We were able to detect significant isolation of all the lineages previously defined by mitochondrial DNA sequencing and also isolation of some populations within those groups. We found that most of the polymorphism of the class IIalpha gene was distributed globally, which indicates ancestral selection; however, in most cases, distinctive allele frequencies and specific haplotypes distinguished each population suggesting that recent selection has also occurred. Although all studied populations showed similar MH class IIalpha polymorphisms, we also found variation in which particular amino acid positions were polymorphic and which were constant in the different populations studied. This variation provides a greater adaptive capacity for the MH class IIalpha receptors in Arctic charr and is yet another illustration of the extraordinary plasticity of the species.

Multiple expressed MHC class II loci in salmonids; details of one non-classical region in Atlantic salmon (Salmo salar)

Background

In teleosts, the Major Histocompatibility Complex (MHC) class I and class II molecules reside on different linkage groups as opposed to tetrapods and shark, where the class I and class II genes reside in one genomic region. Several teleost MHC class I regions have been sequenced and show varying number of class I genes. Salmonids have one major expressed MHC class I locus (UBA) in addition to varying numbers of non-classical genes. Two other more distant lineages are also identifyed denoted L and ZE. For class II, only one major expressed class II alpha (DAA) and beta (DAB) gene has been identified in salmonids so far.

Results

We sequenced a genomic region of 211 kb encompassing divergent MHC class II alpha (Sasa-DBA) and beta (Sasa-DBB) genes in addition to NRGN, TIPRL, TBCEL and TECTA. The region was not linked to the classical class II genes and had some synteny to genomic regions from other teleosts. Two additional divergent and expressed class II sequences denoted DCA and DDA were also identified in both salmon and trout. Expression patterns and lack of polymorphism make these genes non-classical class II analogues. Sasa-DBB, Sasa-DCA and Sasa-DDA had highest expression levels in liver, hindgut and spleen respectively, suggestive of distinctive functions in these tissues. Phylogenetic studies revealed more yet undescribed divergent expressed MHC class II molecules also in other teleosts.

Conclusion

We have characterised one genomic region containing expressed non-classical MHC class II genes in addition to four other genes not involved in immune function. Salmonids contain at least two expressed MHC class II beta genes and four expressed MHC class II alpha genes with properties suggestive of new functions for MHC class II in vertebrates. Collectively, our data suggest that the class II is worthy of more elaborate studies also in other teleost species.

Genomic organization of duplicated major histocompatibility complex class I regions in Atlantic salmon (Salmo salar)

BACKGROUND

We have previously identified associations between major histocompatibility complex (MHC) class I and resistance towards bacterial and viral pathogens in Atlantic salmon. To evaluate if only MHC or also closely linked genes contributed to the observed resistance we ventured into sequencing of the duplicated MHC class I regions of Atlantic salmon.

RESULTS

Nine BACs covering more than 500 kb of the two duplicated MHC class I regions of Atlantic salmon were sequenced and the gene organizations characterized. Both regions contained the proteasome components PSMB8, PSMB9, PSMB9-like and PSMB10 in addition to the transporter for antigen processing TAP2, as well as genes for KIFC1, ZBTB22, DAXX, TAPBP, BRD2, COL11A2, RXRB and SLC39A7. The IA region contained the recently reported MHC class I Sasa-ULA locus residing approximately 50 kb upstream of the major Sasa-UBA locus. The duplicated class IB region contained an MHC class I locus resembling the rainbow trout UCA locus, but although transcribed it was a pseudogene. No other MHC class I-like genes were detected in the two duplicated regions. Two allelic BACs spanning the UBA locus had 99.2% identity over 125 kb, while the IA region showed 82.5% identity over 136 kb to the IB region. The Atlantic salmon IB region had an insert of 220 kb in comparison to the IA region containing three chitin synthase genes.

CONCLUSION

We have characterized the gene organization of more than 500 kb of the two duplicated MHC class I regions in Atlantic salmon. Although Atlantic salmon and rainbow trout are closely related, the gene organization of their IB region has undergone extensive gene rearrangements. The Atlantic salmon has only one class I UCA pseudogene in the IB region while trout contains the four MHC UCA, UDA, UEA and UFA class I loci. The large differences in gene content and most likely function of the salmon and trout class IB region clearly argues that sequencing of salmon will not necessarily provide information relevant for trout and vice versa.