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 γ-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 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.     

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.     

The Mhc class II of the Black grouse (Tetrao tetrix) consists of low numbers of B and Y genes with variable diversity and expression

We found that the Black grouse (Tetrao tetrix) possess low numbers of Mhc class II B (BLB) and Y (YLB) genes with variable diversity and expression. We have therefore shown, for the first time, that another bird species (in this case, a wild lek-breeding galliform) shares several features of the simple Mhc of the domestic chicken (Gallus gallus). The Black grouse BLB genes showed the same level of polymorphism that has been reported in chicken, and we also found indications of balancing selection in the peptide-binding regions. The YLB genes were less variable than the BLB genes, also in accordance with earlier studies in chicken, although their functional significance still remains obscure. We hypothesize that the YLB genes could have been under purifying selection, just as the mammal Mhc-E gene cluster.     

Hitchhiking and recombination in birds: evidence from Mhc-linked and unlinked loci in Red-winged Blackbirds (Agelaius phoeniceus)

Hitchhiking phenomena and genetic recombination have important consequences for a variety of fields for which birds are model species, yet we know virtually nothing about naturally occurring rates of recombination or the extent of linkage disequilibrium in birds. We took advantage of a previously sequenced cosmid clone from Red-winged Blackbirds (Agelaius phoeniceus) bearing a highly polymorphic Mhc class II gene, Agph-DABI, to measure the extent of linkage disequilibrium across approximately 40 kb of genomic DNA and to determine whether non-coding nucleotide diversity was elevated as a result of physical proximity to a target of balancing selection. Application of coalescent theory predicts that the hitchhiking effect is enhanced by the larger effective population size of blackbirds compared with humans, despite the presumably higher rates of recombination in birds. We surveyed sequence polymorphism at three Mhc-linked loci occurring 1.5-40 kb away from Agph-DAB1 and found that nucleotide diversity was indistinguishable from that found at three presumably unlinked, non-coding introns (beta-actin intron 2, beta-fibrinogen intron 7 and rhodopsin intron 2). Linkage disequilibrium as measured by Lewontin's D' was found only across a few hundred base pairs within any given locus, and was not detectable among any Mhc-linked loci. Estimated rates of the per site recombination rate p derived from three different analytical methods suggest that the amounts of recombination in blackbirds are up to two orders of magnitude higher than in humans, a discrepancy that cannot be explained entirely by the higher effective population size of blackbirds relative to humans. In addition, the ratio of the number of estimated recombination events per mutation frequently exceeds 1, as in Drosophila, again much higher than estimates in humans. Although the confidence limits of the blackbird estimates themselves span an order of magnitude, these data suggest that in blackbirds the hitchhiking effect for this region is negligible and may imply that the per site per individual recombination rate is high, resembling those of Drosophila more than those of humans.     

Polymorphism across an exon-intron boundary in an avian Mhc class II B gene.

Twenty-three sequence haplotypes spanning the boundary of the second exon and intron of a red-winged blackbird Mhc class II B gene, Agph-DAB1, are presented. The polymorphism of the exon segment is distributed in two divergent allelic lineages which appear to be maintained by balancing selection. The silent nucleotide diversity of the exon (pi = 0.101) is more than five times that of the intron (pi = 0.018) and decays rapidly across the exon-intron boundary. Additionally, genealogical reconstruction indicates that divergence from a common ancestor in the exon sample is over four times that of the intron. The intron sequences reveal a pattern of polymorphism which is characteristic of directional selection, rather than a pattern expected from linkage to a balanced polymorphism. These results suggest that the evolutionary histories of these two adjacent regions have been disassociated by recombination or gene conversion. The estimated population recombination parameter between the exon and the intron is sufficiently high (4NeC = 8.545) to explain the homogenization of intron sequences. Compatibility analyses estimate that these events primarily occur from the exon-intron boundary to about 20-30 bases into the intron. Additionally, the observation that divergent exon alleles share identical intron sequence supports the conclusion of disassociation of exon and intron evolutionary histories by recombination.     

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.     

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.     

