New beta chain-encoding Mhc class II genes in the carp

Fate of MHCII in salmonids following 4WGD

Major histocompatibility complex (MHC) genes are key players in the adaptive immunity providing a defense against invading pathogens. Although the basic structures are similar when comparing mammalian and teleost MHC class II (MHCII) molecules, there are also clear-cut differences. Based on structural requirements, the teleosts non-classical MHCII molecules do not comply with a function similar to the human HLA-DM and HLA-DO, i.e., assisting in peptide loading and editing of classical MHCII molecules. We have previously studied the evolution of teleost class II genes identifying various lineages and tracing their phylogenetic occurrence back to ancient ray-finned fishes. We found no syntenic MHCII regions shared between cyprinids, salmonids, and neoteleosts, suggesting regional instabilities. Salmonids have experienced a unique whole genome duplication 94 million years ago, providing them with the opportunity to experiment with gene duplicates. Many salmonid genomes have recently become available, and here we set out to investigate how MHCII has evolved in salmonids using Northern pike as a diploid sister phyla, that split from the salmonid lineage prior to the fourth whole genome duplication (4WGD) event. We identified 120 MHCII genes in pike and salmonids, ranging from 11 to 20 genes per species analyzed where DB-group genes had the most expansions. Comparing the MHC of Northern pike with that of Atlantic salmon and other salmonids species provides a tale of gene loss, translocations, and genome rearrangements.

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.

MH II-DAB gene expression in grass carp Ctenopharyngodon idella (Valenciennes) after infection with the ciliate parasite, Ichthyophthirius multifiliis

The grass carp, Ctenopharyngodon idella (Valenciennes), is one of the most extensively aquacultured freshwater fish in China. However, because of the lack of effective control measures and the high-density culture environment, considerable economic losses are caused by infection of C. idella with the parasitic ciliate, Ichthyophthirius multifiliis. The major histocompatibility (MH) DAB gene belongs to antigen-presented genes in the class II genomic region, which is associated with parasite resistance. To understand the relationship of the DAB gene with I. multifiliis infection in grass carp, the expression profiles of MH II-DAB were studied in tissues using real-time quantitative polymerase chain reaction. The results showed that expression of the MH II-DAB gene was up-regulated in head kidney after I. multifiliis infection, and the expression peak appeared earlier in the study (case) group than in the control group. The obvious up-regulation peak of MH II-DAB gene was found at days 2 and 4 in skin; at 12 h to day 4 in spleen; at 12 h and days 1 and 6 in gill; and at day 10 in blood, whereas the MH II-DAB gene was down-regulated in liver and intestines after I. multifiliis infection. These results have implications for better understanding C. idella resistance to I. multifiliis infection.

Isolation and characterization of major histocompatibility class IIβ genes in an endangered North American cyprinid fish, the Rio Grande silvery minnow (Hybognathus amarus)

The major histocompatibility complex (MHC) is a critical component of the adaptive immune response in vertebrates. Due to the role that MHC plays in immunity, absence of variation within these genes may cause species to be vulnerable to emerging diseases. The freshwater fish family Cyprinidae comprises the most diverse and species-rich group of freshwater fish in the world, but some are imperiled. Despite considerable species richness and the long evolutionary history of the family, there are very few reports of MHC sequences (apart from a few model species), and no sequences are reported from endemic North American cyprinids (subfamily Leuciscinae). Here we isolate and characterize the MH Class II beta genes from complementary DNA and genomic DNA of the non-model, endangered Rio Grande silvery minnow (Hybognathus amarus), a North American cyprinid. Phylogenetic reconstruction revealed two groups of divergent MH alleles that are paralogous to previously described loci found in deeply divergent cyprinid taxa including common carp, zebrafish, African large barb and bream. Both groups of alleles were under the influence of diversifying selection yet not all individuals had alleles belonging to both allelic groups. We concluded that the general organization and pattern of variation of MH class II genes in Rio Grande silvery minnow is similar to that identified in other cyprinid fishes studied to date, despite distant evolutionary relationships and evidence of a severe genetic bottleneck.

