Comparative genomic analysis of two avian (quail and chicken) MHC regions

We mapped two different quail Mhc haplotypes and sequenced one of them (haplotype A) for comparative genomic analysis with a previously sequenced haplotype of the chicken Mhc. The quail haplotype A spans 180 kb of genomic sequence, encoding a total of 41 genes compared with only 19 genes within the 92-kb chicken Mhc. Except for two gene families (B30 and tRNA), both species have the same basic set of gene family members that were previously described in the chicken "minimal essential" Mhc. The two Mhc regions have a similar overall organization but differ markedly in that the quail has an expanded number of duplicated genes with 7 class I, 10 class IIB, 4 NK, 6 lectin, and 8 B-G genes. Comparisons between the quail and chicken Mhc class I and class II gene sequences by phylogenetic analysis showed that they were more closely related within species than between species, suggesting that the quail Mhc genes were duplicated after the separation of these two species from their common ancestor. The proteins encoded by the NK and class I genes are known to interact as ligands and receptors, but unlike in the quail and the chicken, the genes encoding these proteins in mammals are found on different chromosomes. The finding of NK-like genes in the quail Mhc strongly suggests an evolutionary connection between the NK C-type lectin-like superfamily and the Mhc, providing support for future studies on the NK, lectin, class I, and class II interaction in birds.     

Variation at the major histocompatibility complex in Savannah sparrows

The class I and class II genes of the major histocompatibility complex (Mhc) encode dimeric glycoproteins responsible for eliciting the adaptive immune response of vertebrates. Recent work with birds suggests that the number, size, and arrangement of these genes can differ markedly across species, although the extent of this variation, and its causes and consequences, are poorly understood. We have used a 157-base-pair (bp) portion of the second exon of a class II B gene to probe the Mhc in a free-living population of Savannah sparrows (Passerculus sandwichensis). Segregation analysis of Mhc bands suggests that class II B genes can be found in two independently assorting clusters, as previously described for domestic chickens (Gallus gallus) and ring-necked pheasants (Phasianus colchicus) but unlike gene organization in mammals. The Mhc in Savannah sparrows appears large (with many class II B genes) and variable; we found 42 unique genotypes among 48 adults breeding on Kent Island, New Brunswick, Canada in 1995. Savannah sparrows are long-distance migrants, and these results support recent predictions that migratory birds should show higher levels of Mhc polymorphism and/or a greater number of genes than sedentary species. Savannah sparrows are also socially polygynous with high levels of extra-pair paternity, suggesting that a history of sexual selection might also influence the size and/or structure of the avian Mhc.     

Evidence that the separation of Mhc class II from class I loci in the zebrafish,Danio rerio,occurred by translocation

In the zebrafish, Danio rerio, and other teleosts, the class I and class II loci of the major histocompatibility complex ( Mhc) reside on different chromosomes. To shed light on the events that might have generated this difference from tetrapods, in which these two types of loci are clustered in a single chromosomal region, the organization of the class II loci in linkage group 8 of the zebrafish was determined by the characterization of contigs of PAC clones. Three contigs were defined: DAB, DCB, and DBB. The 350-kb-long DAB contig contained only four genes: DDB, DAB, SLC7A4, and DAA. The 150-kb-long DCB contig contained the DCB, DCA, and fz10 genes at an undetermined distance from the DAB contig. And the 120-kb-long DBB contig comprised the DBB gene presumably in another linkage group. The low gene density of the linkage group 8 contigs, contrasting with the high gene density of the zebrafish class I region, and the close association with genes [ SLC7A4 coding for an amino acid transporter, and fz10 (frizzled 10) coding for a receptor of the WNT glycoprotein] that are not linked with the tetrapod Mhc, is interpreted to mean that the separation of the class II from class I loci in teleosts occurred by translocation rather than by genomic or chromosomal duplication.     

A contig map of the Mhc class I genomic region in the zebrafish reveals ancient synteny

