Interchromosomal duplication of major histocompatibility complex class I regions in rainbow trout (Oncorhynchus mykiss), a species with a presumably recent tetraploid ancestry

An Ancestral Major Histocompatibility Complex Organization in Cartilaginous Fish: Reconstructing MHC Origin and Evolution

Cartilaginous fish (sharks, rays, and chimeras) comprise the oldest living jawed vertebrates with a mammalian-like adaptive immune system based on immunoglobulins (Ig), T-cell receptors (TCRs), and the major histocompatibility complex (MHC). Here, we show that the cartilaginous fish "adaptive MHC" is highly regimented and compact, containing (i) a classical MHC class Ia (MHC-Ia) region containing antigen processing (antigen peptide transporters and immunoproteasome) and presenting (MHC-Ia) genes, (ii) an MHC class II (MHC-II) region (with alpha and beta genes) with linkage to beta-2-microglobulin (β2m) and bromodomain-containing 2, (iii) nonclassical MHC class Ib (MHC-Ib) regions with 450 million-year-old lineages, and (iv) a complement C4 associated with the MHC-Ia region. No MHC-Ib genes were found outside of the elasmobranch MHC. Our data suggest that both MHC-I and MHC-II genes arose after the second round of whole-genome duplication (2R) on a human chromosome (huchr) 6 precursor. Further analysis of MHC paralogous regions across early branching taxa from all jawed vertebrate lineages revealed that Ig/TCR genes likely arose on a precursor of the huchr9/12/14 MHC paralog. The β2m gene is linked to the Ig/TCR genes in some vertebrates suggesting that it was present at 1R, perhaps as the donor of C1 domain to the primordial MHC gene. In sum, extant cartilaginous fish exhibit a conserved and prototypical MHC genomic organization with features found in various vertebrates, reflecting the ancestral arrangement for the jawed vertebrates.

Analysis of microbiota-host communication mediated by butyrate in Atlantic Salmon

Butyrate is a microbiota-produced metabolite, sensed by host short-chain fatty acid receptors FFAR2 (Gpr43), FFAR3 (Gpr41), HCAR2 (Gpr109A), and Histone deacetylase (HDAC) that promotes microbiota-host crosstalk. Butyrate influences energy uptake, developmental and immune response in mammals. This microbial metabolite is produced by around 79 anaerobic genera present in the mammalian gut, yet little is known about the role of butyrate in the host-microbiota interaction in salmonid fish. To further our knowledge of this interaction, we analyzed the intestinal microbiota and genome of Atlantic salmon (Salmo salar), searching for butyrate-producing genera and host butyrate receptors. We identified Firmicutes, Proteobacteria, and Actinobacteria as the main butyrate-producing bacteria in the salmon gut microbiota. In the Atlantic salmon genome, we identified an expansion of genes orthologous to FFAR2 and HCAR2 receptors, and class I and IIa HDACs that are sensitive to butyrate. In addition, we determined the expression levels of orthologous of HCAR2 in the gut, spleen, and head-kidney, and FFAR2 in RTgutGC cells. The effect of butyrate on the Atlantic salmon immune response was evaluated by analyzing the pro and anti-inflammatory cytokines response in vitro in SHK-1 cells by RT-qPCR. Butyrate decreased the expression of the pro-inflammatory cytokine IL-1β and increased anti-inflammatory IL-10 and TGF-β cytokines. Butyrate also reduced the expression of interferon-alpha, Mx, and PKR, and decreased the viral load at a higher concentration (4 mM) in cells treated with this molecule before the infection with Infectious Pancreatic Necrosis Virus (IPNV) by mechanisms independent of FFAR2, FFAR3 and HCAR2 expression that probably inhibit HDAC. Moreover, butyrate modified phosphorylation of cytoplasmic proteins in RTgutGC cells. Our data allow us to infer that Atlantic salmon have the ability to sense butyrate produced by their gut microbiota via different specific targets, through which butyrate modulates the immune response of pro and anti-inflammatory cytokines and the antiviral response.

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.

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.

An Illumina approach to MHC typing of Atlantic salmon

The IPD-MHC Database represents the official repository for non-human major histocompatibility complex (MHC) sequences, overseen and supported by the Comparative MHC Nomenclature Committee, providing access to curated MHC data and associated analysis tools. IPD-MHC gathers allelic MHC class I and class II sequences from classical and non-classical MHC loci from various non-human animals including pets, farmed and experimental model animals. So far, Atlantic salmon and rainbow trout are the only teleost fish species with MHC class I and class II sequences present. For the remaining teleost or ray-finned species, data on alleles originating from given classical locus is scarce hampering their inclusion in the database. However, a fast expansion of sequenced genomes opens for identification of classical loci where high-throughput sequencing (HTS) will enable typing of allelic variants in a variety of new teleost or ray-finned species. HTS also opens for large-scale studies of salmonid MHC diversity challenging the current database nomenclature and analysis tools. Here we establish an Illumina approach to identify allelic MHC diversity in Atlantic salmon, using animals from an endangered wild population, and alter the salmonid MHC nomenclature to accommodate the expected sequence expansions.

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.

Distribution of ancient α1 and α2 domain lineages between two classical MHC class I genes and their alleles in grass carp

