Unprecedented intraspecific diversity of the MHC class I region of a teleost medaka, Oryzias latipes

The major histocompatibility complex (MHC) is present at a single chromosomal locus of all jawed vertebrate analyzed so far, from sharks to mammals, except for teleosts whose orthologs of the mammalian MHC-encoded genes are dispersed at several chromosomal loci. Even in teleosts, several class IA genes and those genes directly involved in class I antigen presentation preserve their linkage, defining the teleost MHC class I region. We determined the complete nucleotide sequence of the MHC class I region of the inbred HNI strain of medaka, Oryzias latipes (northern Japan population-derived), from four overlapping bacterial artificial chromosome (BAC) clones spanning 540,982 bp, and compared it with the published sequence of the corresponding region of the inbred Hd-rR strain of medaka (425,935 bp, southern Japan population-derived) as the first extensive study of intraspecies polymorphisms of the ectotherm MHC regions. A segment of about 100 kb in the middle of the compared sequences encompassing two class Ia genes and two immunoproteasome subunit genes, PSMB8 and PSMB10, was so divergent between these two inbred strains that a reliable sequence alignment could not be made. The rest of the compared region (about 320 kb) showed a fair correspondence, and an approximately 96% nucleotide identity was observed upon gap-free segmental alignment. These results indicate that the medaka MHC class I region contains an ∼100-kb polymorphic core, which is most probably evolving adaptively by accumulation of point mutations and extensive genetic rearrangements such as insertions, deletions and duplications. The nucleotide sequence data of HNI MHC class I region reported in this paper have been submitted to the DDBJ/EMBL/GenBank and were assigned the accession number AB183488.     

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.     

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

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

Evolution of major histocompatibility complex class I genes in Cetartiodactyls

Previous studies of cattle MHC have suggested the presence of at least four classical class I loci. Analysis of haplotypes showed that any combination of one, two or three genes may be expressed, although no gene is expressed consistently. The aim of this study was to examine the evolutionary relationships among these genes and to study their phylogenetic history in Cetartiodactyl species, including cattle and their close relatives. A secondary aim was to determine whether recombination had occurred between any of the genes. MHC class I data sets were generated from published sequences or by polymerase chain reaction from cDNA. Phylogenetic analysis revealed that MHC class I sequences from Cetartiodactyl species closely related to cattle were distributed among the main cattle gene "groups", while those from more distantly related species were either scattered (sheep, deer) or clustered in a species-specific manner (sitatunga, giraffe). A comparison between gene and species trees showed a poor match, indicating that divergence of the MHC sequences had occurred independently from that of the hosts from which they were obtained. We also found two clear instances of interlocus recombination among the cattle MHC sequences. Finally, positive natural selection was documented at positions throughout the alpha 1 and 2 domains, primarily on those amino acids directly involved in peptide binding, although two positions in the alpha 3 domain, a region generally conserved in other species, were also shown to be undergoing adaptive evolution.     

Evolutionary origin and diversification of the mammalian CD1 antigen genes.

CD1 antigens are cell-surface glycoproteins which have a molecular structure which is similar (consisting of extracellular domains alpha 1, alpha 2, and alpha 3, a transmembrane portion, and a cytoplasmic tail) to that of class I MHC molecules. Phylogenetic analysis of mammalian CD1 DNA sequences revealed that these genes are more closely related to the class I major histocompatibility complex (MHC) than to the class II MHC and that mammalian genes are more closely related to avian class I MHC genes than they are to mammalian class I MHC genes. The CD1 genes form a multigene family with different numbers of genes in different species (five in human, eight in rabbit, and two in mouse). Known CD1 genes are grouped into the following three families, on the basis of evolutionary relationship: (1) the human HCD1B gene and a partial sequence from the domestic rabbit, (2) the human HCD1A and HCD1C genes, and (3) the human HCD1D and HCD1E genes plus the two mouse genes and a sequence from the cottontail rabbit. The alpha 1 and alpha 2 domains of CD1 are much less conserved at the amino acid level than are the corresponding domains of class I MHC molecules, but the alpha 3 domain of CD1 seems to be still more conserved than the well-conserved alpha 3 domain of class I MHC molecules. Furthermore, in the human CD1 gene family, interlocus exon exchange has homogenized alpha 3 domains of all CD1 genes except HCD1C.     

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.     

Rhesus macaque class I duplicon structures, organization, and evolution within the alpha block of the major histocompatibility complex