Molecular characterization and expression analysis of MHC class II alpha and beta genes in large yellow croaker (Pseudosciaena crocea)

MHC class II molecules play an important role in the activation of CD4(+) T cells, which are the central orchestrating cells of an immune response. Here, we report the cloning of MHC class II alpha and beta cDNAs from large yellow croaker (Pscr-DAAs and Pscr-DAB) by expressed sequence tags analysis and RACE-PCR techniques. Three different class II alpha and two class II beta sequences were obtained from spleens of two individual fish. Each of the three class II alpha sequences encodes a polypeptide of 239 amino acids while the two class II beta cDNA sequences encode for a protein of 249 aa. All the characteristic features of MHC class II chain structure could be identified in the deduced proteins of three class II alpha and two class II beta sequences, including the leader peptide, alpha1/beta1 and alpha2/beta2 domains, connecting peptide and transmembrane and cytoplasmic regions, as well as conserved cysteines and N-glycosylation site. RT-PCR analysis showed that MHC class II alpha and beta mRNAs were broadly expressed in various tissues examined, although at different levels. Upon stimulation with inactivated trivalent bacterial vaccine or polyinosinic polycytidylic acid (poly(I:C)), the expression levels of both alpha and beta genes were obviously up-regulated in intestine, kidney and spleen. Real-time PCR analysis demonstrated that the expression levels of class II alpha and beta were quickly up-regulated in spleen, kidney, and intestine at 12 h after induction with poly(I:C), while their expression levels significantly increased at 48 h upon immunization with bacterial vaccine, indicating that the up-regulation of both class II alpha and beta expression was induced by bacterial vaccine or poly(I:C) at the early phase of induction, and that class II alpha and beta transcripts were quicker up-regulated by poly I:C than by bacterial vaccine.

Evolution and trans-species polymorphism of MHC class IIbeta genes in cyprinid fish

The polymorphism of DAB genes encoding MHC IIbeta was investigated in 11 cyprinid species from central Europe. The species belonged to four subfamilies: Cyprininae, Tincinae, Gobioninae and Leuciscinae. Two paralogous groups of sequences, DAB1 and DAB3, were recognised according to the similarity of their nucleotide and amino-acid sequences and from phylogenetic analyses using either partial exon 2 or partial exon 3 sequences. A high allelic variability among species was found for exon 2, indicating extensive MHC polymorphism. Time divergence estimation supports the separation of DAB1 and DAB3 groups predating the separation into fish subfamilies, and a cyprinid origin of the DAB genes. Phylogenetic trees using exon 2 support the hypothesis of trans-species polymorphism, which appears to be limited to the subfamily level, i.e. the presence of sequences from different species in the same allelic group was more often recognised within subfamilies Cyprininae and Leuciscinae than between them. Phylogenetic trees using exon 3 reflect the phylogenetic patterns previously found for Cyprinidae systematics. Specific nucleotides and amino-acids in exon 3 that separate both subfamilies, as well as the species within the Cyprininae subfamily were observed. A lack of segregation in leuciscin species was recognised and the alleles of different leuciscin species tend to share similar motifs in exon 3. This could be explained by the ancient and complicated dispersion history of Cyprininae and the radiation of Leuciscinae. The effects of selective pressures were investigated: (1) within species, (2) among lineages, and (3) among sites. From intraspecific analyses, exon 2 sequences were identified as the targets of diversifying selection, whilst the evolution of exon 3 seems to be under the influence of purifying selection. The analyses among lineages indicate positive selection in many branches when using exon 2, therefore confirming trans-species polymorphism, whilst the DAB lineages of exon 3 are potentially submitted to purifying selection to some extent. Moreover, our results suggest the secondary acquisition of function of DAB1 group after duplication. The analyses among sites reveal that exon 2 exhibits sites under positive selection mostly corresponding to the putative PBR sites involved in the alpha-helix structure of the protein.

Coding sequences of the MHC II beta chain of homozygous rainbow trout (Oncorhynchus mykiss)