Major histocompatibility complex (MHC) class I molecules play a crucial role in the immune response by binding and presenting pathogen-derived peptides to specific CD8⁺ T cells. From cDNA of 20 individuals of wild grass carp (Ctenopharyngodon idellus), we could amplify one or two alleles each of classical MHC class I genes Ctid-UAA and Ctid-UBA. In total, 27 and 22 unique alleles of Ctid-UAA and Ctid-UBA were found. The leader, α1, transmembrane and cytoplasmic regions distinguish between Ctid-UAA and Ctid-UBA, and their encoded α1 domain sequences belong to the ancient lineages α1-V and α1-II, respectively, which separated several hundred million years ago. However, Ctid-UAA and Ctid-UBA share allelic lineage variation in their α2 and α3 sequences, in a pattern suggestive of past interlocus recombination events that transferred α2+α3 fragments. The allelic Ctid-UAA and Ctid-UBA variation involves ancient variation between domain lineages α2-I and α2-II, which in the present study was dated back to before the ancestral separation of teleost fish and spotted gar (> 300 million years ago). This is the first report with compelling evidence that recombination events combining different ancient α1 and α2 domain lineages had a major impact on the allelic variation of two different classical MHC class I genes within the same species.

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.

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.

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.

Chemokine receptors in Atlantic salmon

Teleost sequence data have revealed that many immune genes have evolved differently when compared to other vertebrates. Thus, each gene family needs functional studies to define the biological role of individual members within major species groups. Chemokine receptors, being excellent markers for various leukocyte subpopulations, are one such example where studies are needed to decipher individual gene function. The unique salmonid whole genome duplication that occurred approximately 95 million years ago has provided salmonids with many additional duplicates further adding to the complexity and diversity. Here we have performed a systematic study of these receptors in Atlantic salmon with particular focus on potential inflammatory receptors. Using the preliminary salmon genome data we identified 48 chemokine or chemokine-like receptors including orthologues to the ten receptors previously published in trout. We found expressed support for 40 of the bona fide salmon receptors. Eighteen of the chemokine receptors are duplicated, and when tested against a diploid sister group the majority were shown to be remnants of the 4R whole genome duplication with subsequent high sequence identity. The salmon chemokine receptor repertoire of 40 expressed bona fide genes is comparably larger than that found in humans with 23 receptors. Diversification has been a major driving force for these duplicate genes with the main variability residing in ligand binding and signalling domains.

The genome and linkage map of the northern pike (Esox lucius): conserved synteny revealed between the salmonid sister group and the Neoteleostei

The northern pike is the most frequently studied member of the Esociformes, the closest order to the diverse and economically important Salmoniformes. The ancestor of all salmonids purportedly experienced a whole-genome duplication (WGD) event, making salmonid species ideal for studying the early impacts of genome duplication while complicating their use in wider analyses of teleost evolution. Studies suggest that the Esociformes diverged from the salmonid lineage prior to the WGD, supporting the use of northern pike as a pre-duplication outgroup. Here we present the first genome assembly, reference transcriptome and linkage map for northern pike, and evaluate the suitability of this species to provide a representative pre-duplication genome for future studies of salmonid and teleost evolution. The northern pike genome sequence is composed of 94,267 contigs (N50 = 16,909 bp) contained in 5,688 scaffolds (N50 = 700,535 bp); the total scaffolded genome size is 878 million bases. Multiple lines of evidence suggest that over 96% of the protein-coding genome is present in the genome assembly. The reference transcriptome was constructed from 13 tissues and contains 38,696 transcripts, which are accompanied by normalized expression data in all tissues. Gene-prediction analysis produced a total of 19,601 northern pike-specific gene models. The first-generation linkage map identifies 25 linkage groups, in agreement with northern pike's diploid karyotype of 2N = 50, and facilitates the placement of 46% of assembled bases onto linkage groups. Analyses reveal a high degree of conserved synteny between northern pike and other model teleost genomes. While conservation of gene order is limited to smaller syntenic blocks, the wider conservation of genome organization implies the northern pike exhibits a suitable approximation of a non-duplicated Protacanthopterygiian genome. This dataset will facilitate future studies of esocid biology and empower ongoing examinations of the Atlantic salmon and rainbow trout genomes by facilitating their comparison with other major teleost groups.

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.

Antibodies recognizing both IgM isotypes in Atlantic salmon

