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

Antigen Presentation and Autophagy in Teleost Adaptive Immunity

Infectious diseases are a burden for aquaculture. Antigen processing and presentation (APP) to the immune effector cells that fight pathogens is key in the adaptive immune response. At the core of the adaptive immunity that appeared in lower vertebrates during evolution are the variable genes encoding the major histocompatibility complex (MHC). MHC class I molecules mainly present peptides processed in the cytosol by the proteasome and transported to the cell surface of all cells through secretory compartments. Professional antigen-presenting cells (pAPC) also express MHC class II molecules, which normally present peptides processed from exogenous antigens through lysosomal pathways. Autophagy is an intracellular self-degradation process that is conserved in all eukaryotes and is induced by starvation to contribute to cellular homeostasis. Self-digestion during autophagy mainly occurs by the fusion of autophagosomes, which engulf portions of cytosol and fuse with lysosomes (macroautophagy) or assisted by chaperones (chaperone-mediated autophagy, CMA) that deliver proteins to lysosomes. Thus, during self-degradation, antigens can be processed to be presented by the MHC to immune effector cells, thus, linking autophagy to APP. This review is focused on the essential components of the APP that are conserved in teleost fish and the increasing evidence related to the modulation of APP and autophagy during pathogen infection.

A Highly Complex, MHC-Linked, 350 Million-Year-Old Shark Nonclassical Class I Lineage

Cartilaginous fish, or Chondrichthyes, are the oldest extant vertebrates to possess the MHC and the Ig superfamily-based Ag receptors, the defining genes of the gnathostome adaptive immune system. In this work, we have identified a novel MHC lineage, UEA, a complex multigene nonclassical class I family found in sharks (division Selachii) but not detected in chimaeras (subclass Holocephali) or rays (division Batoidea). This new lineage is distantly related to the previously reported nonclassical class I lineage UCA, which appears to be present only in dogfish sharks (order Squaliformes). UEA lacks conservation of the nine invariant residues in the peptide (ligand)-binding regions (PBR) that bind to the N and C termini of bound peptide in most vertebrate classical class I proteins, which are replaced by relatively hydrophobic residues compared with the classical UAA. In fact, UEA and UCA proteins have the most hydrophobic-predicted PBR of all identified chondrichthyan class I molecules. UEA genes detected in the whale shark and bamboo shark genome projects are MHC linked. Consistent with UEA comprising a very large gene family, we detected weak expression in different tissues of the nurse shark via Northern blotting and RNA sequencing. UEA genes fall into three sublineages with unique characteristics in the PBR. UEA shares structural and genetic features with certain nonclassical class I genes in other vertebrates, such as the highly complex XNC nonclassical class I genes in Xenopus, and we anticipate that each shark gene, or at least each sublineage, will have a unique function, perhaps in bacterial defense.

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.

Major Histocompatibility Complex in Osteichthyes

Based on analysis of available genome sequences, five gene lineages of MHC class I molecules (MHC I-U, -Z, -S, -L and -P) and one gene lineage of MHC class II molecules (MHC II-D) have been identified in Osteichthyes. In the latter lineage, three MHC II molecule sublineages have been identified (MHC II-A, -B and -E). As regards MHC class I molecules in Osteichthyes, it is important to take note of the fact that the lineages U and Z in MHC I genes have been identified in almost all fish species examined so far. Phylogenetic studies into MHC II molecule genes of sublineages A and B suggest that they may be descended from the genes of the sublineage named A/B that have been identified in spotted gar (Lepisosteus oculatus). The sublineage E genes of MHC II molecules, which represent the group of non-polymorphic genes with poor expression in the tissues connected with the immune system, are present in primitive fish, i.e. in paddlefish, sturgeons and spotted gar (Lepisosteus oculatus), as well as in cyprinids (Cyprinidae), Atlantic salmon (Salmo salar), and rainbow trout (Oncorhynchus mykiss). Full elucidation of the details relating to the organisation and functioning of the particular components of the major histocompatibility complex in Osteichthyes can advance the understanding of the evolution of the MHC molecule genes and the immune mechanism.

NLRC5: new cancer buster?

Recent decades, there is significant progress in understanding the mechanisms of tumor progression and immune evasion. The newly discovered protein NLRC5 is demonstrated to participate in regulating cancer immune escape through enhancing MHC class I genes expression in certain tumors. Nevertheless, increasing evidence has revealed that NLRC5 is up-regulated in some other tumors and promote tumor development and progression. The purpose of this review is to describe the role of NLRC5 in tumors and discuss whether NLRC5 can be a potential target in cancer treatment.

The MHC Class Ia Genes in Chenfu's Treefrog (Zhangixalus chenfui) Evolved via Gene Duplication, Recombination, and Selection

Simple Summary