Six lines of homozygous rainbow trout (Oncorhynchus mikiss) from different genetic and geographical backgrounds have been produced as aquatic models for biomedical research by the chromosome set manipulation techniques of androgenesis and gynogenesis. Messenger RNA from spleens was extracted. and the MHC II B cDNA sequences, amplified by RT PCR, were cloned into plasmids. Sequences of the MHC II beta2 domains were highly conserved between the different plasmids from the same and different lines of trout. Most of the variability among sequences was found in the amino terminal half of the beta1 domain, which corresponds with the peptide binding region of the MHC II molecule. This diversity suggests that the different lines of trout may exhibit differences in immune response. Rainbow trout MHC II B sequences were similar to the MHC II B sequences of the Pacific salmon (O. gorbuscha, O. tshawytscha, O. nerka, O. miasou, O. kisutch). Southern blot analysis performed on the restricted DNA of the OSU and Hot Creek trout, and the doubled haploid progeny produced by androgenesis from OSU x Hot Creek hybrids indicates that two distinct genes encode the MHC II B sequences and that these genes are unlinked.

MHC class II B genes in the channel catfish (Ictalurus punctatus)

Two different cDNA sequences for major histocompatibility complex (MHC) class II beta chains from the channel catfish (Ictalurus punctatus) have been identified. Homology between these sequences and those previously identified as MHC class II B genes in other teleosts suggests they represent alleles of the DAB locus. The inferred amino acid sequences show strong evidence for a functional polypeptide chain with a peptide binding region. Southern blot analysis reveals polymorphism in the MHC class II B gene(s) of the channel catfish and suggests the presence of two to four genes.

Evidence for multiple distinct major histocompatibility complex class I lineages in teleostean fish

In the context of studies on the expression of MhcCyca-Z sequences of the common carp, PCR amplifications of exon 4 were performed on cDNA obtained from pooled thymi of 20 carp F1 individuals. Five recombinant clones (Cyca-TC3, -TC13, -TC15, -TC17 and -TC18) were found to be 96% similar to the exon 4 region of Cyca-ZA1. Each of the five sequences was unique, and differed in a few positions in both the nucleotide and the derived amino acid sequences from any of the Cyca-Z sequences known to date. These data suggest that multiple Z genes per locus are present in the carp, which are transcribed in the thymus. In the course of analysing the amplified Cyca-Z sequences, serendipity yielded a clone, Cyca-TC16, containing a class I-like sequence substantially different from any other carp class I sequence. The predicted amino acid sequence of Cyca-TC16 was most similar to the class I genes (Lach-U) from the coelacanth (42-46% amino acid identity). Cyca-TC16 contains three conserved beta 2-microglobulin contact residues, and the secondary structure was predicted by computer algorithms to be similar to that of the alpha 3 domain of HLA-A2. Phylogenetic analysis shows that carp class I sequences reside in four distinct clusters: (i) Cyca-Z, Cyca-TC3, -TC13, -TC15, -TC17 and -TC18 together with Caau-Z from ginbuna crucian carp, (ii) Cyca-U with Bree-U (zebrafish) and Sasa-p30 (Atlantic salmon), (iii) Cyca-TC16 with Lach-U (coelacanth), and (iv) Cyca-C4.

Trans-species polymorphism of class II Mhc loci in danio fishes

A characteristic feature of the major histocompatibility complex (Mhc) polymorphism in mammals is the existence of allelic lineages shared by related species. This trans-species polymorphism has thus far been documented only in primates, rodents, and artiodactyls. In this communication we provide evidence that it also exists in cyprinid (bony) fishes at the class II A and B loci coding for the alpha and beta polypeptide chains of the class II alpha:beta heterodimers. The study has focused on three species of the family Cyprinidae, subfamily Rasborinae: the zebrafish (Danio rerio), the giant danio (D. malabaricus), and the pearl danio (D. albolineatus). The polymerase chain reaction was used to amplify and then sequence intron 1 and exon 2 of the class II B loci and exon 2 of the class II A loci in these species. Phylogenetic analysis of the sequences revealed the existence of allelic lineages whose divergence predates the divergence of the three species at both the A and B loci. The lineages at the B locus in particular are separated by large genetic distances. The polymorphism is concentrated in the peptide-binding region sites and is apparently maintained by balancing selection. Sharing of this unique Mhc feature by both bony fishes and mammals suggests that the main function of the Mhc (presentation of peptides to T lymphocytes) has not changed during the last 400 million years of its evolution.