Identification and characterization of subpopulations of cells involved in immunological reactions against invading organisms are essential for understanding defense mechanisms against disease. In lower vertebrates like teleost fish, as opposed to mammals, immune cell subsets are still poorly defined, mostly due to the lack of appropriate working tools like antibodies and functional assays. Membrane bound molecules like immunoglobulins (Ig) serve as cell surface markers for specific cell subsets and the identification of cells relies upon the production of specific antibodies towards these molecules. The present study aimed at identifying tools to separate IgM positive (IgM(+)) B cells from IgM negative (IgM(-)) non-B cell populations using flow cytometry. Several monoclonal antibodies (mAbs), and one polyclonal antibody (pAb) to both rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) IgM, either commercially available or locally produced were tested for their recognition of Atlantic salmon IgM(+) cells. Leukocytes were isolated from peripheral blood (PB), spleen (S) and head kidney (HK) and stained with all mAbs and the pAb, to possibly verify the approximate number of IgM(+) cells in the respective tissues in salmon. To our surprise, this seemingly simple task did not reveal similar staining patterns for all antibodies as expected, but rather large differences in the number of positively stained cells were discovered. In short, positively stained cells by each antibody ranged from below 5% to above 80% with similar ratios between the antibodies in each tissue. The three most used mAbs, 4c10, N2 and 1.14; originally produced towards rainbow trout IgM, recognize only a fraction of salmon B cells as previously shown for the 4c10 mAb binding exclusively to the IgM-A isotype. In comparison, our three novel mAbs, IgF1-3, -18 and -19, bind to both IgM-A and -B isotypes as shown using intracellular staining of 293T cells transfected with both IgM-A and -B constructs. Based on binding percentages, one of three commercially available Abs, IgH FITC from Cedarlane, may also identify both isotypes. The three new IgF1-3, -18 and -19 mAbs and potentially IgH FITC from Cedarlane, provide us with great tools enabling complete depletion or enrichment of IgM(+) B cells and/or IgM(-) T cells in Atlantic salmon.

Long-lived dichotomous lineages of the proteasome subunit beta type 8 (PSMB8) gene surviving more than 500 million years as alleles or paralogs

On an evolutionary time scale, polymorphic alleles are believed to have a short life, persisting at most tens of millions of years even under long-term balancing selection. Here, we report highly diverged trans-species dimorphism of the proteasome subunit beta type 8 (PSMB8) gene, which encodes a catalytic subunit of the immunoproteasome responsible for the generation of peptides presented by major histocompatibility complex (MHC) class I molecules, in lower teleosts including Cypriniformes (zebrafish and loach) and Salmoniformes (trout and salmon), whose last common ancestor dates to 300 Ma. Moreover, phylogenetic analyses indicated that these dimorphic alleles share lineages with two shark paralogous genes, suggesting that these two lineages have been maintained for more than 500 My either as alleles or as paralogs, and that conversion between alleles and paralogs has occurred at least once during vertebrate evolution. Two lineages termed PSMB8A and PSMB8F show an A(31)F substitution that would probably affect their cleaving specificity, and whereas the PSMB8A lineage has been retained by all analyzed jawed vertebrates, the PSMB8F lineage has been lost by most jawed vertebrates except for cartilaginous fish and basal teleosts. However, a possible functional equivalent of the PSMB8F lineage has been revived as alleles within the PSMB8A lineage at least twice during vertebrate evolution in the amphibian Xenopus and teleostean Oryzias species. Dynamic evolution of the PSMB8 polymorphism through long-term persistence, loss, and regaining of dimorphism and conversion between alleles and paralogs implies the presence of strong selective pressure for functional polymorphism of this gene.

Isolation and molecular characterization of Rem2 isoforms in the rainbow trout (Oncorhynchus mykiss): Tissue and central nervous system expression

REM2 is a member of the REM, RAD, and GEM/KIR (RGK) subfamily of RAS superfamily proteins and plays an important role in brain development and function. In this study, two Rem2 isoforms were isolated from the rainbow trout (Oncorhynchus mykiss). The two genes, designated O. mykiss rem2a and rem2b, both encode 304 amino acid proteins with 61% and 62% identities to zebrafish (Danio rerio) Rem2, respectively, and each with 43% identity to mammalian (human) REM2. To our knowledge, this is the first incidence of Rem2 isoforms in a species that are the result of gene duplication. Both isoforms possessed similar tissue expression profiles with the highest levels in the brain. The rem2a gene has significantly higher expression levels than rem2b in all tissues assayed except the brain and head kidney. In the central nervous system, both isoforms showed similar expression levels with the highest levels occurring in the olfactory bulb, cerebrum, and midbrain, though rem2a expression is significantly higher in the spinal cord. Based on known functional roles of Rem2 in synapse development and stem cell proliferation, the characterization of Rem2 in rainbow trout could shed light on its role in adult vertebrate neurogenesis and brain regeneration.

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.

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.

Primordial linkage of β2-microglobulin to the MHC

β2-Microglobulin (β2M) is believed to have arisen in a basal jawed vertebrate (gnathostome) and is the essential L chain that associates with most MHC class I molecules. It contains a distinctive molecular structure called a constant-1 Ig superfamily domain, which is shared with other adaptive immune molecules including MHC class I and class II. Despite its structural similarity to class I and class II and its conserved function, β2M is encoded outside the MHC in all examined species from bony fish to mammals, but it is assumed to have translocated from its original location within the MHC early in gnathostome evolution. We screened a nurse shark bacterial artificial chromosome library and isolated clones containing β2M genes. A gene present in the MHC of all other vertebrates (ring3) was found in the bacterial artificial chromosome clone, and the close linkage of ring3 and β2M to MHC class I and class II genes was determined by single-strand conformational polymorphism and allele-specific PCR. This study satisfies the long-held conjecture that β2M was linked to the primordial MHC (Ur MHC); furthermore, the apparent stability of the shark genome may yield other genes predicted to have had a primordial association with the MHC specifically and with immunity in general.

