Characterization and polymorphic analysis of 4.5 kb genomic full-length HLA-C in the Chinese Han population

This study used long-range polymerase chain reaction to sequence 4.5 or 4.3 kb of genomic DNA covering human leukocyte antigen C (HLA-C) and its flanks in 45 Chinese Han subjects to better characterize variation in the gene in a single population. Sequences of 35 HLA-C alleles were obtained from the population, including major alleles of 13 lineages of HLA-C. Four novel alleles, C*03:04:01:02, C*04:01:01:03, C*08:22, and C*17:01:01:02, were identified, and complete full-length sequences of 18 HLA-C alleles were obtained for the first time. All sequences herein reported also represent extensions through the promoter region and the 3'-untranslated region. Fourteen 5'-nucleotide sequences and 14 3'-nucleotide sequences were detected outside the coding region. In total, 316 single-nucleotide polymorphisms unequally distributed in HLA-C subregions were observed. In addition to exons 2 and 3, nucleotide variability was found to be particularly high in exon 5, which encodes the transmembrane region. The differentiation of the C*07 and C*17 lineages in this region accounts for the high variability. The congruence of phylogeny across most regions of the gene suggests that gene conversion or recombination has not markedly influenced divergence between lineages in the evolution of HLA-C.

The surface structure of trypanosomes in relation to their molecular phylogeny

Molecular phylogenetic analysis using genes coding for ribosomal RNA and proteins suggests that trypanosomes are monophyletic. Salivarian trypanosomes showing antigenic variation of the variant surface glycoprotein (VSG) diverged from non-Salivarian trypanosomes some 200-300 million years ago. Representatives of the non-Salivarian group, the mammalian parasite, Trypanosoma cruzi, and the fresh-water fish trypanosome, T. carassii, are characterised by surfaces dominated by carbohydrate-rich mucin-like glycoproteins, which are not subject to antigenetic variation. It is suggested that this latter surface structure is typical for non-Salivarian trypanosomes as well as members of the other Kinetoplastid suborder, the Bodonina. This would imply that at some point in time in the evolution of the Salivaria the highly abundant and comparatively poorly immunogenetic mucin-like molecules must have been replaced for equally abundant but highly immunogenic VSG-like molecules. While the selective advantage for such a unique transition is difficult to imagine, the subsequent diversification of VSG genes/molecules may have been comparatively straightforward because even the most limited form of antigenic variation would have extended the duration of infection in the vertebrate and thus would have increased the chance for transfer to the vector.

On the origin of Darwin's finches

Darwin's finches comprise a group of 15 species endemic to the Galápagos (14 species) and Cocos (1 species) Islands in the Pacific Ocean. The group is monophyletic and originated from an ancestral species that reached the Galápagos Archipelago from Central or South America. Descendants of this ancestor on the Archipelago then colonized Cocos Island. In the present study, we used sequences of two mitochondrial (mt) DNA segments (922 bp of the cytochrome b gene and 1,082 bp of the control region), as well as two nuclear markers (830 bp of numt2, consisting of 140 bp of mtDNA control region and 690 bp of flanking nuclear DNA; and 740 bp of numt3, consisting of 420 bp of mt cytochrome b sequence flanked by 320 bp of nuclear DNA) to identify the species group most closely related to the Darwin's finches. To this end, we analyzed the sequences of 28 species representing the main groups (tribes) of the family Fringillidae, as well as 2 outgroup species and 13 species of Darwin's finches. In addition, we used mtDNA cytochrome b sequences of some 180 additional Fringillidae species from the database for phylogeny reconstruction by maximum-parsimony, maximum-likelihood, minimum-evolution, and neighbor-joining methods. The study identifies the grassquit genus Tiaris, and specifically the species Tiaris obscura, as the nearest living relative of Darwin's finches among the species surveyed. Darwin's finches diverged from the Tiaris group shortly after the various extant species of Tiaris diverged from one another. The initial adaptive radiation of the Tiaris group apparently occurred on the Caribbean islands and then spread to Central and South America, from where the ancestors of Darwin's finches departed for the Galápagos Islands approximately 2.3 MYA, at the time of the dramatic climatic changes associated with the closure of the Panamanian isthmus and the onset of Pleistocene glaciation.

Independent origin of functional MHC class II genes in humans and New World monkeys

In previous studies, major histocompatibility complex (MHC) class II DP, DQ, and DR families of genes were characterized in different primate species mostly on the basis of their second exon sequences. Resemblances were found between Old World monkey (OWM) and New World monkey (NWM) genes and were interpreted as being the result of transspecies evolution. Subsequent analysis of intron sequences of catarrhine and platyrrhine DRB genes, however, revealed that the amplifiable genes were not, in fact, orthologous. To test other DRB genes and other families of the class II region Southern blot hybridizations were carried out with tamarin genomic DNA using probes specific for the third exons of the tamarin DQA, DQB, DPB, and DRB genes. The hybridizing bands were extracted from the gel and the third exons of the genes were amplified by PCR, cloned, and sequenced. With two exceptions, all NWM class II genes were found to group separately from the human sequences. Only the sequences of one nonfunctional DQB locus appeared to be more closely related to human genes than to other platyrrhine DQB genes. In the DRB family one gene was found that grouped with sheep and strepsirhine DRB sequences and might represent an old gene lineage. To extend the sequences to the second exon, long PCRs were performed on tamarin genomic DNA. This approach was successful for five of the ten third exon sequences. From these data, we conclude that at least the functional MHC class II genes have expanded independently in catarrhines and platyrrhines.

The implications of intergenic polymorphism for major histocompatibility complex evolution