Functional characterization of a chicken major histocompatibility complex class II B gene promoter

 A 0.7 kilobase (kb) DNA fragment from the 5′ flanking region of a chicken major histocompatibility complex (MHC) class II B gene was cloned into chloramphenicol acetyltransferase (CAT) reporter vectors and was transfected into a chicken macrophage cell line that expresses a low level of MHC class II antigens. Positive orientation-dependent promoter activity of the chicken DNA was evident in a reporter construct containing an SV40 enhancer. Deletion analysis of this 0.7 kb DNA fragment revealed a short fragment in the 3′ end that was crucial for the promoter function and negative regulatory elements (NRE) located further upstream. The conserved MHC class II X and Y boxes did not have a significant effect on promoter activity. Sequence analysis of the 0.7 kb class II B gene upstream region suggests possible involvement of interferon (IFN), E twenty-six specific (ETS)-related proteins, and other factors in regulating this promoter. A chicken T-cell line culture supernatant increased surface expression of MHC class II antigens, as well as class II promoter activity, in this macrophage cell line. This first functional characterization of a chicken MHC class II B gene promoter will aid in understanding the regulatory mechanisms that control the expression of these genes. Received: 9 July 1996 / Revised: 7 October 1996     

Structural characterization of the chicken SF-1/Ad4BP gene

SF-1/Ad4BP was identified as a master regulator controlling steroidogenic P-450 genes and belongs to the steroid hormone receptor superfamily. It is expressed in the adrenal cortex, gonads, and pituitary gonadotroph. Targeted disruption of the mouse SF-1/Ad4BP gene showed that it plays a critical role in the development of the steroidogenic tissues and pituitary gonadotroph. We have recently cloned the chicken SF-1/Ad4BP cDNA and have now cloned the chicken SF-1/Ad4BP gene and analyzed its promoter activity. This gene consists of seven exons as well as mammalian counterparts and spans about 15 kb. In mice, the gene encodes another protein, ELP, but we could not find the open reading frame of ELP in the chicken SF-1/Ad4BP gene. The promoter of this gene included five putative cis elements (E, CCAAT, GC and TATA boxes and a GA-rich element), although no TATA box has been found in mammalian counterparts. The E and CCAAT boxes moderately affected promoter activity and the GA-rich element and TATA box were essential for the expression of the chicken SF-1/Ad4BP gene.     

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.     

Conserved regulatory elements in the promoter region of the N-CAM gene.

Genomic clones containing 5'-flanking sequences, the first exon, and the entire first intron from the chicken N-CAM gene were characterized by restriction mapping and DNA sequencing. A > 600-bp segment that includes the first exon is very G + C-rich and contains a large proportion of CpG dinucleotides, suggesting that it represents a CpG island. SP-1 and AP-1 consensus elements are present, but no TATA- or CCAAT-like elements were found within 300 bp upstream of the first exon. Comparison of the chicken promoter region sequence with similar regions of the human, rat, and mouse N-CAM genes revealed that some potential regulatory elements including a "purine box" seen in mouse and rat N-CAM genes, one of two homeodomain binding regions seen in mammalian N-CAM genes, and several potential SP-1 sites are not conserved within this region. In contrast, high CpG content, a homeodomain binding sequence, an SP-1 element, an octomer element, and an AP-1 element are conserved in all four genes. The first intron of the chicken gene is 38 kb, substantially smaller than the corresponding intron from mammalian N-CAM genes. Together with previous studies, this work completes the cloning of the chicken N-CAM gene, which contains at least 26 exons distributed over 85 kb.     

Molecular evolution at Mhc genes in two populations of chinook salmon Oncorhynchus tshawytscha