Evolution of duplicated IgH loci in Atlantic salmon, Salmo salar

BACKGROUND

The Atlantic salmon (Salmo salar) immunoglobulin heavy chain (IgH) locus possesses two parallel IgH isoloci (IGH-A and IGH-B), that are related to the genomic duplication event in the family Salmonidae. These duplicated IgH loci in Atlantic salmon provide a unique opportunity to examine the mechanisms of genome diversity and genome evolution of the IgH loci in vertebrates. In this study, we defined the structure of these loci in Atlantic salmon, and sequenced 24 bacterial artificial chromosome (BAC) clones that were assembled into the IGH-A (1.1 Mb) and IGH-B (0.9 Mb) loci. In addition, over 7,000 cDNA clones from the IgH variable (VH) region have been sequenced and analyzed.

RESULTS

The present study shows that the genomic organization of the duplicated IgH loci in Atlantic salmon differs from that in other teleosts and other vertebrates. The loci possess multiple Cτ genes upstream of the Cμ region, with three of the Cτ genes being functional. Moreover, the duplicated loci possess over 300 VH segments which could be classified into 18 families. This is the largest number of VH families currently defined in any vertebrate. There were significant structural differences between the two loci, indicating that both IGH-A and -B loci have evolved independently in the short time after the recent genome duplication approximately 60 mya.

CONCLUSIONS

Our results indicate that the duplication of the IgH loci in Atlantic salmon significantly contributes to the increased diversity of the antibody repertoire, as compared with the single IgH locus in other vertebrates.

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.

Identification, characterization and genetic mapping of TLR7, TLR8a1 and TLR8a2 genes in rainbow trout (Oncorhynchus mykiss)

Induction of the innate immune pathways is critical for early anti-viral defense but there is limited understanding of how teleost fish recognize viral molecules and activate these pathways. In mammals, Toll-like receptors (TLR) 7 and 8 bind single-stranded RNA of viral origin and are activated by synthetic anti-viral imidazoquinoline compounds. Herein, we identify and describe the rainbow trout (Oncorhynchus mykiss) TLR7 and TLR8 gene orthologs and their mRNA expression. Two TLR7/8 loci were identified from a rainbow trout bacterial artificial chromosome (BAC) library using DNA fingerprinting and genetic linkage analyses. Direct sequencing of two representative BACs revealed intact omTLR7 and omTLR8a1 open reading frames (ORFs) located on chromosome 3 and a second locus on chromosome 22 that contains an omTLR8a2 ORF and a putative TLR7 pseudogene. We used the omTLR8a1/2 nomenclature for the two trout TLR8 genes as phylogenetic analysis revealed that they and all the other teleost TLR8 genes sequenced to date are similar to the zebrafish TLR8a, but are distinct from the zebrafish TLR8b. The duplicated trout loci exhibit conserved synteny with other fish genomes extending beyond the tandem of TLR7/8 genes. The trout TLR7 and 8a1/2 genes are composed of a single large exon similar to all other described TLR7/8 genes. The omTLR7 ORF is predicted to encode a 1049 amino acid (aa) protein with 84% similarity to the Fugu TLR7 and a conserved pattern of predicted leucine-rich repeats (LRR). The omTLR8a1 and omTLR8a2 are predicted to encode 1035- and 1034-aa proteins, respectively, and have 86% similarity to each other. omTLR8a1 is likely the ortholog of the only Atlantic salmon TLR8 gene described to date as they have 95% aa sequence similarity. The tissue expression profiles of omTLR7, omTLR8a1 and omTLR8a2 in healthy trout were highest in spleen tissue followed by anterior and then posterior kidney tissues. Rainbow trout anterior kidney leukocytes produced elevated levels of pro-inflammatory and type I interferon cytokines mRNA in response to stimulation with the human TLR7/8 agonist R848 or the TLR3 agonist poly I:C. Only poly I:C-induced IFN2 transcription was significantly suppressed in the presence of chloroquine, a compound known to block endosomal acidification and inhibit endosomal maturation. The effect of chloroquine on R848-induced cytokine expression was equivocal and so it remains questionable whether rainbow trout recognition of R848 requires endosomal maturation. TLR7 and TLR8a1 expression levels in rainbow trout anterior kidney leukocytes were not affected by poly I:C or R848 treatments, but surprisingly, TLR8a2 expression was moderately down-regulated by R848. The down-regulation of omTLR8a2 may imply that this gene has evolved to a new or altered function in rainbow trout, as often occurs when the two duplicated genes remain active.

[Genome similarity of Baikal omul and sig]

Two members of the Baikal sig family, a lake sig (Coregonus lavaretus baicalensis Dybovsky) and omul (C. autumnalis migratorius Georgi), are close relatives that diverged from the same ancestor 10-20 thousand years ago. In this work, we studied genomic polymorphism of these two fish species. The method of subtraction hybridization (SH) did not reveal the presence of extended sequences in the sig genome and their absence in the omul genome. All the fragments found by SH corresponded to polymorphous noncoding genome regions varying in mononucleotide substitutions and short deletions. Many of them are mapped close to genes of the immune system and have regions identical to the Tc-1-like transposons abundant among fish, whose transcription activity may affect the expression of adjacent genes. Thus, we showed for the first time that genetic differences between Baikal sig family members are extremely small and cannot be revealed by the SH method. This is another endorsement of the hypothesis on the close relationship between Baikal sig and omul and their evolutionarily recent divergence from a common ancestor.