A systematic survey of six intergenic regions flanking the human HLA-B locus in eight haplotypes reveals the regions to be up to 20 times more polymorphic than the reported average degree of human neutral polymorphism. Furthermore, the extent of polymorphism is directly related to the proximity to the HLA-B locus. Apparently linkage to HLA-B locus alleles, which are under balancing selection, maintains the neutral polymorphism of adjacent regions. For these linked polymorphisms to persist, recombination in the 200-kb interval from HLA-B to TNF must occur at a low frequency. The high degree of polymorphism found distal to HLA-B suggests that recombination is uncommon on both sides of the HLA-B locus. The least-squares estimate is 0.15% per megabase with an estimated range from 0.02 to 0.54%. These findings place strong restrictions on possible recombinational mechanisms for the generation of diversity at the HLA-B.

Biglycan-like extracellular matrix genes of agnathans and teleosts

Biglycan and decorin are two members of a family of small extracellular matrix proteoglycans characterized by the presence of 10 leucine-rich repeats and one or two attachment sites for glucosaminoglycans. Both have thus far been described only from tetrapod species, mainly mammals. Because the extracellular matrix has played an important part in the evolution of Metazoa, the phylogeny of its components is of considerable interest. In this study, biglycan-like (BGL) cDNA sequences have been obtained from two teleost (Oreochromis cichlid and zebrafish) and two lamprey species. The analysis of the sequences suggests that, like tetrapods, the lampreys possess two types of proteoglycans, both of which are biglycan-like; decorin-like proteoglycans could not be identified in these species. The genes specifying these two types apparently arose by duplication in the lamprey lineage after its divergence from gnathostomes. The two teleost species possess a BGL proteoglycan and a bona fide decorin. The BGL proteoglycan is highly divergent from the tetrapod biglycan and related to the BGL proteoglycans of the lamprey. Hence, although the duplication generating the ancestors of biglycan and decorin genes occurred after the divergence of agnathans but before the emergence of teleosts, only decorin acquired its characteristic properties in the bony fishes. The BGL gene presumably turned into a typical biglycan only in the tetrapod lineages. The presumed acquisitions of new functions appear to have been accompanied by changes in the evolutionary rate.

Do lampreys have lymphocytes? The Spi evidence

It is generally accepted that living jawless vertebrates (lampreys and hagfishes) lack the capability of mounting an adaptive immune response. At the same time, however, there are reports describing histological evidence for the presence in agnathan tissues of lymphocytes, the key players in adaptive immunity. The question therefore arises whether the cells identified morphologically as lymphocytes are true lymphocytes in terms of their genetic developmental program. In this study, evidence is provided that the lampreys express a member of the purine box 1 (PU.1)/spleen focus-forming virus integration B (Spi-B) gene family known to be critically and specifically involved in the differentiation of lymphocytes in jawed vertebrates. The lamprey gene is expressed in the lymphocyte-like cells of the digestive tract and inexplicably also in the ovary.

Mhc class II B gene evolution in East African cichlid fishes

A distinctive feature of essential major histocompatibility complex (Mhc) loci is their polymorphism characterized by large genetic distances between alleles and long persistence times of allelic lineages. Since the lineages often span several successive speciations, we investigated the behavior of the Mhc alleles during or close to the speciation phase. We sequenced exon 2 of the class II B locus 4 from 232 East African cichlid fishes representing 32 related species. The divergence times of the (sub)species ranged from 6,000 to 8.4 million years. Two types of evolutionary analysis were used to elucidate the pattern of exon 2 sequence divergence. First, phylogenetic methods were applied to reconstruct the most likely evolutionary pathways leading from the last common ancestor of the set to the extant sequences, and to assess the probable mechanisms involved in allelic diversification. Second, pairwise comparisons of sequences were carried out to detect differences seemingly incompatible with origin by nonparallel point mutations. The analysis revealed point mutations to be the most important mechanism behind allelic divergences, with recombination playing only an auxiliary part. Comparison of sequences from related species revealed evidence of random allelic (lineage) losses apparently associated with speciation. Sharing of identical alleles could be demonstrated between species that diverged 2 million years ago. The phylogeny of the exon was incongruent with that of the flanking introns, indicating either a high degree of convergent evolution at the peptide-binding region-encoding sites, or intron homogenization.

Alu elements support independent origin of prosimian, platyrrhine, and catarrhine Mhc-DRB genes

The primate major histocompatibility complex (Mhc) genes fall into two classes and each of the classes into several families. Of the class II families, the DRB family has a long and complex evolutionary history marked by gene turnover, rearrangement, and molecular convergence. Because the history is not easily decipherable from sequences alone, Alu element insertions were used as cladistic markers to support the surmised phylogenetic relationships among the DRB genes. Intron 1 segments of 24 DRB genes from five platyrrhine species and five DRB genes from three prosimian species were amplified by PCR and cloned, and the amplification products were sequenced or PCR-typed for Alu repeats. Three Alu elements were identified in the platyrrhine and four in the prosimian DRB genes. One of the platyrrhine elements (Alu50J) is also found in the Catarrhini, whereas the other two (Alu62Sc, Alu63Sc) are restricted to the New World monkeys. Similarly, the four prosimian elements are found only in this taxon. This distribution of Alu elements is consistent with the phylogeny of the DRB genes as determined from their intron 1 sequences in an earlier and the present study. It contradicts the exon 2-based phylogeny and thus corroborates the conclusion that the evolution of DRB exon 2 sequences is, to some extent, shaped by molecular convergence. Taken together, the data indicate that each of the assemblages of DRB genes in prosimians, platyrrhines, and catarrhines is derived from a separate ancestral gene.

Convergent evolution of major histocompatibility complex molecules in humans and New World monkeys