The alpha block of the human and chimpanzee major histocompatibility complex (MHC) class I genomic region contains 10 to 11 duplicated MHC class I genes, including the HLA/Patr-A, -G, and -F genes. In comparison, the alpha block of the rhesus macaque (Macaca mulatta, Mamu) has an additional 20 MHC class I genes within this orthologous region. The present study describes the identification and analysis of the duplicated segmental genomic structures (duplicons) and genomic markers within the alpha block of the rhesus macaque and their use to reconstruct the duplication history of the genes within this region. A variety of MHC class I genes, pseudogenes, transposons, and retrotransposons, such as Alu and ERV16, were used to categorize the 28 duplicons into four distinct structural categories. The phylogenetic relationship of MHC class I genes, Alu, and LTR16B sequences within the duplicons was examined by use of the Neighbor-Joining (NJ) method. Two single-duplicon tandem duplications, two polyduplicon tandem duplications with an accompanying inversion product per duplication, eight polyduplicon tandem duplications steps, 12 deletions, and at least two recombinations were reconstructed to explain the highly complex organization and evolution of the 28 duplicons (nine inversions) within the Mamu alpha block. On the basis of the phylogenetic evidence and the reconstructed tandem duplication history of the 28 duplicons, the Mamu/Patr/HLA-F ortholog was the first MHC class I gene to have been fixed without further duplication within the alpha block of primates. Assuming that the rhesus macaque and the chimpanzee/human lineages had started with the same number of MHC class I duplicons at the time of their divergence approximately 24 to 31 MYA, then the number of genes within the alpha block have been duplicated at an approximately three times greater rate in the rhesus macaque than in either the human or chimpanzee.     

Structure of a cDNA for the pro alpha 2 chain of human type I procollagen. Comparison with chick cDNA for pro alpha 2(I) identifies structurally conserved features of the protein and the gene

Nucleotide sequences were determined for cloned cDNAs encoding for more than half of the pro alpha 2 chain of type I procollagen from man. Comparisons with previously published data on homologous cDNAs from chick embryos made it possible to examine evolution of the gene in two species which have diverged for 250-300 million years. The amino acid sequence of the alpha-chain domain supported previous indications that there is a strong selective pressure to maintain glycine as every third amino acid and to maintain a prescribed distribution of charged amino acids. However, there is little apparent selective pressure on other amino acids. The amino acid sequence of the C-propeptide domain showed less divergence than the alpha-chain domain. The 5' end or N terminus of the human C-propeptide, however, contained an insert of 12 bases coding for 4 amino acids not found in the chick C-propeptide. About 100 amino acid residues from the N terminus, two residues found in the chick sequence were missing from the human. In the second half of the C-propeptide, there was complete conservation of a 37 amino acid sequence and conservation of 50 out of 51 amino acids in the same region, an observation which suggested that the region serves some special purpose such as directing the association of one pro alpha 2(I) C-propeptide with two pro alpha 1(I) C-propeptides so as to produce the heteropolymeric structure of type I procollagen. In addition, comparison of human and chick DNAs for pro alpha 2(I) revealed three different classes of conservation of nucleotide sequence which have no apparent effect on the structure of the protein: a preference for U on the third base position of codons for glycine, proline, and alanine; a high degree of nucleotide conservation in the 51 amino acid highly conserved region of the C-propeptide; a high degree of nucleotide conservation in the 3'-noncoding region. These three classes of nucleotide conservation may reflect unusual features of collagen genes, such as their high GC content or their highly repetitive coding sequences.     

Identification of novel major histocompatibility complex class I sequences in a marsupial,the brushtail possum (<Emphasis Type="Italic">Trichosurus vulpecula</Emphasis>)

The major histocompatibility complex (MHC) is an essential part of the vertebrate immune response. MHC genes may be classified as classical, non-classical or non-functional pseudogenes. We have investigated the diversity of class I MHC genes in the brushtail possum, a marsupial native to Australia and an introduced pest in New Zealand. The MHC of marsupials is poorly characterised compared to eutherian mammal species. Comparisons between marsupials and eutherians may enhance understanding of the evolution and functions of this important genetic region. We found a high level of diversity in possum class I MHC genes. Twenty novel sequences were identified using polymerase chain reaction (PCR) primers designed from existing marsupial class I MHC genes. Eleven of these sequences shared a high level of homology with the only previously identified possum MHC class I gene TrvuUB and appear to be alleles at a single locus. Another seven sequences are also similar to TrvuUB but have frame-shift mutations or stop codons early in their sequence, suggesting they are non-functional alleles of a pseudogene locus. The remaining sequences are highly divergent from other possum sequences and clusters with American marsupials in phylogenetic analysis, indicating they may have changed little since the separation of Australian and American marsupials.     

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.     

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.Subject terms: Molecular evolution, Immunogenetics     

Fugu orthologues of human major histocompatibility complex genes: a genome survey

The major histocompatibility complex (MHC) region in fish has been subjected to piecemeal analysis centering on the in-depth characterization of single genes. The emphasis has been on those genes proven to be involved in the immune response such as the class I and class II antigen presenting genes and the complement genes. The Fugu genome data presents the opportunity to examine the short-range linkage of potentially all the human MHC orthologues and examine conserved synteny with the human and, to a more limited extent, zebrafish genomes. Analysis confirms the existence of a limited MHC locus in Fugu comprising the MHC class Ia genes and associated class II region genes involved in class I antigen presentation. Identification of additional human MHC orthologues indicates the completely dispersed nature of this region in fish, with a maximum of six MHC genes maintained within close proximity in any one contig. The majority of the other genes are present in the genome data as either singletons or pairs. Comparison with zebrafish substantiates previously observed linkages between class III region orthologues and hints at an ancient conserved class III region.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. MHC Class I genes in the Tuatara (Sphenodon spp.): evolution of the MHC in an ancient reptilian order