A review of the need and possible uses for genetically standardized Atlantic salmon (Salmo salar) in research

Large numbers of Atlantic salmon ( Salmo salar) are used as research animals in basic research and to solve challenges related to the fish-farming industry. Most of this research is performed on farmed animals provided by local breeders or national breeding companies. The genetic constitution of these animals is usually unknown and highly variable. As a result, large numbers of fish are often needed to produce significant results, and results from one study are often impossible to reproduce in another facility. The production of standardized salmon could in many cases reduce the number of animals used in research and at the same time provide more reproducible results. This paper provides an overview of the methods available for the production of standardized Atlantic salmon, and discusses the pros and cons of each technique. The use of zebrafish and other well-defined laboratory fish species as a model for salmon is also discussed. Access to genetically defined fish would greatly benefit the scientific community, in the same way as genetically defined lines of rodents have revolutionized mammalian research.

Dichotomous haplotypic lineages of the immunoproteasome subunit genes, PSMB8 and PSMB10, in the MHC class I region of a Teleost Medaka, Oryzias latipes

Sequence comparison of the medaka, Oryzias latipes, major histocompatibility complex (MHC) class I region between two inbred strains, the HNI (derived from the Northern Population) and the Hd-rR (from the Southern Population), revealed a approximately 100 kb highly divergent segment encompassing two MHC class IA genes, Orla-UAA and Orla-UBA, and two immunoproteasome beta subunit genes, PSMB8 and PSMB10. To elucidate the genetic diversity of this region, we analyzed polymorphisms of the PSMB8 and PSMB10 genes using wild populations of medaka from three genetically different groups: the Northern Population, the Southern Population, and the China-West Korean Population. A total of 1,245 specimens from 10 localities were analyzed, and all the PSMB8 and PSMB10 alleles were classified into the N (fixed in the HNI strain) or the d (fixed in the Hd-rR strain) lineage. Polymerase chain reaction analysis of the region from PSMB8 to PSMB10 indicated that the two allelic lineages of these genes are segregating together constituting dichotomous haplotypic lineages. Both haplotypic lineages were identified in all three groups, although the frequency of d haplotypic lineage (73-100%) was much higher than that of N haplotypic lineage (0-27%) in all analyzed populations. The two allelic lineages of the PSMB8 gene showed curious substitutions at the 31st and 53rd residues of the mature peptide, which are likely involved in formation of the S1 pocket, suggesting that these alleles have a functional difference in cleavage specificity. These results indicate that the two medaka MHC haplotypic lineages encompassing the PSMB8 and PSMB10 genes are maintained in wild populations by a balancing selection.

Suggestive association of major histocompatibility IB genetic markers with resistance to bacterial cold water disease in rainbow trout (Oncorhynchus mykiss)

Genes within the major histocompatibility complex (MHC) are important for both innate and adaptive immune responses in mammals; however, much less is known regarding their contribution in teleost fishes. We examined the involvement of four major histocompatibility (MH) genomic regions in rainbow trout in resistance to the causative agent of bacterial coldwater disease (BCWD), Flavobacterium psychrophilum. Fish from the 2005 NCCCWA brood-year (71 full-sib families) were challenged with F. psychrophilum strain CSF 259-93. The overall mortality rate was 70%, with large variation in mortality between families. Disease resistance was quantified as post-challenge days to death. Phenotypic variation and additive genetic variation were estimated using mixed models of survival analysis. To examine association, eight microsatellite markers were isolated from MH gene-containing BAC clones and mapped onto the rainbow trout genetic linkage map. The parents and grandparents of the 2005 brood-year families were genotyped with these eight markers and another two markers tightly linked to the MH-IB region to assess the extent of linkage disequilibrium (LD) of MH genomic regions MH-IA, MH-IB, TAP1, and MH-II with survival post-challenge. MH-IB and MH-II markers were linked to BCWD survivability when data were analyzed by family. Tests for disease association at the population level substantiated the involvement of MH-IB, but not MH-II, with disease resistance. The impact of selective breeding for disease resistance on MH sequence variation is discussed in the context of aquaculture production.

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.

Bursts and horizontal evolution of DNA transposons in the speciation of pseudotetraploid salmonids

Background

Several genome duplications have occurred in the evolutionary history of teleost fish. In returning to a stable diploid state, the polyploid genome reorganized, and large portions are lost, while the fish lines evolved to numerous species. Large scale transposon movement has been postulated to play an important role in the genome reorganization process. We analyzed the DNA sequence of several large loci in Salmo salar and other species for the presence of DNA transposon families.

Results