Detection of MHC class II transcripts in lymphoid tissues of the common carp (Cyprinus carpio L.)

In all vertebrates studied to date, the expression of MHC class II genes is known to be restricted to a limited number of tissues and cell types. In order to have a better understanding of the function of the equivalent genes in teleost fish, the distribution of MHC class II beta transcripts (Cyca-DAB) in the common carp (Cyprinus carpio L.) was investigated. RNA was isolated from tissues and leucocytes, cDNA was produced, and amplification of the Cyca-DAB genes was carried out by PCR. Of the organs with known immunological function, the highest level of Cyca-DAB transcription was found in the thymus. Despite their expected different cellular organization, total blood, head kidney, spleen and the second segment of the gut had similar Cyca-DAB expression levels. No class II transcripts were detected in the skeletal muscle. The studies carried out with leucocytes isolated from the lymphoid tissues point to a direct correlation between the levels of expression and the amounts of surface immunoglobulin positive (sIg+) cells present in the different cell fractions. However, thymus leucocytes did not follow this correlation since the highest level of class II expression was found in a thymocyte fraction that contained very low numbers of Ig+ cells. In PBL the Ig+ cells were highly positive whereas the Ig- were weakly positive. Adherent leucocytes shown to be class II positive, although adherent cells from PBL show a lower level of expression compared to those from the spleen and head kidney.

Investigations into the ubiquitous nature of high or low immune responsiveness after divergent selection for antibody production in common carp (Cyprinus carpio L.)

This paper reports on the selection of individual carp with a high or low antibody response, in combination with reproduction by gynogenesis, in order to develop well-characterised inbred carp lines consisting of practically unlimited numbers of carp with the same genotype. Two homozygous progenies, previously characterised as having a high or low immune response to dinitrophenyl keyhole limpet haemocyanin (DNP-KLH), were immunised with either a T-dependent (DNP-human serum albumin (DNP-HSA)) or T-independent (trinitrophenyl lipopolysaccharide (TNP-LPS)) hapten-carrier complex. In comparison with the antibody response after DNP-KLH immunisation, the response to DNP-HSA was observed to be highly variable and did not differ between the divergently selected progenies. This suggests that the divergent selection for antibody production to DNP-KLH has been carrier-specific. Immunisation with T-independent TNP-LPS induced a very rapid response which differed between the high and low responders, and likely measured changes in the DNP-specific precursor pool of B cells caused by the selection. A number of selected individuals with a high immune response to DNP-KLH were infected with Trypanoplasma borreli, a haemoflagellate parasite of carp, to examine a possible relationship between the increase in immune responsiveness and disease resistance, but no change could be detected. However, individual homozygous carp were able to escape inbreeding depression and survive the infection. Such carp would be likely candidates for gynogenetic reproduction to obtain viable inbred carp lines.

Exon-intron organization of Xenopus MHC class II beta chain genes