In both Old World and New World monkeys Mhc-DRB sequences have been found which resemble human DRB1*03 and DRB3 genes in their second exon. The resemblance is shared sequence motifs and clustering of the genes or the encoded proteins in phylogenetic trees. This similarity could be due to common ancestry, convergence at the molecular level, or chance. To test which of these three explanations applies, we sequenced segments of New World monkey and macaque genes which encompass the entire second exon and large parts of both flanking introns. The test strongly supports the monophyly of New World monkey DRB intron sequences. The phylogenies of introns 1 and 2 from DRB1*03-like and DRB3-like genes are congruent, but both are incongruent with the exon 2-based phylogeny. The matching of intron 1- and intron 2-based phylogenies with each other suggests that reciprocal recombination has not played a major role in exon 2 evolution. Statistical comparisons of exon 2 from different DRB1*03 and DRB3 lineages indicate that it was neither gene conversion (descent), nor chance, but molecular convergence that has shaped their characteristic motifs. The demonstration of convergence in anthropoid Mhc-DRB genes has implications for the classification, age, and mechanism of generation of DRB allelic lineages.

Expression of the dentin matrix protein 1 gene in birds

The emergence of jawed vertebrates was predicated on the appearance of several innovations, including tooth formation. The development of teeth requires the participation of several specialized genes, in particular, those necessary for the formation of hard tissues--dentin, enamel, and cementum. Some vertebrates, most conspicuously birds, secondarily lost the tooth-forming ability. To determine the fate of some of the tooth-forming genes in the birds, we tested a domestic fowl cDNA library for the expression of the dentin matrix protein 1 (DMP1) gene. The library was prepared from the poly(A+) RNA isolated from the jaws of 11- to 13-day-old embryos and the testing was carried out by the polymerase chain reaction with degenerate primers designed on the basis of the available mammalian and reptile sequences. A chicken homologue of the DMP1 gene identified by this approach was shown to be expressed in the jaws and long bones, the same two tissues as in mammals. The chicken DMP1 gene has an exon/ intron organization similar to that of its mammalian and reptile counterparts. The chicken gene contains three short highly conserved segments, the rest of the gene being poorly alignable or not alignable with its mammalian or reptilian homologues. The distribution of similarities and dissimilarities along the gene is indicative of a mode of evolution in which only short segments are kept constant, while the rest of the gene is relatively free to vary as long as the proportion of certain amino acid residues is retained in the encoded polypeptide. The DMP1 gene may have been retained in birds because of its involvement in bone formation.

Characterization of dentin matrix protein 1 gene in crocodilia

Several clones containing DMP1 cDNA were isolated from a caiman tooth library by screening with a platypus DMP1 probe. The caiman DMP1 shows little amino acid sequence similarity to mammalian DMP1s for much of its length. A few highly conserved regions can, however, be identified that correspond to the slowly evolving parts of the corresponding mammalian genes. Southern blot analysis using probes comprising either conserved regions or longer segments of the gene indicates that only a single DMP1 locus exists. In coding regions, exon-intron boundaries and reading frames are shared by caiman and mammalian genes with the exception of exons 1 and 5, which are longer in the caiman. The repetitive sequence of the last exon is shared by mammals and caiman as are the high Ser content and acidity due to a high proportion of Asp and Glu residues. The conserved mammalian cell-attachment signal Arg-Gly-Asp is absent in the caiman DMP1. In contrast to the amelogenin gene, the DMP1 gene appears to evolve rapidly in vertebrates.

Origin of gene overlap: the case of TCP1 and ACAT2

The human acetyl-CoA acetyltransferase 2 gene, ACAT2, codes for a thiolase, an enzyme involved in lipid metabolism. The human T-complex protein 1 gene, TCP1, encodes a molecular chaperone of the chaperonin family. The two genes overlap by their 3'-untranslated regions, their coding sequences being located on opposite DNA strands in a tail-to-tail orientation. To find out how the overlap might have arisen in evolution, the homologous genes of the zebrafish, the African toad, caiman, platypus, opossum, and wallaby were identified. In each species, standard or long polymerase chain reactions were used to determine whether the ACAT2 and TCP1 homologs are closely linked and, if so, whether they overlap. The results reveal that the overlap apparently arose during the transition from therapsid reptiles to mammals and has been retained for >200 million years. Part of the overlapping untranslated region shows remarkable sequence conservation. The overlap presumably arose during the chromosomal rearrangement that brought the two unrelated and previously separated genes together. One or both of the transposed genes found by chance signals that are necessary for the processing of their transcripts to be present on the noncoding strand of the partner gene.

Phylogeny of Darwin's finches as revealed by mtDNA sequences

Darwin's finches comprise a group of passerine birds first collected by Charles Darwin during his visit to the Galápagos Archipelago. The group, a textbook example of adaptive radiation (the diversification of a founding population into an array of species differentially adapted to diverse environmental niches), encompasses 14 currently recognized species, of which 13 live on the Galápagos Islands and one on the Cocos Island in the Pacific Ocean. Although Darwin's finches have been studied extensively by morphologists, ecologists, and ethologists, their phylogenetic relationships remain uncertain. Here, sequences of two mtDNA segments, the cytochrome b and the control region, have been used to infer the evolutionary history of the group. The data reveal the Darwin's finches to be a monophyletic group with the warbler finch being the species closest to the founding stock, followed by the vegetarian finch, and then by two sister groups, the ground and the tree finches. The Cocos finch is related to the tree finches of the Galápagos Islands. The traditional classification of ground finches into six species and tree finches into five species is not reflected in the molecular data. In these two groups, ancestral polymorphisms have not, as yet, been sorted out among the cross-hybridizing species.

The dentin matrix protein 1 gene of prototherian and metatherian mammals

Mineralization of tooth dentin (the deposition of hydroxyapatite crystals in and around collagen type I fibers of the extracellular matrix) requires the involvement of several genes, among them the gene coding for the dentin matrix protein 1, DMP1. We determined the exon-intron organization of the cattle DMP1 gene and used this information to amplify by the polymerase chain reaction homologous gene fragments from the genomic DNA of two species of metatherian (marsupial) mammals and one prototherian (monotreme) species. The translated proto- and metatherian protein sequences are highly divergent from the eutherian sequences but retain the general characteristics of the DMP1 (high acidity, serine-richness, multiple glycosylation sites, and the presence of the RGD cell attachment tripeptide). They therefore appear to be functional even though, evolutionarily, teeth are in a regression phase in prototherians. It is possible, therefore, that DMP1 is also involved in other functions besides dentinogenesis. The DMP1 gene appears to evolve rapidly and apparently tolerates non-frame-shifting insertions/deletions throughout the coding sequence.