We have identified bursts of activity of 14 families of DNA transposons (12 Tc1-like and 2 piggyBac-like families, including 11 novel ones) in genome sequences of Salmo salar. Several of these families have similar sequences in a number of closely and distantly related fish, lamprey, and frog species as well as in the parasite Schistosoma japonicum. Analysis of sequence similarities between copies within the families of these bursts demonstrates several waves of transposition activities coinciding with salmonid species divergence. Tc1-like families show a master gene-like copying process, illustrated by extensive but short burst of copying activity, while the piggyBac-like families show a more random copying pattern. Recent families may include copies with an open reading frame for an active transposase enzyme.

Conclusion

We have identified defined bursts of transposon activity that make use of master-slave and random mechanisms. The bursts occur well after hypothesized polyploidy events and coincide with speciation events. Parasite-mediated lateral transfer of transposons are implicated.

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.

Major histocompatibility polymorphism associated with resistance towards amoebic gill disease in Atlantic salmon (Salmo salar L.)

The association between major histocompatibility (MH) polymorphism and the severity of infection by amoebic gill disease (AGD) was investigated across 30 full sibling families of Atlantic salmon. Individuals were challenged with AGD for 19days and then their severity of infection scored by histopathological examination of the gills. Fish were then genotyped for the MH class I (Sasa-UBA) and MH class II alpha (Sasa-DAA) genes using polymorphic repeats embedded within the 3' untranslated regions of the Sasa-UBA and Sasa-DAA genes. High variation in the severity of infection was observed across the sample material, ranging from 0% to 85% gill filaments infected. In total, seven Sasa-DAA-3UTR and ten Sasa-UBA-3UTR marker alleles were identified across the 30 families. A significant association between the marker allele Sasa-DAA-3UTR 239 and a reduction in AGD severity was detected. There was also a significant association found between AGD severity and the presence of two Sasa-DAA-3UTR genotypes. While the associations between MH allele/genotypes and AGD severity reported herein may be statistically significant, the small sample sizes observed for some alleles and genotypes means these associations should be considered as suggestive and future research is required to verify their biological significance.

A third broad lineage of major histocompatibility complex (MHC) class I in teleost fish; MHC class II linkage and processed genes

Most of the previously studied teleost MHC class I molecules can be classified into two broad lineages: "U" and "Z/ZE." However, database reports on genes in cyprinid and salmonid fishes show that there is a third major lineage, which lacks detailed analysis so far. We designated this lineage "L" because of an intriguing linkage characteristic. Namely, one zebrafish L locus is closely linked with MHC class II loci, despite the extensively documented nonlinkage of teleost class I with class II. The L lineage consists of highly variable, nonclassical MHC class I genes, and has no apparent orthologues outside teleost fishes. Characteristics that distinguish the L lineage from most other MHC class I are (1) absence of two otherwise highly conserved tryptophan residues W51 and W60 in the alpha1 domain, (2) a low GC content of the alpha1 and alpha2 exons, and (3) an HINLTL motif including a possible glycosylation site in the alpha3 domain. In rainbow trout (Oncorhynchus mykiss) we analyzed several intact L genes in detail, including their genomic organization and transcription pattern. The gene Onmy-LAA is quite different from the genes Onmy-LBA, Onmy-LCA, Onmy-LDA, and Onmy-LEA, while the latter four are similar and categorized as "Onmy-LBA-like." Whereas the Onmy-LAA gene is organized like a canonical MHC class I gene, the Onmy-LBA-like genes are processed and lack all introns except intron 1. Onmy-LAA is predominantly expressed in the intestine, while the Onmy-LBA-like transcripts display a rather homogeneous tissue distribution. To our knowledge, this is the first description of an MHC class I lineage with multiple copies of processed genes, which are intact and transcribed. The present study significantly improves the knowledge of MHC class I variation in teleosts.

Contrasting mode of evolution between the MHC class I genomic region and class II region in the three-spined stickleback (Gasterosteus aculeatus L.; Gasterosteidae: Teleostei)

Major histocompatibility complex (MHC) class I molecules display peptides on cell surfaces for subsequent T-cell recognition and are involved in the immune response against intracellular pathogens. In this study, a BAC library was created from a single three-spined stickleback and screened for clones containing MHC class I genes. In a 163.2-kb genomic sequence segment of a single clone, we identified three MHC class I genes in the same transcriptional orientation. Two class I genes are potentially expressed and functional. In one class I gene, the transmembrane region is missing and could therefore present a pseudogene. Alternatively, it presents a functional gene that encodes a soluble MHC class Ib molecule. Despite genomic similarities to the MHC class II region, which is characterized by interlocus recombination, we did not find any evidence for this kind of recombination in the class I genes. It thus seems that interlocus recombination may play a rather minor role in generating class I diversity in stickleback and that the class I region displays a higher genomic stability (i.e., lower local recombination rate). In addition, two non-MHC genes (Oct-2 beta and Na(+),K(+)-ATPasealpha3) have been identified in the analyzed class I region. The Oct-2 beta gene is a transcription factor that is expressed primarily in B lymphocytes, in activated T-cells, and in neuronal cells. The Na(+),K(+)-ATPasealpha3 gene is primarily expressed in the brain and heart and mediates catalytic activities. Both genes are located on the same linkage group together with the MHC class I genes in the zebra fish. In humans, however, homologues of Oct-2 beta and ATPasealpha3 lie outside the MHC region, which indicates that the concentration of immune genes found in mammalian genomes is a derived state.