The amphibian Xenopus laevis is the most primitive vertebrate in which the major histocompatibility complex (MHC) has been defined at the biochemical, functional, and molecular genetic levels. We previously described the isolation and characterization of cDNA clones encoding X. laevis MHC class II beta chains. In the present study, genomic clones encoding class II beta chains were isolated from X. laevis homozygous for the MHC f haplotype. Three class II beta chain genes, designated Xela-DAB, Xela-DBB, and Xela-DCB, were identified. Sequence analysis of these genes showed that Xela-DBB and Xela-DCB correspond to the previously characterized cDNA clones F3 and F8, respectively, whereas Xela-DAB encodes a third, hitherto unidentified class II beta chain of the MHC f haplotype. As a representative of X. laevis class II beta chain genes, the Xela-DAB gene underwent detailed structural analysis. In addition, the nucleotide sequence of Xela-DABf cDNA clones was determined. The Xela-DAB gene is made up of at least six exons, with an exon-intron organization similar to that of a typical mammalian class II beta chain gene. The 5'-flanking region of the Xela-DAB gene contains transcriptional control elements known as X1, X2, and Y, but lacks typical TATA or CCAAT boxes. A notable feature of the X. laevis class II beta chain genes is that the sizes of the introns are larger than those of their mammalian counterparts. As assessed by northern blot analysis, the three class II beta chain genes had similar expression patterns, with the highest level of transcription detected in the intestine. Identification of the Xela-DAB, -DBB, and -DCB genes is consistent with our previous observations, which suggested that the MHC of the tetraploid frog X. laevis is diploidized at the genomic level and contains three class II beta chain genes per haplotype that cross-hybridize to one another under reduced stringency conditions.

Fish major histocompatibility complex genes: an expansion

The advent of polymerase chain reaction technology has provoked a large amount of progress in the field of fish major histocompatibility complex (MHC) research. Many new teleost sequences have been reported in the last four years, including representatives of all classes of MHC genes. While the intron-exon structure of teleost MHC genes is now becoming clear, the organisation of the genes within the teleost MHC is still unclear. The sequences reported to date have been used for phylogenetic analysis and, due to their evolutionary position, are discussed in relation to hypotheses regarding the origin of the MHC. Teleost MHC gene sequences are also examined to see if conserved features of the both the nucleotide and amino acid sequences of higher vertebrate MHC genes are present. Differences in these features will reflect functional differences between teleost and mammalian MHC genes and may also have evolutionary implications.

Identification of major histocompatibility complex genes in the guppy, Poecilia reticulata

The guppy, Poecilia reticulata, a teleostean fish of the order Cyprinodontiformes, has been used extensively in studies of host-parasite interactions, courtship behavior, and mating preference, as well as in ecological and evolutionary genetics. A related species was among the first poikilotherm vertebrates to be used in the study of histocompatibility genes. All these studies could benefit from the identification and characterization of the guppy major histocompatibility complex (Mhc) genes. Here, both class I and class II genes of the guppy are described. The number of expressed loci, as determined by representation of clones in a cDNA library, sequencing, and Southern blot analysis, may be low in both Mhc classes: combined evidence suggests that there may be one expressed class II locus only and one or two expressed class I loci. The variability of aquaristic guppy stocks is very low: only three and two genes have been detected at the class I and class II loci, respectively, in the stocks examined. This genetic paucity is most likely the consequence of breeding practices employed by aquarists and commercial establishments. Limited sampling of wild guppy populations revealed extensive Mhc polymorphism at loci of both classes in nature. Comparison of guppy Mhc sequences with those of other vertebrates has revealed the existence of a set of insertions/deletions which can be used as characters in cladistic analysis to infer phylogenetic relationships among vertebrate taxa and the Mhc genes themselves. These indels are particularly frequent in the regions coding for the loops of alpha 1 and alpha 2 domains of class I proteins.

Divergent selection for antibody production in common carp (Cyprinus carpio L.) using gynogenesis

A base population (n = 101) of carp, consisting of a single hybrid cross, was immunized with the hapten-carrier complex DNP-KLH, to perform a divergent selection for antibody response. Measurement of the DNP-specific antibody response at 12 and 21 days postimmunization, allowed the classification of a low number of individual carp as early/high (10%) or late/low (13%) responders. Three individuals defined as early/high and three defined as late/low responding, were gynogenetically reproduced to obtain corresponding homozygous progenies within one generation only. Upon immunization with DNP-KLH, the antibody response was found to be significantly higher in the early/high responder homozygous offspring. Although the homozygosity of the offspring apparently caused a (s)lower antibody response (compared with the base population), the differences between the high and low responder offspring to indicate a genetic influence on the antibody response. The realized heritability (h2) for antibody production was estimated at 0.37 +/- 0.36. The present study provides the basis for a divergent selection of homozygous inbred carp lines with a genetically controlled difference in antibody response. These inbred lines will allow us to investigate relationship(s) between immune responsiveness and resistance to infectious diseases in fish.