Amphibians, the first terrestrial vertebrates, provide materials for adaptive evolutionary studies, such as the evolution of the major histocompatibility complex (MHC). To date, various MHC evolutionary mechanisms have been identified in frogs, but more research is needed to determine the evolutionary mechanisms of the frog MHC. The main purpose of this study was to evaluate polymorphisms in the MHC class Ia genes of the Chenfu’s Treefrog. The MHC class Ia genes of the Chenfu’s Treefrog have high polymorphism. The mechanisms responsible for the formation of the polymorphisms include gene duplication, recombination, and selection.

Abstract

The molecular mechanisms underlying the evolution of adaptive immunity-related proteins can be deduced by a thorough examination of the major histocompatibility complex (MHC). Currently, in vertebrates, there is a relatively large amount of research on MHCs in mammals and birds. However, research related to amphibian MHC genes and knowledge about the evolutionary patterns is limited. This study aimed to isolate the MHC class I genes from Chenfu’s Treefrog (Zhangixalus chenfui) and reveal the underlying evolutionary processes. A total of 23 alleles spanning the coding region of MHC class Ia genes were identified in 13 individual samples. Multiple approaches were used to test and identify recombination from the 23 alleles. Amphibian MHC class Ia alleles, from NCBI, were used to construct the phylogenetic relationships in MEGA. Additionally, the partition strategy was adopted to construct phylogenetic relationships using MrBayes and MEGA. The sites of positive selection were identified by FEL, PAML, and MEME. In Chenfu’s Treefrog, we found that: (1) recombination usually takes place between whole exons of MHC class Ia genes; (2) there are at least 3 loci for MHC class Ia, and (3) the diversity of genes in MHC class Ia can be attributed to recombination, gene duplication, and positive selection. We characterized the evolutionary mechanisms underlying MHC class Ia genes in Chenfu’s Treefrog, and in so doing, broadened the knowledge of amphibian MHC systems.

Major histocompatibility complex class I (MHC Iα) of Japanese flounder (Paralichthys olivaceus) plays a critical role in defense against intracellular pathogen infection

The major histocompatibility complex (MHC) is a highly polymorphic region of the vertebrate genome that plays a critical role in initiating immune responses towards invading pathogens. It is well known that MHC I molecules play a central role in the immune response to viruses. However, rare literatures were reported the role of MHC I in the resistance to intracellular bacteria. Sequences of MHC Iα were identified in multiple teleost species, including Japanese flounder (Paralichthys olivaceus), however, the immunological function of MHC Iα remain largely unknown. In this study, we examined the expression profile and biological activity of an MHC Iα homologue, PoMHC Iα, from P. olivaceus. Structural analysis showed that PoMHC Iα possesses conserved structural characteristics of MHC Iα proteins, including MHC_I domain, IGc1 domain, transmembrane region. Expression of PoMHC Iα was upregulated in a time-dependent manner by extracellular and intracellular bacterial pathogens and viral pathogen infection. Different expression patterns were exhibited in response to the infection of different types of microbial pathogens in different immune tissues. Recombinant PoMHC Iα increased the capability of host cells to defense against intracellular pathogen Edwardsiella tarda infection and enhanced the expression of immune related genes. The knockdown of PoMHC Iα attenuated the ability of cells to eliminate E. tarda, which was sustained by the in vivo results that overexpression of PoMHC Iα promoted the host defense against invading E. tarda. Antigen uptake assay indicated PoMHC Iα participated in cells antigen presentation. Collectively, this study is the first report that MHC Iα plays an important role in immune defense against intracellular bacterial pathogen in teleost. Taken together, these findings add new insights into the biological function of teleost MHC Iα and emphasize the importance of MHC I gene products for the control of E. tarda infection.

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.

Comparative MHC nomenclature: report from the ISAG/IUIS-VIC committee 2018

Significant progress has been made over the last decade in defining major histocompatibility complex (MHC) diversity at the nucleotide, allele, haplotype, diplotype, and population levels in many non-human species. Much of this progress has been driven by the increased availability and reduced costs associated with nucleotide sequencing technologies. This report provides an update on the activities of the comparative MHC nomenclature committee which is a standing committee of both the International Society for Animal Genetics (ISAG) and the International Union of Immunological Societies (IUIS) where it operates under the umbrella of the Veterinary Immunology Committee (VIC). A previous report from this committee in 2006 defined the role of the committee in providing guidance in the development of a standardized nomenclature for genes and alleles at MHC loci in non-human species. It described the establishment of the Immuno Polymorphism Database, IPD-MHC, which continues to provide public access to high quality MHC sequence data across a range of species. In this report, guidelines for the continued development of a universal MHC nomenclature framework are described, summarizing the continued development of each species section within the IPD-MHC project.

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.

Whole genome duplications have provided teleosts with many roads to peptide loaded MHC class I molecules

Background