Cloning and expression analysis of an Atlantic salmon (Salmo salar L.) tapasin gene

Loading of the major histocompatibility complex (MHC) class I molecule with peptide is mediated by the multimeric peptide-loading complex in the ER where the glycoprotein tapasin (TAPBP) is required for stabilization of the complex and for control of peptide loading onto MHC class I. To expand our knowledge on antigen presentation genes in Atlantic salmon, we isolated a full-length salmon tapasin cDNA sequence (Sasa-TAPBP). It encoded a 443 bp amino acid sequence with two N-glycosylation sites, two conserved mammalian tapasin signature motifs, two Ig superfamily (IgSf) domains, a transmembrane (TM) domain and an ER-retention KK motif at the C-terminal end, indicative of a similar function as mammalian tapasins. We analysed the regulation of Sasa-TAPBP under immunostimulatory conditions and found an mRNA-upregulation during early infectious salmon anemia virus (ISAV) infection and poly I:C stimulation in vivo and in vitro, in line with our previous findings for other MHC class I pathway genes.

Evolutionary history of the ABCB2 genomic region in teleosts

Gene duplication, silencing and translocation have all been implicated in shaping the unique genomic architecture of the teleost MH regions. Previously, we demonstrated that trout possess five unlinked regions encoding MH genes. One of these regions harbors ABCB2 which in all other vertebrate classes is found in the MHC class II region. In this study, we sequenced a BAC contig for the trout ABCB2 region. Analysis of this region revealed the presence of genes homologous to those located in the human class II (ABCB2, BRD2, psiDAA), extended class II (RGL2, PHF1, SYGP1) and class III (PBX2, Notch-L) regions. The organization and syntenic relationships of this region were then compared to similar regions in humans, Tetraodon and zebrafish to learn more about the evolutionary history of this region. Our analysis indicates that this region was generated during the teleost-specific duplication event while also providing insight about potential MH paralogous regions in teleosts.

The major histocompatibility complex (Mhc) class IIB region has greater genomic structural flexibility and diversity in the quail than the chicken

Background

The quail and chicken major histocompatibility complex (Mhc) genomic regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated class I, class IIB, natural killer (NK)-receptor-like, lectin-like and BG genes. Therefore, the elucidation of genetic factors that contribute to the greater Mhc diversity in the quail would help to establish it as a model experimental animal in the investigation of avian Mhc associated diseases.

Aims and approaches

The main aim here was to characterize the genetic and genomic features of the transcribed major quail MhcIIB (CojaIIB) region that is located between the Tapasin and BRD2 genes, and to compare our findings to the available information for the chicken MhcIIB (BLB). We used four approaches in the study of the quail MhcIIB region, (1) haplotype analyses with polymorphic loci, (2) cloning and sequencing of the RT-PCR CojaIIB products from individuals with different haplotypes, (3) genomic sequencing of the CojaIIB region from the individuals with the different haplotypes, and (4) phylogenetic and duplication analysis to explain the variability of the region between the quail and the chicken.

Results

Our results show that the Tapasin-BRD2 segment of the quail Mhc is highly variable in length and in gene transcription intensity and content. Haplotypic sequences were found to vary in length between 4 to 11 kb. Tapasin-BRD2 segments contain one or two major transcribed CojaIIBs that were probably generated by segmental duplications involving c-type lectin-like genes and NK receptor-like genes, gene fusions between two CojaIIBs and transpositions between the major and minor CojaIIB segments. The relative evolutionary speed for generating the MhcIIBs genomic structures from the ancestral BLB2 was estimated to be two times faster in the quail than in the chicken after their separation from a common ancestor. Four types of genomic rearrangement elements (GRE), composed of simple tandem repeats (STR), were identified in the MhcIIB genomic segment located between the Tapasin-BRD2 genes. The GREs have many more STR numbers in the quail than in the chicken that displays strong linkage disequilibrium.

Conclusion

This study suggests that the Mhc classIIB region has a flexible genomic structure generated by rearrangement elements and rapid SNP accumulation probably as a consequence of the quail adapting to environmental conditions and pathogens during its migratory history after its divergence from the chicken.

A polymerase chain reaction screening method for rapid detection of microsatellites in bacterial artificial chromosomes

Standard protocols aimed at identifying subclones of interest from bacterial artificial chromosomes (BACs) include the use of hybridization methods that are time consuming and often require the use of radioactive isotopes. Through our efforts to identify microsatellites in BACs from rainbow trout (Oncorhynchus mykiss) we have developed a nonradioactive polymerase chain reaction (PCR)-based screening technique to select microsatellites containing subclones for marker development. Two BACs were subcloned and screened by PCR using a vector-specific primer and a mix of microsatellite repeat primers. The subclones were then sequenced to evaluate the efficiency of the PCR screening method. Correlation between positive PCR amplification and presence of microsatellites varied between the two BACs (21.9% and 71.4%), but still a sufficient number of subclones were identified to enable design and optimization of microsatellite markers.

The salmonid MHC class I: more ancient loci uncovered