In contrast to the human and mouse Mhc, in which the clusters of class I and class II loci reside in close vicinity to one another, in the zebrafish, Danio rerio, they are found in different linkage groups. Chromosome walking using BAC (bacterial artificial chromosome) and PAC (P1 artificial chromosome) clones reveals the zebrafish class I region to occupy a segment of approximately 450 kb and to encompass at least 19 loci. These include three class I (Dare-UDA, -UEA, -UFA), five proteasome subunit beta (PSMB8, -9A, -9C, -11, -12), two TAPs (TAP2A, TAP2B), and one TAP binding protein (TAPBP). This arrangement contrasts with the arrangements found in human and mouse Mhc, in which the orthologues of the PSMB, TAP, and TAPBP loci reside within the class II region. In addition to this main zebrafish class I contig, a shorter contig of about 150 kb contains two additional class I (UBA, UCA) and at least five other loci. It probably represents a different haplotype of part of the class I region. The previously identified UAA gene shares an identical 5' part with UEA, but the two genes differ in their 3' parts. One of them is probably the result of an unequal crossing over. The described organization has implications for the persistence of syntenic relationships, coevolution of loci, and interpretation of the origin of the human/mouse Mhc organization.     

Pinpointing a selective sweep to the chimpanzee MHC class I region by comparative genomics

Chimpanzees experienced a reduction of the allelic repertoire at the major histocompatibility complex (MHC) class I A and B loci, which may have been caused by a retrovirus belonging to the simian immunodeficiency virus (SIV) family. Extended MHC haplotypes were defined in a pedigreed chimpanzee colony. Comparison of genetic variation at microsatellite markers mapping inside and outside the Mhc region was carried out in humans and chimpanzees to investigate the genomic extent of the repertoire reduction. Multilocus demographic analyses underscored that chimpanzees indeed experienced a selective sweep that mainly targeted the chromosomal segment carrying the Mhc class I region. Probably due to genetic linkage, the sweep also affected other polymorphic loci, mapping in the close vicinity of the Mhc class I region genes. Nevertheless, although the allelic repertoire at particular Mhc class I and II loci appears to be limited, naturally occurring recombination events allowed the establishment of haplotype diversity after the sweep. However, recombination did not have sufficient time to erase the signal of the selective sweep.     

Songbird genomics: analysis of 45 kb upstream of a polymorphic Mhc class II gene in red-winged blackbirds (Agelaius phoeniceus)

Here we present the sequence of a 45 kb cosmid containing a previously characterized poly-morphic Mhc class II B gene (Agph-DAB1) from the red-winged blackbird (Agelaius phoeniceus). We compared it with a previously sequenced cosmid from this species, revealing two regions of 7.5 kb and 13.0 kb that averaged greater than 97% similarity to each another, indicating a very recent shared duplication. We found 12 retroelements, including two chicken repeat 1 (CR1) elements, constituting 6.4% of the sequence and indicating a lower frequency of retroelements than that found in mammalian genomic DNA. Agph-DAB3, a new class II B gene discovered in the cosmid, showed a low rate of polymorphism and may be functional. In addition, we found a Mhc class II B gene fragment and three genes likely to be functional (encoding activin receptor type II, a zinc finger, and a putative gamma-filamin). Phylogenetic analysis of exon 2 alleles of all three known blackbird Mhc genes indicated strong clustering of alleles by locus, implying that large amounts of interlocus gene conversion have not occurred since these genes have been diverging. Despite this, interspecific comparisons indicate that all three blackbird Mhc genes diverged from one another less than 35 million years ago and are subject to concerted evolution in the long term. Comparison of blackbird and chicken Mhc promoter regions revealed songbird promoter elements for the first time. The high gene density of this cosmid confirms similar findings for the chicken Mhc, but the segment duplications and diversity of retroelements resembles mammalian sequences.     

A high‐resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen‐presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene‐dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle‐derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH₅₀₀₀), we generated a high‐resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR₅₀₀₀) was aligned with the cattle MHC RH map (cR₁₂₀₀₀) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two‐point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single‐ and double‐color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups.     

Evolutionary expansion of Mhc class I loci in the mole-rat, Spalax ehrenbergi

A cosmid genomic library was prepared from a single individual of the rodent Spalax ehrenbergi, the mole rat, captured in Israel. The library was screened with a mouse probe hybridizing with all mouse class I major-histocompatibility-complex (Mhc) genes; the cross-hybridizing clones were isolated; and their restriction maps were prepared using five enzymes. A total of 93 class I-bearing clones could be identified in the library. Forty-five of these clones showed partial overlaps and could be arranged into 14 clusters. Eleven of these clusters could be shown to contain two class I genes each; the remaining clusters, as well as most of the non-overlapping clones, each contained one class I gene. After the elimination of clones with possible cloning artifacts and of clones that may carry allelic forms of a given gene in the heterozygous animal, the total number of class I loci identified in Spalax is approximately 65. The high number of loci probably arose from the duplication of either the entire class I set or the different class I families. The high number of gene copies might represent a means of selecting different functional genes from the family in different mammalian orders. Three of the approximately 65 Spalax class I genes cross-hybridize with a probe specific for the mouse K, D, and L genes; two of these genes are in the same cluster. These three elements might possibly be the functional class I genes of the mole rat.     