In sharks, chickens, rats, frogs, medaka and zebrafish there is haplotypic variation in MHC class I and closely linked genes involved in antigen processing, peptide translocation and peptide loading. At least in chicken, such MHCIa haplotypes of MHCIa, TAP2 and Tapasin are shown to influence the repertoire of pathogen epitopes being presented to CD8+ T-cells with subsequent effect on cell-mediated immune responses.

Results

Examining MHCI haplotype variation in Atlantic salmon using transcriptome and genome resources we found little evidence for polymorphism in antigen processing genes closely linked to the classical MHCIa genes. Looking at other genes involved in MHCI assembly and antigen processing we found retention of functional gene duplicates originating from the second vertebrate genome duplication event providing cyprinids, salmonids, and neoteleosts with the potential of several different peptide-loading complexes. One of these gene duplications has also been retained in the tetrapod lineage with orthologs in frogs, birds and opossum.

Conclusion

We postulate that the unique salmonid whole genome duplication (SGD) is responsible for eliminating haplotypic content in the paralog MHCIa regions possibly due to frequent recombination and reorganization events at early stages after the SGD. In return, multiple rounds of whole genome duplications has provided Atlantic salmon, other teleosts and even lower vertebrates with alternative peptide loading complexes. How this affects antigen presentation remains to be established.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1138-9) contains supplementary material, which is available to authorized users.

Multiple major histocompatibility complex class I genes in Asian anurans: Ontogeny and phylogeny

Amphibians, as the first terrestrial vertebrates, offer a window into early major histocompatibility complex (MHC) evolution. We characterized the MHC class I of two Korean amphibians, the Asiatic toad (Bufo gargarizans) and the Japanese tree frog (Hyla japonica). We found at least four transcribed MHC class I (MHC I) loci, the highest number confirmed in any anuran to date. Furthermore, we identified MHC I transcripts in terrestrial adults, and possibly in aquatic larvae, of both species. We conducted a phylogenetic analysis based on MHC I sequence data and found that B. gargarizans and H. japonica cluster together in the superfamily Nobleobatrachia. We further identified three supertypes shared by the two species. Our results reveal substantial variation in the number of MHC I loci in anurans and suggest that certain supertypes have particular physiochemical properties that may confer pathogen resistance.

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.

The MHC class I genes of zebrafish

Major histocompatibility complex (MHC) molecules play a central role in the immune response and in the recognition of non-self. Found in all jawed vertebrate species, including zebrafish and other teleosts, MHC genes are considered the most polymorphic of all genes. In this review we focus on the multi-faceted diversity of zebrafish MHC class I genes, which are classified into three sequence lineages: U, Z, and L. We examine the polygenic, polymorphic, and haplotypic diversity of the zebrafish MHC class I genes, discussing known and postulated functional differences between the different class I lineages. In addition, we provide the first comprehensive nomenclature for the L lineage genes in zebrafish, encompassing at least 15 genes, and characterize their sequence properties. Finally, we discuss how recent findings have shed new light on the remarkably diverse MHC loci of this species.

Antiviral functions of CD8(+) cytotoxic T cells in teleost fish

Cytotoxic T-cells (CTLs) play a pivotal role in eliminating viruses in mammalian adaptive immune system. Many recent studies on T-cell immunity of fish have suggested that teleost CTLs are also important for antiviral immunity. Cellular functional studies using clonal ginbuan crucian carp and rainbow trout have provided in vivo and in vitro evidence that in many respects, virus-specific CTLs of fish have functions similar to those of mammalian CTLs. In addition, mRNA expression profiles of CTL-related molecules, such as CD8, TCR and MHC class I, have shown that in a wide range of fish species, CTLs are involved in antiviral adaptive immunity. These findings are a basis to formulate possible vaccination strategies to trigger effective antiviral CTL responses in teleost fish. This review describes recent advances in our understanding of antiviral CTL functions in teleost fish and discusses vaccination strategies for efficiently inducing CTL activities.

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.

Characterization of the Z lineage Major histocompatability complex class I genes in zebrafish

Zebrafish (Danio rerio) are a valuable model for studying immunity, infection, and hematopoietic disease and have recently been employed for transplantation assays. However, the lack of syngeneic zebrafish creates challenges with identifying immune-matched individuals. The MHC class I genes, which mediate allogeneic recognition in mammals, have been grouped into three broad lineages in zebrafish: the classical U genes on chromosome 19, the Z genes which have been reported to map to chromosome 1, and the L genes that map to multiple loci. Transplantations between individual zebrafish that are matched at the U locus fail to consistently engraft suggesting that additional loci contribute to allogeneic recognition. Although two full-length zebrafish Z transcripts have been described, the genomic organization and diversity of these genes have not been reported. Herein we define ten Z genes on chromosomes 1 and 3 and on an unplaced genomic scaffold. We report that neither of the Z transcripts previously described match the current genome assembly and classify these transcripts as additional gene loci. We characterize full-length transcripts for 9 of these 12 genes. We demonstrate a high level of expression variation of the Z genes between individual zebrafish suggestive of haplotypic variation. We report low level sequence variation for individual Z genes between individual zebrafish reflecting a possible nonclassical function, although these molecules may still contribute to allogeneic recognition. Finally, we present a gene nomenclature system for the Z genes consistent with MHC nomenclature in other species and with the zebrafish gene nomenclature guidelines.