The major histocompatibility complex (MHC) is an extremely dynamic region of the genome, characterized by high polymorphism and frequent gene duplications and rearrangements. This has resulted in considerable differences in MHC organization and evolution among vertebrate lineages, particularly between birds and mammals. As nonavian reptiles are ancestral to both mammals and birds, they occupy an important phylogenetic position for understanding these differences. However, little is known about reptile MHC genes. To address this, we have characterized MHC class I sequences from the tuatara (Sphenodon spp.), the last survivor of an ancient order of reptiles, Sphenodontia. We isolated two different class I cDNA sequences, which share 93% sequence similarity with each other but are highly divergent from other vertebrate MHC genes. Southern blotting and polymerase chain reaction amplification of class I sequences from seven adult tuatara plus a family group indicate that these sequences represent at least two to three loci. Preliminary analysis of variation among individuals from an island population of tuatara indicates that these loci are highly polymorphic. Maximum likelihood analysis of reptile MHC class I sequences indicates that gene duplication has occurred within reptilian orders. However, the evolutionary relationships among sequences from different reptilian orders cannot be resolved, reflecting the antiquity of the major reptile lineages.     

Comparative genomic analysis of the major histocompatibility complex class I region in the teleost genus <Emphasis Type="Italic">Oryzias</Emphasis>

The major histocompatibility complex (MHC) class I region of teleosts harbors a tight cluster of the class IA genes and several other genes directly involved in class I antigen presentation. Moreover, the dichotomous haplotypic lineages (termed d- and N- lineages) of the proteasome subunit beta genes, PSMB8 and PSMB10, are present in this region of the medaka, Oryzias latipes. To understand the evolution of the Oryzias MHC class I region at the nucleotide sequence level, we analyzed bacterial artificial chromosome clones covering the MHC class I region containing the d- lineage of Oryzias luzonensis and the d- and N- lineages of Oryzias dancena. Comparison among these three elucidated sequences and the published sequences of the d- and N- lineages of O. latipes indicated that the order and orientation of the encoded genes were completely conserved among these five genomic regions, except for the class IA genes, which showed species-specific variation in copy number. The PSMB8 and PSMB10 genes showed trans-species dimorphism. The remaining regions flanking the PSMB10, PSMB8, and class IA genes showed high degrees of sequence conservation at interspecies as well as intraspecies levels. Thus, the three independent evolutionary patterns under apparently distinctive selective pressures are recognized in the Oryzias MHC class I region. Electronic Supplementary Material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Fractal landscapes and molecular evolution: modeling the myosin heavy chain gene family.

Mapping nucleotide sequences onto a "DNA walk" produces a novel representation of DNA that can then be studied quantitatively using techniques derived from fractal landscape analysis. We used this method to analyze 11 complete genomic and cDNA myosin heavy chain (MHC) sequences belonging to 8 different species. Our analysis suggests an increase in fractal complexity for MHC genes with evolution with vertebrate > invertebrate > yeast. The increase in complexity is measured by the presence of long-range power-law correlations, which are quantified by the scaling exponent alpha. We develop a simple iterative model, based on known properties of polymeric sequences, that generates long-range nucleotide correlations from an initially noncorrelated coding region. This new model-as well as the DNA walk analysis-both support the intron-late theory of gene evolution.     

Unusually limited nucleotide sequence variation of the expressed major histocompatibility complex class I genes of a New World primate species (Saguinus oedipus)

Although major histocompatibility complex (MHC) class I molecules are, as a rule, highly polymorphic in mammalian species, those of the New World primate Saguinus oedipus (cotton-top tamarin) exhibit limited polymorphism. We have cloned and sequenced twelve MHC class I cDNAs from this species. Since cloned cotton-top tamarin cell lines express three to six MHC class I molecules, this species must have at least three functional MHC class I loci. There was, however, no evidence of locus-specific substitutions in the tamarin cDNAs. Unlike all other species studied, tamarin MHC class I cDNAs displayed limited nucleotide sequence variation. The sequence similarity between the two most divergent tamarin cDNAs was 95%. To ensure that the polymerase chain reaction (PCR) primers employed in these studies had amplified all of the tamarins' expressed MHC class I genes, we used another set of primers to amplify only exons 2 and 3 from RNA and DNA. PCR of genomic DNA resulted in the amplification of six distinct clones, of which only three were well expressed. Two of these nonexpressed genes were pseudogenes and the other was a nonclassical gene. Southern blot analysis demonstrated that the tamarin has 8–11 MHC class I genes, suggesting we had indeed cloned the majority of these genes. Cotton-top tamarins are, therefore, unique among mammalian species studied to date in that they express MHC class I molecules with limited nucleotide sequence variation.The nucleotide sequence data reported in this paper have been submitted to the GenBank nucleotide sequence database and have been assigned the accession numbers M38403-15.