Isolation of mhc class II DMA and DMB cDNA sequences in a marsupial: the gray short-tailed opossum (Monodelphis domestica)

We report the cDNA sequences for the DMA and DMB family of Mhc genes of the gray short-tailed opossum. Until now DM sequences were available only in eutherian mammals. The marsupial sequences indicate that both members of the family are old and probably diverged from other classical class II families about the time of the radiation of jawed vertebrates some 450 million years ago. We examine the evolutionary rates of equivalent sets of classical and nonclassical genes to check for rate heterogeneity. We find the alpha-1 domain of the DR genes to be untypically conservative in its evolutionary mode. The DM genes appear to evolve at rates typical of other class II genes, indicating that their placement at the root of class II gene evolutionary trees may be justified.

Identification and characterization of amelogenin genes in monotremes, reptiles, and amphibians

Two features make the tooth an excellent model in the study of evolutionary innovations: the relative simplicity of its structure and the fact that the major tooth-forming genes have been identified in eutherian mammals. To understand the nature of the innovation at the molecular level, it is necessary to identify the homologs of tooth-forming genes in other vertebrates. As a first step toward this goal, homologs of the eutherian amelogenin gene have been cloned and characterized in selected species of monotremes (platypus and echidna), reptiles (caiman), and amphibians (African clawed toad). Comparisons of the homologs reveal that the amelogenin gene evolves quickly in the repeat region, in which numerous insertions and deletions have obliterated any similarity among the genes, and slowly in other regions. The gene organization, the distribution of hydrophobic and hydrophilic segments in the encoded protein, and several other features have been conserved throughout the evolution of the tetrapod amelogenin gene. Clones corresponding to one locus only were found in caiman, whereas the clawed toad possesses at least two amelogenin-encoding loci.

Cloning of major histocompatibility complex (Mhc) genes from threespine stickleback, Gasterosteus aculeatus

The threespine stickleback Gasterosteus aculeatus is an important model in evolutionary and ethologic studies. Its utility would greatly be increased by the availability of molecular markers distinguishing individuals and populations. Such markers can be provided by the major histocompatibility complex (Mhc) genes, which are well known for their extensive polymorphism. In the present study, both class I and class II B Mhc genes have been identified and sequenced. Fifteen distinct class I exon 2 and exon 3 sequences were obtained and assigned to 12 loci on the basis of intron 2 length differences. Some of the loci appear to be related to class I loci identified previously in cichlid fish. The intron 2 sequences and insertions/deletions in exon 2 group the loci into three families (with one family divided further into two subfamilies) derived from different ancestral genes. The ancestors presumably diverged from one another before the divergence of Gasterosteiformes from Perciformes. The 12 distinct class II B sequences may be derived from six loci, which are, however, closely related to one another in both exonic and intronic parts and may have diverged from a single common ancestor after the divergence of Gasterosteiformes from Perciformes. The intron 2 of some of the class I genes contains two microsatellites that can be used as markers, in addition to the polymorphism of the Mhc genes in their exonic regions.

The molecular phylogeny of trypanosomes: evidence for an early divergence of the Salivaria

Chronic infections with trypanosomes dwelling extracellularly in the blood and tissues of their hosts are observed in all vertebrate classes. We present here a molecular phylogenetic reconstruction of trypanosome evolution based on nucleotide sequences of small subunit rRNA genes. The evolutionary tree suggests an ancient split into one branch containing all Salivarian trypanosomes and a branch containing all non-Salivarian lineages. The latter branch splits into a clade containing bird, reptilian and Stercorarian trypanosomes infecting mammals and a clade with a branch of fish trypanosomes and a branch of reptilian/amphibian lineages. The branching order of the non-Salivarian trypanosomes supports host-parasite cospeciation scenarios, but also suggests host switches, e.g. between bird and reptilian trypanosomes. The tree is discussed in relation to the modes of adaptation that allow trypanosomes to infect immunocompetent vertebrates. Most importantly, the early divergence of the Salivarian lineages suggests that the presence of a dense proteinaceous surface coat that is subject to antigenic variation is a unique invention of this group of parasites.

Polymorphism of the HLA class II loci in Siberian populations

The populations that colonized Siberia diverged from one another in the Paleolithic and evolved in isolation until today. These populations are therefore a rich source of information about the conditions under which the initial divergence of modern humans occurred. In the present study we used the HLA system, first, to investigate the evolution of the human major histocompatibility complex (MHC) itself, and second, to reveal the relationships among Siberian populations. We determined allelic frequencies at five HLA class II loci (DRB1, DQA1, DQB1, DPA1, and DPB1) in seven Siberian populations (Ket, Evenk, Koryak, Chukchi, Nivkh, Udege, and Siberian Eskimo) by the combination of single-stranded conformational polymorphism and DNA sequencing analysis. We then used the gene frequency data to deduce the HLA class II haplotypes and their frequencies. Despite high polymorphism at four of the five loci, no new alleles could be detected. This finding is consistent with a conserved evolution of human class II MHC genes. We found a high number of HLA class II haplotypes in Siberian populations. More haplotypes have been found in Siberia than in any other population. Some of the haplotypes are shared with non-Siberian populations, but most of them are new, and some represent "forbidden" combinations of DQA1 and DQB1 alleles. We suggest that a set of "public" haplotypes was brought to Siberia with the colonizers but that most of the new haplotypes were generated in Siberia by recombination and are part of a haplotype pool that is turning over rapidly. The allelic frequencies at the DRB1 locus divide the Siberian populations into eastern and central Siberian branches; only the former shows a clear genealogical relationship to Amerinds.