An unprecedented level of sequence diversity has been maintained in the salmonid major histocompatibility complex (MHC) class I UBA gene, with between lineage AA sequence identities as low as 34%. The derivation of deep allelic lineages may have occurred through interlocus exon shuffling or convergence of ancient loci with the UBA locus, but until recently, no such ancient loci were uncovered. Herein, we document the existence of eight additional MHC class I loci in salmon (UCA, UDA, UEA, UFA, UGA, UHA, ULA, and ZE), six of which share exon 2 and 3 lineages with UBA, and three of which have not been described elsewhere. Half of the UBA exon 2 lineages and all UBA exon 3 lineages are shared with other loci. Two loci, UGA and UEA, share only a single exon lineage with UBA, likely generated through exon shuffling. Based on sequence homologies, we hypothesize that most exchanges and duplications occurred before or during tetraploidization (50 to 100 Ma). Novel loci that share no relationship with other salmonid loci are also identified (UHA and ZE). Each locus is evaluated for its potential to function as a class Ia gene based on gene expression, conserved residues and polymorphism. UBA is the only locus that can indisputably be classified as a class Ia gene, although three of the eight loci (ZE, UCA, and ULA) conform in three out of four measures. We hypothesize that these additional loci are in varying states of degradation to class Ib genes.

Polymorphism of two very similar MHC class Ib loci in rainbow trout (Oncorhynchus mykiss)

As part of an ongoing elucidation of rainbow trout major histocompatibility complex (MHC) class I, the polymorphism of two MHC class Ib loci was analyzed. These loci, Onmy-UCA and Onmy-UDA, are situated head-to-tail and share more than 89% nucleotide identity in their open reading frames. They share 80% identity with some trout Ia alleles. The deduced amino acid sequences suggest that the UCA and UDA molecules are transported to endosomal compartments and may bind peptides in their binding groove. Our survey revealed seven UCA and eight UDA alleles. Similarity indices overlap when comparing within and between UCA and UDA alleles and some cross-locus motif variation is observed. In most trout both UCA and UDA transcripts were found. However, there probably is functional redundancy, because some trout lacked transcription of one of the two loci. Furthermore, for some UCA and UDA alleles, splicing deficiencies, early stop codons, and upstream start codons were found, which may interfere with efficient protein expression. The present study is the first extensive report on MHC class Ib polymorphism assigned to locus in ectotherm species.

Ancestral Organization of the MHC Revealed in the Amphibian Xenopus

With the advent of the Xenopus tropicalis genome project, we analyzed scaffolds containing MHC genes. On eight scaffolds encompassing 3.65 Mbp, 122 MHC genes were found of which 110 genes were annotated. Expressed sequence tag database screening showed that most of these genes are expressed. In the extended class II and class III regions the genomic organization, excluding several block inversions, is remarkably similar to that of the human MHC. Genes in the human extended class I region are also well conserved in Xenopus, excluding the class I genes themselves. As expected from previous work on the Xenopus MHC, the single classical class I gene is tightly linked to immunoproteasome and transporter genes, defining the true class I region, present in all nonmammalian jawed vertebrates studied to date. Surprisingly, the immunoproteasome gene PSMB10 is found in the class III region rather than in the class I region, likely reflecting the ancestral condition. Xenopus DMalpha, DMbeta, and C2 genes were identified, which are not present or not clearly identifiable in the genomes of any teleosts. Of great interest are novel V-type Ig superfamily (Igsf) genes in the class III region, some of which have inhibitory motifs (ITIM) in their cytoplasmic domains. Our analysis indicates that the vertebrate MHC experienced a vigorous rearrangement in the bony fish and bird lineages, and a translocation and expansion of the class I genes in the mammalian lineage. Thus, the amphibian MHC is the most evolutionary conserved MHC so far analyzed.

Identification and regulatory analysis of rainbow trout tapasin and tapasin-related genes

Tapasin (TAPBP) is a key member of MHC class Ia antigen-loading complexes, bridging the class Ia molecule to the transporter associated with antigen presentation (TAP). As part of an ongoing study of MHC genomics in rainbow trout, we have identified two rainbow trout TAPBP genes (Onmy-TAPBP.a and .b) and a similar but distinct TAPBP-related gene (Onmy-TAPBP-R) that had previously only been described in mammals. Physical and genetic mapping indicate that Onmy-TAPBP.a is on chromosome 18 in the MHC class Ia region and that Onmy-TAPBP.b resides on chromosome 14 in the MHC class Ib region. There are also at least two copies of TAPBP-R, Onmy-TAPBP-R.a and Onmy-TAPBP-R.b, located on chromosomes 2 and 3, respectively. Due to the central role of TAPBP expression during acute viral infection, we have characterized the transcriptional profile and regulatory regions for both Onmy-TAPBP and Onmy-TAPBP-R. Transcription of both genes increased during acute infection with infectious hematapoeitic necrosis virus (IHNV) in a fashion indicative of interferon-mediated regulation. Promoter-reporter assays in STE-137 cells demonstrate that the trout TAPBP and TAPBP-R promoters respond to interferon regulatory factors, Onmy-IRF1 and Onmy-IRF2. Overall, TAPBP is expressed at higher levels than TAPBP-R in naïve tissues and TAPBP transcription is more responsive to viral infection and IRF1 and 2 binding.