Mhc diversity in two passerine birds: no evidence for a minimal essential Mhc

Humans express an array of Mhc genes, while the chicken has an Mhc that is relatively small and compact with fewer expressed genes. Here we ask whether the "minimal essential Mhc" of the chicken is representative for birds. We investigated the RFLP genotypes in 55 great reed warblers Acrocephalus arundinaceus and 10 willow warblers Phylloscopus trochilus to obtain an overview of the number of class II B genes. There were 13-17 bands per individual in the great reed warblers and 25-30 in the willow warblers, and every individual had a unique RFLP genotype. The high number of RFLP bands indicates that both species have a large number of class II B genes although some may be pseudogenes. Seven different class II B sequences were detected in a great reed warbler cDNA library. There was considerable sequence divergence between the cDNA sequences in exon 2 (peptide-binding region, PBR), whereas they were very similar in exon 3. The cDNA sequences were easily alignable to a classical chicken class II B sequence, and balancing selection was acting in the PBR. One of the cDNA sequences had two deletions and is likely nonfunctional. Finally, the polymorphic class I and class II B RFLP fragments seemed to be linked in the five studied great reed warbler families. These and previous results suggest that birds of the order Passeriformes in general have more Mhc class I and II B genes than birds of the order Galliformes. This difference could be caused by their phylogenetic past, and/or by variance in the selection pressure for maintaining a high number of Mhc genes.     

Mhc-DRB genes of the pigtail macaque (Macaca nemestrina): implications for the evolution of human DRB genes.

The DRB family of human class II major histocompatibility complex (Mhc) loci is unusual in that individuals differ in the number and combination of genes (haplotypes) they carry. Indications are that both the allelic and haplotype polymorphisms of the DRB loci predate speciation. Searching for the evolutionary origins of these polymorphisms, we have sequenced five DRB clones isolated from a cDNA library of a pigtail macaque (Macaca nemestrina) B lymphocyte line. The clones represent five different genes which we designate Mane-DRB*01-Mane-DRB*05. The genes appears to be approximately equidistant from each other, so that allelic relationships between them cannot be established on the basis of the sequence data alone. If positions coding for the peptide-binding region of the class II beta chains are eliminated from sequence comparisons, the Mane-DRB genes appear to be most closely related to the human (HLA) DRB1 genes of the DRw52 group. We interpret this finding to indicate that the ancestral gene of the DRw52 group of human DRB1 alleles separated from the rest of the HLA-DRB1 alleles before the separation of the Old World monkeys (Cercopithecoidea) from the apes (Hominoidea) in the early Oligocene. After this separation, the ancestral DRB1 gene of the DRw52 group duplicated in the Old World monkey lineage to give rise to genes at three loci at least, while in the ape lineage this gene may have remained single and diverged into a number of alleles instead. These findings suggest that some of the polymorphism currently present at the DRB1 locus is greater than 35 Myr old.     

Mhc class I genes of swordtail fishes, Xiphophorus: variation in the number of loci and existence of ancient gene families

Swordtail fishes and platies in the genus Xiphophorus (order Cyprinodontiformes, Teleostei) encompass 22 closely related species which are the products of a recent adaptive radiation in the streams of Central America. To investigate the evolution of the major histocompatibility complex (Mhc) genes in the period immediately following speciation, the class I genes from 20 of the 22 species were cloned and characterized by sequencing. The analysis revealed the existence of multiple loci (at least seven in some individuals) whose numbers vary among the different species and probably also among individuals of the same species. The variation does not seem to bear any relationship to the taxonomy of the genus. Genes at the different loci are distinguished by their intron sequences and by the presence of characteristic motifs in exons 2 and 3. The variation in copy number of loci may have been effected in part by unequal crossing over occurring between introns of misaligned closely related genes. The sequences of the genes fall into two groups, A and B, which represent ancient lineages. The groups define two families of loci, which diverged from each other an estimated 85 million years ago, before the separation of the Acanthopterygii from the Paracanthopterygii of the advanced bony fishes. Evolution of the genes within each family can be explained by the birth-and-death process driven by gene duplications and mutational differentiation.     

Conserved haplotype blocks within the sheep MHC and low SNP heterozygosity in the Class IIa subregion