The DNA sequences of four exons of the MHC (major histocompatibilty complex) were examined in chinook salmon ( Oncorhynchus tshawytscha ) from an interior (Nechako River) and a coastal (Harrison River) population in the Fraser River drainage of British Columbia. Mhc class I A1, A2 and A3 sequences and a class II B1 sequence were obtained by PCR from each of 16–20 salmon from each population. The class I A1 and a pair of linked A2–A3 exons were derived from two different classical salmonid class I genes, Sasa-A and Onmy-UA , respectively. Allelic variation for B1, A1 and A2 was characterized by the high levels of nonsynonymous substitution indicative of the effects of natural selection on Mhc domains that contain peptide binding regions. The number of alleles detected at each of the four exons ranged from three ( B1 ) to 22 ( A1 ), but levels of nucleotide sequence divergence at all four exons were low relative to classical mammalian Mhc genes. The nucleotide similarity among alleles ranged between 89 and 99% over all exons, and all four domains possessed only two major sequence motifs. Allelic distributions at B1, A1 and A3 confirmed the genetic distinctiveness of the Harrison and Nechako chinook salmon populations revealed in previous studies. The two major allelic motifs of B1 and A1 segregated strongly between the populations. In spite of evidence that allelic diversity at these chinook salmon Mhc exons has been generated by selection, the level and distribution of diversity in the two salmon populations strongly reflected the demographic history of the species, which has been characterized by repeated bottlenecks and isolation-by-distance in glacial refugia.     

Analysis of part of the chicken<Emphasis Type="Italic"> Rfp-Y</Emphasis> region reveals two novel lectin genes,the first complete genomic sequence of a class I α-chain gene,a truncated class II β-chain gene,and a large CR1 repeat

The Rfp-Y region lies on the same microchromosome as the B-F/B-L region of the B complex, yet in contrast to the latter it is poorly characterised. To date it has been shown to contain at least two class I alpha-chain ( Y-F) genes, a class II B-chain gene and a C-type lectin-like gene. We describe the sequencing and analysis of some 20 kb of the Rfp-Y region, and identify several new genes. These include two novel C-type lectin-like genes ( Y-Lec1 and Y-Lec2) that differ strongly from the previously described C-type lectin-like gene found in the Rfp-Y region. We describe a complete genomic sequence of a class I alpha-chain ( Y-F) gene and its promoter from the Rfp-Y region. The predicted cDNA from this gene has high homology to the previously reported Y-F cDNAs. The promoter contains an altered enhancer A element. This portion of the Rfp-Y region also contains a truncated class II B-chain ( Y-LB) gene, as well as a large chicken repeat 1 (CR1) element.     

Natural selection of the major histocompatibility complex (Mhc) in Hawaiian honeycreepers (Drepanidinae)

The native Hawaiian honeycreepers represent a classic example of adaptive radiation and speciation, but currently face one the highest extinction rates in the world. Although multiple factors have likely influenced the fate of Hawaiian birds, the relatively recent introduction of avian malaria is thought to be a major factor limiting honeycreeper distribution and abundance. We have initiated genetic analyses of class II beta chain Mhc genes in four species of honeycreepers using methods that eliminate the possibility of sequencing mosaic variants formed by cloning heteroduplexed polymerase chain reaction products. Phylogenetic analyses group the honeycreeper Mhc sequences into two distinct clusters. Variation within one cluster is high, with dN > dS and levels of diversity similar to other studies of Mhc (B system) genes in birds. The second cluster is nearly invariant and includes sequences from honeycreepers (Fringillidae), a sparrow (Emberizidae) and a blackbird (Emberizidae). This highly conserved cluster appears reminiscent of the independently segregating Rfp-Y system of genes defined in chickens. The notion that balancing selection operates at the Mhc in the honeycreepers is supported by transpecies polymorphism and strikingly high dN/dS ratios at codons putatively involved in peptide interaction. Mitochondrial DNA control region sequences were invariant in the i'iwi, but were highly variable in the 'amakihi. By contrast, levels of variability of class II beta chain Mhc sequence codons that are hypothesized to be directly involved in peptide interactions appear comparable between i'iwi and 'amakihi. In the i'iwi, natural selection may have maintained variation within the Mhc, even in the face of what appears to a genetic bottleneck.     

The chicken proinflammatory cytokines interleukin-1beta and interleukin-6: differences in gene structure and genetic location compared with their mammalian orthologues