Evolution of MHC class I in the order Crocodylia

The major histocompatibility complex (MHC) is a dynamic genomic region with an essential role in the adaptive immunity of jawed vertebrates. The evolution of the MHC has been dominated by gene duplication and gene loss, commonly known as the birth-and-death process. Evolutionary studies of the MHC have mostly focused on model species. However, the investigation of this region in non-avian reptiles is still in its infancy. To provide insights into the evolutionary mechanisms that have shaped the diversity of this region in the Order Crocodylia, we investigated MHC class I exon 3, intron 3, and exon 4 across 20 species of the families Alligatoridae and Crocodilidae. We generated 124 DNA sequences and identified 31 putative functional variants as well as 14 null variants. Phylogenetic analyses revealed three gene groups, all of which were present in Crocodilidae but only one in Alligatoridae. Within these groups, variants generally appear to cluster at the genus or family level rather than in species-specific groups. In addition, we found variation in gene copy number and some indication of interlocus recombination. These results suggest that MHC class I in Crocodylia underwent independent events of gene duplication, particularly in Crocodilidae. These findings enhance our understanding of MHC class I evolution and provide a preliminary framework for comparative studies of other non-avian reptiles as well as diversity assessment within Crocodylia.

Comprehensive analysis of medaka major histocompatibility complex (MHC) class II genes: implications for evolution in teleosts

The major histocompatibility complex (MHC) class II molecules play central roles in adaptive immunity by regulating immune response via the activation of CD4 T cells. The full complement of the MHC class II genes has been elucidated only in mammalian species to date. To understand the evolution of these genes, we performed their first comprehensive analysis in nonmammalian species using a teleost, medaka (Oryzias latipes). Based on a database search, cDNA cloning, and genomic PCR, medaka was shown to possess five pairs of expressed class II genes, comprising one IIA and one IIB gene. Each pair was located on a different chromosome and was not linked to the class I genes. Only one pair showed a high degree of polymorphism and was considered to be classical class II genes, whereas the other four pairs were nonclassical. Phylogenetic analysis of all medaka class II genes and most reported teleost class II genes revealed that the IIA and IIB genes formed separate clades, each containing three well-corresponding lineages. One lineage contained three medaka genes and all known classical class II genes of Ostariophysi and Euteleostei and was presumed to be an original lineage of the teleost MHC class II genes. The other two lineages contained one nonclassical medaka gene each and some Euteleostei genes. These results indicate that multiple lineages of the teleost MHC class II genes have been conserved for hundreds of millions of years and that the tightly linked IIA and IIB genes have undergone concerted evolution.

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.

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.

MHC class I of saltwater crocodiles (Crocodylus porosus): polymorphism and balancing selection

Saltwater crocodiles are in high demand for the production of luxury fashion items. However, their susceptibility to disease incurs substantial losses and it is hoped to be able to genetically select these animals for disease resistance. So far, this has only been enabled by phenotypic selection. Investigating the major histocompatibility complex (MHC) could provide insight into the ability of an individual to respond to pathogens acting as a selective pressure on the host. Here, we assessed genetic diversity and a role of selection in shaping the diversity of MHC class I exon 3 among 42 saltwater crocodiles from nine river basins in the Northern Territory, Australia. We generated 640 sequences using cloning and sequencing methods and identified 43 MHC variants among them. Phylogenetic analyses clustered these variants into two major clades, which may suggest two gene lineages. We found the number of variants within an individual varying between one and seven, indicating that there are at least four gene loci in this species. Selection detection analyses revealed an elevated ratio of nonsynonymous to synonymous substitutions (mean = 1.152 per codon), suggesting balancing selection. Population differentiation analyses revealed that the MHC did not show structuring among the river basins, and there were some shared variants among them. This may be a result of possible gene flow and/or similar selection pressures among populations. These findings provide background knowledge to identify potential MHC markers, which could be used for selecting genetically variable individuals for future disease associations. All MHC class I exon 3 sequences reported in this paper were submitted to the GenBank database with following accession numbers: HQ008785-HQ008789, HQ008791-HQ008798, HQ008808-HQ008815, HQ008824, HQ008826-HQ008830, HQ008835, HQ008839, HQ008842-HQ008850, and JX023536-JX023540.