This report describes single-nucleotide polymorphisms (SNPs) in the sheep major histocompatibility complex (MHC) class II and class III regions and provides insights into the internal structure of this important genomic complex. MHC haplotypes were deduced from sheep family trios based on genotypes from 20 novel SNPs representative of the class II region and 10 previously described SNPs spanning the class III region. All 30 SNPs exhibited Hardy-Weinberg proportions in the sheep population studied. Recombination within an extended sire haplotype was observed within the class II region for 4 of 20 sheep chromosomes, thereby supporting the presence of separated IIa and IIb subregions similar to those present in cattle. SNP heterozygosity varied across the class II and III regions. One segment of the class IIa subregion manifested very low heterozygosity for several SNPs spanning approximately 120 Kbp. This feature corresponds to a subregion within the human MHC class II region previously described as a 'SNP desert' because of its paucity of SNPs. Linkage disequilibrium (LD) was reduced at the junction separating the putative class IIb and IIa subregions and also between the class IIa and the class III subregions. The latter observation is consistent with either an unmapped physical separation at this location or more likely a boundary characterized by more frequent recombination between two conserved subregions, each manifesting high within-block LD. These results identify internal blocks of loci in the sheep MHC, within which recombination is relatively rare.     

Major histocompatibility complex of the mole-rat

The major histocompatibility complex (Mhc) is a group of loci coding for lymphocyte membrane glycoproteins that provide the context for the recognition of foreign antigens in the initial phase of the immune response. The complex contains a large number of loci, some of which are highly polymorphic. The complexity and polymorphism pose a number of questions concerning the evolution of the Mhc. In an attempt to answer some of these questions, we have begun to study the Mhc of the mole-rat, Spalax ehrenbergi, a rodent representing a complex of sibling species occupying ecologically and geographically clearly delineated regions within the borders of Israel. In an earlier publication we identified the Spalax major histocompatibility (Smh) complex serologically and biochemically. Here, we analyze the Smh by Southern blotting of DNA fragments produced by restriction enzyme digestion. The fragments were hybridized to mouse probes specific for class I, class II, and C4 genes. The analysis has revealed that the Smh complex contains as many class I genes as the mouse does and that these genes are polymorphic. The number of class II genes could not be determined with certainty, but it is probably not greater than in the mouse. Polymorphism was also detected at the loci coding for the complement component 4 (C4), which are probably closely linked to the Smh complex. The polymorphism of mole-rat class I loci contrasts with the reported monomorphism of these loci in the Syrian hamster. Since the mole-rat leads a solitary, subterranean life, as the Syrian hamster does, ecology cannot be an explanation for the lack of class I polymorphism in the latter species.On leave from the Department of Physiology, University of Zagreb Medical Faculty, Zagreb, Yugoslavia.     

A 39-kb sequence around a blackbird Mhc class II gene: ghost of selection past and songbird genome architecture

To gain an understanding of the evolution and genomic context of avian major histocompatibility complex (Mhc) genes, we sequenced a 38.8-kb Mhc-bearing cosmid insert from a red-winged blackbird (Agelaius phoeniceus). The DNA sequence, the longest yet retrieved from a bird other than a chicken, provides a detailed view of the process of gene duplication, divergence, and degeneration ("birth and death") in the avian Mhc, as well as a glimpse into major noncoding features of a songbird genome. The peptide-binding region (PBR) of the single Mhc class II B gene in this region, Agph-DAB2, is almost devoid of polymorphism, and a still-segregating single-base-pair deletion and other features suggest that it is nonfunctional. Agph-DAB2 is estimated to have diverged about 40 MYA from a previously characterized and highly polymorphic blackbird Mhc gene, Aph-DAB1, and is therefore younger than most mammalian Mhc paralogs and arose relatively late in avian evolution. Despite its nonfunctionality, Agph-DAB2 shows very high levels of nonsynonymous divergence from Agph-DAB1 and from reconstructed ancestral sequences in antigen-binding PBR codons-a strong indication of a period of adaptive divergence preceding loss of function. We also found that the region sequenced contains very few other unambiguous genes, a partial Mhc- class II gene fragment, and a paucity of simple-sequence and other repeats. Thus, this sequence exhibits some of the genomic streamlining expected for avian as compared with mammalian genomes, but is not as densely packed with functional genes as is the chicken Mhc.     