Variability in an MHCMosa class II beta chain-encoding gene in striped bass (Morone saxatilis)

The major histocompatibility complex (MHC) class II B locus of the striped bass (Morone saxatilis) was found to contain multiple forms of the class II B gene. Seven complete MHC class II B cDNA clones were isolated and sequenced, identifying five unique allelic forms of a MHC class II B gene. Among three specimens, each representing a geographically distinct population (Chesapeake Bay, MD; Roanoke River, NC; and Santee-Cooper Reservoir, SC) extensive variability was detected in the beta 1 encoding domain, which corresponds with the functional peptide-binding region (PBR) of known MHC class II molecules. The location of variable amino acid residues in the beta 1 domains corresponds with polymorphic sites observed in other teleosts and higher vertebrates. The amino acid translated beta 2 domain encoding regions, transmembrane regions, and cytoplasmic regions of the five clones correlated well with those of known vertebrate MHC class II proteins. Seventy-one percent of the variability found within the presumed PBR encoded at the MHCMosa class II B locus corresponded with that of the PBR of a human MHC class II B gene. Overall, the Mosa sequences showed greatest similarity to the MHC class II B genes of cichlid fishes, as expected from phylogenetic relationships.

A study of variability in the MHC class II beta 1 and class I alpha 2 domain exons of Atlantic salmon, Salmo salar L

Variability in the most extracellular exons of Atlantic salmon MHC-Sasa class I and class II was studied by polymerase chain reaction (PCR) amplification followed by sequencing. The domains studied were class I alpha 2 and class II beta 1. The material used was genomic DNA of fish, mainly derived from the major Norwegian breeding pool, supplemented with some material from a minor breeder and a local river strain. The analysis revealed extensive variation, most individuals being heterozygous with at least two variants.

Composite origin of major histocompatibility complex genes

Major histocompatibility complex (MHC) genes have now been cloned from representatives of all vertebrate classes except Agnatha. The recent accumulation of sequence data has given great insight into the course of evolution of these genes. Although the primary structure of the MHC genes varies greatly from class to class and also within the individual classes, the general features of the tertiary and quaternary structure have been conserved remarkably well during more than 400 million years of evolution. The ancestral MHC genes may have been assembled from at least three structural elements derived from different gene families. Class II MHC genes appear to have been assembled first, and then to have given rise to class I genes.

Exon-intron organization of fish major histocompatibility complex class II B genes

Major histocompatibility complex (Mhc) molecules bind self and foreign peptides and present them to lymphocytes for recognition. Activation of lymphocytes by Mhc-bound foreign peptides leads to specific immune response against parasites. The Mhc genes have been studied extensively in mammals and birds but much less in other vertebrate classes. In this communication we provide the first description of the exon-intron organization of class II beta-chain-encoding genes from the teleost fish Aulonocara hansbaenschi, family Cichlidae. Each of the genes consists of six exons, E1 through E6, encoding the leader peptide (E1), beta 1 domain (E1+E2), beta 2 domain (E3+E4), connecting peptide (E5), transmembrane region (E5), cytoplasmic domain (E5+E6), and the 3' untranslated region (E6). The exons are separated by relatively short introns, the length of the longest intron being 1.3 kilobase pairs. An important difference between these and all other known class II B genes is that the beta 2 domain-encoding exon is split by an intron 97 base pairs in length. The intron is absent in other teleost fishes such as Brachydanio rerio. A change in the 3' splice site of intron 4 in some of the genes of A. hansbaenschi and of another cichlid fish, Cyphotilapia frontosa, has produced two extra codons at the 5' end of exon 5. Comparison of the A. hansbaenschi coding sequences with those of C. frontosa has revealed a concentration of variability in exon 2 and part of exon 3. Taken together, these observations provide evidence for the existence in cichlid fishes of at least two class II B loci which are functionally equivalent to the corresponding loci in mammals. The exon-intron organization and sequence similarities indicate that the two loci arose by duplication from a common ancestral gene.