Evidence for convergent evolution of A and B blood group antigens in primates

To determine whether convergent or trans-specific evolution is responsible for the persistence of the ABO polymorphism in apes, we have sequenced segments of introns 5 and 6 of the ABO gene. Four substitutions and one insertion or deletion group human A, B, and O alleles together, separate from their chimpanzee A and gorilla B counterparts. No shared substitutions support a trans-species mode of evolution for any of the alleles examined. We conclude that the A and B antigens of the chimpanzee and gorilla, respectively, have arisen by convergent evolution. Phylogenetic analysis suggests that the human A and B alleles are ancient, having diverged at least 3 million years ago. These alleles must have therefore been trans-specifically inherited within the genus Homo.

Polymorphism of the HLA-DRB1 locus in Colombian, Ecuadorian, and Chilean Amerinds

We have characterized the DRB1 genotypes in a sample of 64 South American Indians drawn from populations in Chile, Colombia, and Ecuador. No novel DRB1 alleles were found in the total of 17 different alleles characterized, indicating that rapid allelic generation does not occur at the DRB1 loci, in contrast to HLA-B. Comparison between Chilean and Colombian/Ecuadorian samples revealed no major differences in their allelic frequencies. In the combined Amerind sample the HLA-DRB1*0407 and HLA-DRB1*1402 alleles occurred in the highest frequencies (38% and 22%, respectively). Genetic distance measurement showed the HLA-DRB1 frequencies reported here to agree with findings in other Amerind groups. The high frequencies of both HLA-DRB1*0407 and HLA-DRB1*1602 alleles, in conjunction with their absence in Siberian samples, suggest that migratory groups other than Siberians may have been involved in the peopling of the Americas.

Class I mhc genes of cichlid fishes: identification, expression, and polymorphism

Cichlid fishes of the East African Rift Valley lakes constitute an important model of adaptive radiation. Explosive speciation in the Great Lakes, in some cases as recently as 12 400 years ago, generated large species flocks that have been the focus of evolutionary studies for some time. The studies have, however, been hampered by the paucity of biochemical markers for phylogenetic reconstruction. Here, we describe a set of markers which should help to alleviate this problem. They are the class I genes of the major histocompatibility complex. We provide evidence for the existence of at least 17 class I loci in cichlid fishes, and for extensive polymorphism of three of these loci. Since the polymorphism has a trans-species character, it will be possible to use it in investigating the founding events of the individual species. The sequences of the cichlid class I fishes support the monophyly of actinopterygian fish on the one hand, and of tetrapods on the other.

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.

Molecular cloning of major histocompatibility complex class II B gene cDNA from the Bengalese finch Lonchura striata

The only avian major histocompatibility complex (Mhc) genes thus far identified are from species of the relatively small order of Galliformes, while by far the largest order of Passeriformes (songbirds), containing some 60% of extant bird species, has not been studied at all in this regard. The Galliformes emerged more than 55 million years (my) ago, the Passeriformes some 25 my ago. Because of the potential for the use of Mhc genes as markers in the study of songbird populations, an attempt was made to clone class II B genes of a passeriform species, the Bengalese finch Lonchura striata acuticauda. Using a set of primers designed on the basis of known sequences, a probe corresponding to part of exon II was obtained by the polymerase chain reaction. The probe was then used to screen a Bengalese finch cDNA library and to isolate and sequence two nearly full-length clones. The sequences reveal the presence of one presumably functional class II B locus in this bird species.

Cloning and characterization of class I Mhc genes of the zebrafish, Brachydanio rerio

The zebrafish (Brachydanio rerio) offers many advantages for immunological and immunogenetic research and has the potential for becoming one of the most important nonmammalian vertebrate research models. With this in mind, we initiated a systematic study of the zebrafish major histocompatibility complex (Mhc) genes. In this report, we describe the cloning and characteristics of the zebrafish class I A genes coding for the alpha chains of the alpha beta heterodimer and thus complete the identification of all four classes and subclasses of the Mhc in this species. We describe the full class I alpha cDNA sequence as well as the exon-intron organization of the class I A genes, including intron sequences. We identify three families of class I A genes which we designate Brre-UAA, -UBA, and -UCA. The three families originated about the time of the divergence of cyprinid and salmonid fishes. All three families are members of an ancient lineage that diverged from another, older lineage also represented in cyprinid fishes before the radiation of teleost orders. The fish class I A genes therefore evolve differently from mammalian class I A genes, in which the establishment of lineages and families mostly postdates the divergence of orders.

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.

Intensity of natural selection at the major histocompatibility complex loci

Long persistence of allelic lineages, prevalence of nonsynonymous over synonymous substitutions in the peptide-binding region (PBR), and deviation from neutrality of the expected gene identity parameter F all indicate indirectly that balancing selection is operating at functional major histocompatibility complex (MHC) loci. Direct demonstrations of the existence of balancing selection at MHC loci are, however, either lacking or not fully convincing. To define the conditions under which balancing selection could be demonstrated, we estimated its intensity from the mean number of nonsynonymous substitutions, KB, at the PBR and the mutation rate mu. We compared the five available methods for estimating KB by computer simulation and chose the most reliable ones for estimation of selection intensity. For the human MHC, the selection coefficients of the HLA-A, -B, -C, -DRB1, -DQB1, -DQA1, and -DPB1 loci are 0.015, 0.042, 0.0026, 0.019, 0.0085, 0.0028, and 0.0007, respectively. This low selection intensity places severe restrictions on the possibility of measuring selection directly in vertebrate populations.

Alu elements of the primate major histocompatibility complex