Analysis of zebrafish Mhc using BAC clones

 Although major histocompatibility complex (Mhc) genes have been identified in a number of species, little is yet known about their organization in species other than human and mouse. The zebrafish, Danio rerio, is a good candidate for full elucidation of the organization of its Mhc. As a step toward achieving this goal, a commercially available zebrafish BAC library was screened with probes specific for previously identified zebrafish class I and class II genes, as well as for genes controlling the proteasome subunits LMP7 and LMP2. Restriction maps of the individual positive clones were prepared and the Mhc (LMP7) genes localized to specific fragments. The total length of genomic DNA fragments with Mhc genes was approximately 1700 kilobases (kb) (200 kb of fragments bearing class I loci and 1500 kb of fragments bearing class II loci). One of the two class I loci (Dare-UCA) is closely associated with the LMP7 locus; the second class I locus (Dare-UAA) is more than 50 kb distant from the UCA locus and has no LMP genes associated with it. None of the class II genes are linked to the class I or the LMP genes. All six of the previously identified class II B genes and one of the three class II A genes were found to be present in the BAC clones; no new Mhc loci could be identified in the library. Each of the six previously identified class II B loci was found to be borne by a separate group of BAC clones. The Dare-DAB and -DAA loci were found on the same clone, approximately 15 kb apart from each other. An expansion of DCB and DDB loci was detected: the zebrafish genome may contain at least five closely related DCB and two closely related DDB loci which are presumably the products of relatively recent tandem duplication. These results are consistent with linkage studies and indicate that in the zebrafish, the class I and class II loci are on different chromosomes, and the class II loci are in three different regions, at least two of which are on different chromosomes. Received: 14 August 1997 / Revised: 16 September 1997     

Comparative analysis of the bovine MHC class IIb sequence identifies inversion breakpoints and three unexpected genes

The bovine major histocompatibility complex (MHC) or BoLA is organized differently from typical mammalian MHCs in that a large portion of the class II region, called class IIb, has been transposed to a position near the centromere on bovine chromosome 23. Gene mapping indicated that the rearrangement resulted from a single inversion, but the boundaries and gene content of the inverted segment have not been fully determined. Here, we report the genomic sequence of BoLA IIb. Comparative sequence analysis with the human MHC revealed that the proximal inversion breakpoint occurred approximately 2.5 kb from the 3' end of the glutamate-cysteine ligase, catalytic subunit (GCLC) locus and that the distal breakpoint occurred about 2 kb from the 5' end from a divergent class IIDRbeta-like sequence designated DSB. Gene content, order and orientation of BoLA IIb are consistent with the single inversion hypothesis when compared with the corresponding region of the human class II MHC (HLA class II). Differences with HLA include the presence of a single histone H2B gene located between the proteasome subunit, beta type, 9 (PSMB9) and DMB loci and a duplicated TAP2 with a variant splice site. BoLA IIb spans approximately 450 kb DNA, with 20 apparently intact genes and no obvious pseudogenes. The region contains 227 simple sequence repeats (SSRs) and approximately 167 kb of retroviral-related repetitive DNA. Nineteen of the 20 genes identified in silico are supported by bovine EST data indicating that the functional gene content of BoLA IIb has not been diminished because it has been transposed from the remainder of BoLA genes.     

Persistence of Mhc heterozygosity in homozygous clonal killifish,Rivulus marmoratus: implications for the origin of hermaphroditism

The mangrove killifish Rivulus marmoratus, a neotropical fish in the order Cyprinodontiformes, is the only known obligatorily selfing, synchronous hermaphroditic vertebrate. To shed light on its population structure and the origin of hermaphroditism, major histocompatibility complex (Mhc) class I genes of the killifish from seven different localities in Florida, Belize, and the Bahamas were cloned and sequenced. Thirteen loci and their alleles were identified and classified into eight groups. The loci apparently arose approximately 20 million years ago (MYA) by gene duplications from a single common progenitor in the ancestors of R. marmoratus and its closest relatives. Distinct loci were found to be restricted to different populations and different individuals in the same population. Up to 44% of the fish were heterozygotes at Mhc loci, as compared to near homozygosity at non-Mhc loci. Large genetic distances between some of the Mhc alleles revealed the presence of ancestral allelic lineages. Computer simulation designed to explain these findings indicated that selfing is incomplete in R. marmoratus populations, that Mhc allelic lineages must have diverged before the onset of selfing, and that the hermaphroditism arose in a population containing multiple ancestral Mhc lineages. A model is proposed in which hermaphroditism arose stage-wise by mutations, each of which spread through the entire population and was fixed independently in the emerging clones.