The genes encoding the chicken proinflammatory cytokines interleukin (IL)-1B and IL-6 were cloned, sequenced and mapped. The exon:intron structure of the coding region of chicken IL1B corresponds almost exactly to those of mammalian IL1B. As yet, we have no evidence for a 5'-UTR non-coding exon equivalent to that found in mammalian IL1B. The exon:intron structure of chicken IL6 differs from those of mammalian IL6, having one exon fewer (the first two exons in mammalian IL6 genes appear to be fused in the chicken gene). We were unable to clone or sequence the promoter of chicken IL1B. The chicken IL6 promoter shares a number of potential regulatory sequences similar to those found in the human IL6 promoter. These putative elements include (5'-3') a glucocorticoid response element (GRE), an AP-1 binding site, an NF-IL-6 binding site (albeit in the reverse orientation), an NF-kappaB binding site, a second AP-1 binding site and a TATAAA box. A further GRE, a cAMP response element and regions with homology to c-fos serum responsive elements or retinoblastoma control elements were absent. Promoter sequence polymorphisms were not identified in eight different inbred chicken lines. A restriction single-stranded conformational polymorphism was identified which enabled chicken IL1B to be genetically mapped to one end of chromosome 2. Chicken IL6 was mapped by fluorescent in situ hybridization also to chromosome 2, at an FLpter of 0.26.     

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.     

Comparative anatomy of the primate major histocompatibility complex DR subregion: evidence for combinations of DRB genes conserved across species.

The class II region of the human major histocompatibility complex (HLA) is made up of three major subregions designated DR, DQ, and DP. With the aim of gaining an insight into the evolution and stability of DR haplotypes, a total of 63 cosmid clones were isolated from the DR subregion (Gogo-DR) of a western lowland gorilla. All but one of these cosmid clones were found to fall into two clusters. The larger cluster, A, was defined by 41 overlapping cosmid clones and contained a DRB gene segment made up of exons 4 through 6 and four DRB genes, designated Gogo-DRB6, Gogo-DRB5*01, Gogo-DRB8, and Gogo-DRB3*01. The total length of this cluster was approximately 180 kb. The second cluster, B, encompassed a contiguous DNA stretch of approximately 145 kb and was composed of 21 overlapping cosmid clones. Cluster B contained three DRB genes, designated Gogo-DRB1*08, Gogo-DRB2, and Gogo-DRB3*02. One cosmid clone (WP1-9) containing a DRB pseudogene could not be linked to either cluster A or B. Neither the organization of cluster A nor that of cluster B was identical to that of known HLA-DR haplotypes. However, two gorilla DRB genes, Gogo-DRB6 and Gogo-DRB5*01, the human counterparts of which are linked in the HLA-DR2 haplotype, were found to be located next to each other in cluster A. The arrangement of the Gogo-DRB genes in cluster B, which is presumed to be the gorilla DR8 haplotype, was similar to that of HLA-DR3/DR5/DR6 haplotypes and to that of the presumed ancestral HLA-DR8 haplotype.(ABSTRACT TRUNCATED AT 250 WORDS)     

Ancestral Polymorphism of Mhc Class II Genes in Mice: Implications for Balancing Selection and the Mammalian Molecular Clock

To investigate the evolutionary dynamics at Mhc class II DR genes of mice (genus Mus), we sequenced the peptide binding regions (PBRs) of 41 DRB (=Eβ) genes and eight DRA (=Eα) genes from 15 strains representing eight species. As expected trees of these PBR sequences imply extensive maintenance of ancestral DRB alleles across species. We use a coalescent simulation model to show that the number of interspecific coalescent events (c) observed on these trees was higher than the number expected for neutral genealogies and similar sample sizes and is more consistent with balancing selection than with neutrality. Patterns of ancestral polymorphism in mouse DRB alleles were also used to examine the tempo of synonymous substitution in the PBR of mouse class II genes. Both absolute and relative rate tests on DRA and DRB genes imply increased substitution rates at two- and fourfold degnerate sites of mice and rats relative to primates, and decreased rates for the DRB genes of primates relative to ungulate and carnivore relatives. Thus rates of synonymous substitution at Mhc DR genes in mammals appear to be subject to generation time effects in ways similar to those found at other mammalian genes.