The chromosomal region constituting the major histocompatibility complex (MHC) has undergone complex evolution that is often difficult to decipher. An important aid in the elucidation of the MHC evolution is the presence of Alu elements (repeats) which serve as markers for tracing chromosomal rearrangements. As the first step toward the establishment of sets of evolutionary markers for the MHC, Alu elements present in selected MHC haplotypes of the human species, the gorilla, and the chimpanzee were identified. Restriction fragments of cosmid clones from the libraries of the three species were hybridized with Alu-specific probes, Alu elements were amplified by the polymerase chain reaction, and the amplification products were sequenced. In some cases, sequences of the regions flanking the Alu elements were also obtained. Altogether, 31 new Alu elements were identified, representing six Alu subfamilies. The average density of Alu elements in the MHC is one element per four kilobases (kb) of sequence. Alu elements have apparently been inserted steadily into the MHC over the last 65 million years (my). On average, one Alu element is inserted into the primate MHC every 4 my. Analysis of the human DR3 haplotype supports its origin by duplication from an ancestral haplotype consisting of DRB1 and DRB2 genes. The sharing of an old Alu element by the DRB1 and DRB2 genes, in turn, supports their divergence from a common ancestor more than 55 my ago.

The origin of the primate Mhc-DRB genes and allelic lineages as deduced from the study of prosimians

MHC class II genes of the DRB family were partially sequenced from 10 individuals representing six species of prosimians: Galago senegalensis, G. moholi, Otolemur garnetti, Loris tardigradus, Petterus (Lemur) fulvus, and Lemur catta. Altogether, 41 different genes were discerned, all distinct from genes identified previously. Comparative analysis of the sequences has led to the following conclusions. First, the DRB loci present in human populations diverged from one another before the divergence of prosimian and anthropoid primates. Second, major allelic lineages of the DRB1 locus, such as DRB1*03 (DRB1*13) and DRB1*04, were established more than 85 million years ago. Third, the DRB6 gene was inactivated before the separation of prosimians and anthropoids, and has remained a pseudogene for more than 85 million years. Fourth, the primate DRB region is structurally and functionally unstable. In Lemur catta, for example, all DRB genes have apparently been lost and their function taken over by DOB and/or DPB genes. DRB genes are, however, present in a related species, Petterus (Lemur) fulvus. Fifth, the prosimian DRB3 genes are all inactive; their function seems to have been taken over by new genes. Sixth, several of the prosimian DRB genes and pseudogenes have recently been duplicated. In Otolemur garnetti, for example, one chromosome carries at least three copies of the DRB3 pseudogene.

The origin of the primate Mhc-DRB genes and allelic lineages as deduced from the study of prosimians

MHC class II genes of the DRB family were partially sequenced from 10 individuals representing six species of prosimians: Galago senegalensis, G. moholi, Otolemur garnetti, Loris tardigradus, Petterus (Lemur) fulvus, and Lemur catta. Altogether, 41 different genes were discerned, all distinct from genes identified previously. Comparative analysis of the sequences has led to the following conclusions. First, the DRB loci present in human populations diverged from one another before the divergence of prosimian and anthropoid primates. Second, major allelic lineages of the DRB1 locus, such as DRB1*03 (DRB1*13) and DRB1*04, were established more than 85 million years ago. Third, the DRB6 gene was inactivated before the separation of prosimians and anthropoids, and has remained a pseudogene for more than 85 million years. Fourth, the primate DRB region is structurally and functionally unstable. In Lemur catta, for example, all DRB genes have apparently been lost and their function taken over by DOB and/or DPB genes. DRB genes are, however, present in a related species, Petterus (Lemur) fulvus. Fifth, the prosimian DRB3 genes are all inactive; their function seems to have been taken over by new genes. Sixth, several of the prosimian DRB genes and pseudogenes have recently been duplicated. In Otolemur garnetti, for example, one chromosome carries at least three copies of the DRB3 pseudogene.

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.

Resolution of the HLA-DRB6 puzzle: a case of grafting a de novo-generated exon on an existing gene

HLA-DRB6, one of the human major histocompatibility complex genes, lacks exon 1, which normally codes for the leader and the first four amino acid residues of the mature protein. Because it also lacks the HLA promoter, it was surprising to find that the gene is transcribed at a low level in a chimpanzee B-lymphoblastoma cell line, in which the DRB6 homolog is truncated as in humans. The study designed to resolve the paradox has revealed that a retrovirus related to the mouse mammary tumor viruses was inserted into intron 1 of DRB6 > 23 million years ago. The insertion was either accompanied or followed by the deletion of exon 1 and the promoter region of DRB6. In the 3' long terminal repeat of the retrovirus, however, an open reading frame for a new exon arose, which codes for a sequence of mostly hydrophobic amino acid residues; the sequence could function as a leader for the truncated DRB6 gene. This new exon has a functional donor splice site at its 3' end, which enables it to be spliced in register with DRB6 exon 2. Upstream from the new exon is a promoter enabling transcription of the DRB6 gene. Besides providing an example of a de novo generation of an exon, the study suggests a potential mechanism for generating new genes through the replacement of old exons with newly generated ones.

Major-histocompatibility-complex variation in two species of cichlid fishes from Lake Malawi

Lake Malawi in eastern Africa harbors > 500 endemic species of cichlid fishes, all of which are believed to have emerged from a single founding population in the past 2 Myr. Molecular characterization of differences among the species could provide important information about the nature of speciation in the period of adaptive radiation. Because of the close relationship, however, molecular variation among the species has been difficult to ascertain. In this communication, we provide evidence for extensive differences, in major-histocompatibility-complex (Mhc) class II genes, between two related species, Pseudotropheus zebra and Melanochromis auratus. We used specific primers to amplify and sequence intron 1 and exon 2 of the class II genes from 18 individuals. Although we found 20 different sequences among the 42 that we produced, there was not a single sequence shared by the two species. Thus the study suggests that different cichlid species of Lake Malawi have different profiles of class II alleles, presumably because the polymorphism present in the ancestral founding population segregated differentially into the various species. These results make Mhc genes an important tool for elucidating speciation.

Zebrafish Mhc class II alpha chain-encoding genes: polymorphism, expression, and function

Its small size and short generation time renders the zebrafish (Brachydanio rerio) an ideal vertebrate for immunological research involving large populations. A prerequisite for this is the identification of the molecules critical for an immune response in this species. In earlier studies, we cloned the zebrafish genes coding for the beta chains of the class I and class II major histocompatibility complex (Mhc) molecules. Here, we describe the cloning of the zebrafish alpha chain-encoding class II gene, which represents the first identification of a class II A gene in teleost fishes. The gene, which is less than 3 kilobases (kb) distant from one of the beta chain-encoding genes, is approximately 1.2 kb long and consists of four exons interrupted by very short (< 200 base pairs) introns. Its organization is similar to that of the mammalian class II A genes, but its sequence differs greatly from the sequence of the latter (36% sequence similarity). Among the most conserved parts is the promoter region, which contains X, Y, and TATA boxes with high sequence similarity to the corresponding mammalian boxes. The observed striking conservation of the promoter region suggests that the regulatory system of the class II genes was established more than 400 million years ago and has, principally, remained the same ever since. Like the DMA, but unlike all other mammalian class II A genes, the zebrafish gene codes for two cysteine residues which might potentially be involved in the formation of a disulfide bond in the alpha 1 domain. The primary transcript of the gene is 1196 nucleotides long and contains 708 nucleotides of coding sequence. The gene is expressed in tissues with a high content of lymphoid/myeloid cells (spleen, pronephros, hepatopancreas, and intestine). The analyzed genomic and cDNA sequences are probably derived from different loci (their overall sequence similarity in the coding region is 73% and their 3' untranslated regions are highly divergent from each other). The genes are apparently functional. Comparison of genes from different zebrafish populations reveals high exon 2 variability concentrated in positions coding for the putative peptide-binding region. Phylogenetic analysis suggests that the zebrafish class II A genes stem from a different ancestor than the mammalian class II A genes and the recently cloned shark class II A gene.

The synonymous substitution rate of the major histocompatibility complex loci in primates

Because the divergence of many allelic lineages at the major histocompatibility complex (MHC) loci predates species divergence, standard methods of calculating synonymous substitution rates are not applicable to this system. We used three alternative methods of rate estimation: one based on the minimum number of substitutions (Dm), another on the nucleotide difference (Dxy), and the third on the net nucleotide difference (Dn). We applied these methods to the protein-encoding sequences of primate MHC class I (A, B, and C) and class II (DRB1) genes. To determine the reliability of the different estimates, we carried out computer simulation. The distribution of the estimates based on Dxy or Dn is generally much broader than that based on Dm. More importantly, the Dm-based method nearly always has the highest probability of recovering true rates, provided that Dm is not smaller than 5. Because of its desirable statistical properties, we used the Dm-based method to estimate the rate of synonymous substitutions. The rate is 1.37 +/- 0.61 for A, 1.84 +/- 0.40 for B, 3.87 +/- 1.05 for C, and 1.18 +/- 0.36 for DRB1 loci, always per site per 10(9) years. Hence despite the extraordinary polymorphism, the mutation rate at the primate MHC loci is no higher than that of other loci.

Extensive MHC variability in cichlid fishes of Lake Malawi

Lake Malawi in East Africa harbours 500-1,000 endemic species of cichlid fishes, all presumably derived by adaptive radiation from a single founding population within the past two million years. The species of this 'flock' differ strikingly in their ecology and behaviour, moderately in their external morphology and very little in their molecular characteristics. Here we describe high sequence variability of class II major histocompatibility complex genes in a sample of species from Lake Malawi. The variability provides a set of molecular markers for studying adaptive radiation and should be useful for estimating the size of the population that founded the species flock.

Characterization of two novel Ly-6 genes. Protein sequence and potential structural similarity to alpha-bungarotoxin and other neurotoxins

Genomic clones cross-hybridizing with Ly-6A.2 cDNA were isolated and characterized for functional Ly-6-related genes. Two new Ly-6 genes, designated Ly-6F.1 and Ly-6G.1, were found to have high nucleotide homology (> or = 70%) and the characteristic four exon gene organization of Ly-6A/E and Ly-6C. By a PCR-based assay, Ly-6G.1 mRNA was readily found in bone marrow, whereas Ly-6F.1 mRNA was not detected in lymphoid tissues. Thus, Ly-6G.1 represents an additional Ly-6 gene with apparent selective expression in hematopoietic cells distinct from Ly-6A/E and Ly-6C. Using the available deduced protein sequence data for mature Ly-6 proteins, searches of the database uncovered an evolutionary relationship of Ly-6 proteins with neurotoxins isolated from snake venoms. The protein sequence conservation between the two groups was selective for, but not limited to, residues in neurotoxins that have been found to be important for their tertiary structures. From this relationship, we propose a neurotoxin-like structure for Ly-6 and Ly-6-related proteins, such as CD59.

Characterization of two novel Ly-6 genes. Protein sequence and potential structural similarity to alpha-bungarotoxin and other neurotoxins

Genomic clones cross-hybridizing with Ly-6A.2 cDNA were isolated and characterized for functional Ly-6-related genes. Two new Ly-6 genes, designated Ly-6F.1 and Ly-6G.1, were found to have high nucleotide homology (&gt; or = 70%) and the characteristic four exon gene organization of Ly-6A/E and Ly-6C. By a PCR-based assay, Ly-6G.1 mRNA was readily found in bone marrow, whereas Ly-6F.1 mRNA was not detected in lymphoid tissues. Thus, Ly-6G.1 represents an additional Ly-6 gene with apparent selective expression in hematopoietic cells distinct from Ly-6A/E and Ly-6C. Using the available deduced protein sequence data for mature Ly-6 proteins, searches of the database uncovered an evolutionary relationship of Ly-6 proteins with neurotoxins isolated from snake venoms. The protein sequence conservation between the two groups was selective for, but not limited to, residues in neurotoxins that have been found to be important for their tertiary structures. From this relationship, we propose a neurotoxin-like structure for Ly-6 and Ly-6-related proteins, such as CD59.

Dating the primigenial C4-CYP21 duplication in primates

C4 and CYP21 are two adjacent, but functionally unrelated genes residing in the middle of the mammalian major histocompatibility complex (Mhc). The C4 gene codes for the fourth component of the complement cascade, whereas the CYP21 gene specifies an enzyme (cytochrome P450c21) of the glucocorticoid and mineralocorticoid pathways. The genes occur frequently in multiple copies on a single chromosome arranged in the order C4 ... CYP21 ... C4 ... CYP21. The unit of duplication (a module) is the C4-CYP21 gene pair. We sequenced the flanking regions of the C4-CYP21 modules and the intermodular regions of the chimpanzee, gorilla, and orangutan, as well as the intermodular region of an Old World monkey, the pigtail macaque. By aligning the sequences, we could identify the duplication breakpoints in these species. The breakpoint turned out to be at exactly the same position as that found previously in humans. The sequences flanking paralogous genes in the same species were found to be more similar to one another than sequences flanking orthologous genes in different species. We interpret these results as indicating that the original (primigenial) duplication occurred before the separation of apes from Old World monkeys more than 23 million years ago. The nature of the sequence at the breakpoint suggests that the duplication occurred by nonhomologous recombination. Since then, the C4-CYP21 haplotypes have been expanding and contracting by homologous crossing over which has homogenized the sequences in each species.(ABSTRACT TRUNCATED AT 250 WORDS)

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.

Mhc-DRB genes of platyrrhine primates

The two infraorders of anthropoid primates, Platyrrhini (New World monkeys) and Catarrhini (Old World monkeys and the hominoids) are estimated to have diverged from a common ancestor 37 million years ago. The major histocompatibility complex class II DRB gene and haplotype polymorphism of the Catarrhini has been characterized in several recent studies. The present study was undertaken to obtain information on the DRB polymorphism of the Platyrrhini. Fifty-five complete exon 2 DRB sequences were obtained from six species of Platyrrhini representing both the Callitrichidae and the Cebidae families. Combined with the results of a parallel contig mapping study, our data indicate that at least three loci (DRB1*03, DRB3, and DRB5) are shared by the Catarrhini and the Platyrrhini. However, the three loci are occupied by functional genes in the former infraorder and mostly by pseudogenes in the latter. Instead of the pseudogenes, the Platyrrhini have evolved a new set of apparently functional genes-DRB11 and DRB*W12 through DRB*W19, which have thus far not been found in the Catarrhini. The DRB*W13, *W14, *W15, *W17, *W18, and *W19 genes seem to be restricted to the Cebidae family, whereas the DRB*W16 locus has so far been documented in the Callitrichidae family only. The DRB alleles of the cotton-top tamarin, and perhaps also those of the common marmoset (both members of the family Callitrichidae), are characterized by low nucleotide diversity, possibly indicating that they diverged from a common ancestral gene relatively recently.

The molecular descent of the major histocompatibility complex

In the last few years, more than 500 primate major histocompatibility complex (Mhc) genes or parts thereof have been sequenced. The extraordinary sequence information is used here to draw conclusions about the manner of Mhc evolution. The Mhc genes are found to evolve at a relatively slow rate with the regularity of a clock. It takes from 1 to 6 million years for a new mutation to be incorporated into an Mhc allele, and the mutation rate is comparable to that of most other primate genes. The nonsynonymous sites coding for the peptide-binding region (PBR) are under relatively weak positive selection pressure (selection coefficient of a few percent only); the nonsynonymous non-PBR sites are under moderate negative selection pressure. The positive pressure is probably provided by parasites and is responsible for the trans-species persistence of allelic lineages at functional Mhc loci for more than 40 million years.

Cloning of the beta 2-microglobulin gene in the zebrafish

The beta 2-microglobulin (beta 2m) is a protein found in the serum in a free form and on the cell surface in a form noncovalently associated with the alpha chain of the class I major histocompatibility complex (Mhc) molecules. In mammals, the beta 2m-encoding gene (B2m) is found on a chromosome different from the Mhc proper. We have isolated and characterized the B2m gene of the zebrafish, Brachydanio rerio, family Cyprinidae. We obtained both cDNA and genomic clones of the Brre-B2m gene. The cDNA clones contained the entire coding sequence, the entire 3' untranslated (UT) region, and at least part of the 5'UT region. The genomic clone contained the entire Brre-B2m gene. The coding sequence specifies 97 amino acid residues of the mature protein so that the zebrafish beta 2m is two residues shorter than human and one residue shorter than cattle, fowl, or turkey beta 2m (codons at positions 85 and 86 have been deleted in the Brre-B2m gene). The amino acid and nucleotide sequence similarities between zebrafish and human beta 2m (B2m) are 45% and 59%, respectively. Approximately 24% of the positions are invariant and an additional 9% show only conservative substitutions in comparisons which include all known beta 2m sequences (fish, avian, and mammalian). Most of the conserved positions are in the beta strands (some 47% of the beta-strand positions are conserved in the three vertebrate classes). The Brre-B2m gene consists of four exons separated by three introns. All of the introns are considerably shorter than the corresponding introns in the mammalian B2m genes. The coding sequences of the cDNA and the genomic clones are almost identical but the sequences of the 3'UT regions differ at 1.7% of the sites, suggesting that the genes borne by these clones might have diverged at least 0.7 million years (my) ago. In contrast to the human B2m gene, the Brre-B2m gene shows no bias in the distribution of the CpG dinucleotides: the dinucleotides are distributed evenly along the entire available sequence. The haploid genome of the zebrafish contains only one copy of the B2m gene.