Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution.

BACKGROUND: Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations. RESULTS: Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations. CONCLUSIONS: We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

Human and Rhesus Macaque KIR Haplotypes Defined by Their Transcriptomes.

The killer-cell Ig-like receptors (KIRs) play a central role in the immune recognition in infection, pregnancy, and transplantation through their interactions with MHC class I molecules. KIR genes display abundant copy number variation as well as high levels of polymorphism. As a result, it is challenging to characterize this structurally dynamic region. KIR haplotypes have been analyzed in different species using conventional characterization methods, such as Sanger sequencing and Roche/454 pyrosequencing. However, these methods are time-consuming and often failed to define complete haplotypes, or do not reach allele-level resolution. In addition, most analyses were performed on genomic DNA, and thus were lacking substantial information about transcription and its corresponding modifications. In this paper, we present a single-molecule real-time sequencing approach, using Pacific Biosciences Sequel platform to characterize the KIR transcriptomes in human and rhesus macaque (Macaca mulatta) families. This high-resolution approach allowed the identification of novel Mamu-KIR alleles, the extension of reported allele sequences, and the determination of human and macaque KIR haplotypes. In addition, multiple recombinant KIR genes were discovered, all located on contracted haplotypes, which were likely the result of chromosomal rearrangements. The relatively high number of contracted haplotypes discovered might be indicative of selection on small KIR repertoires and/or novel fusion gene products. This next-generation method provides an improved high-resolution characterization of the KIR cluster in humans and macaques, which eventually may aid in a better understanding and interpretation of KIR allele-associated diseases, as well as the immune response in transplantation and reproduction.

Analysis of macaque BTN3A genes and transcripts in the extended MHC: conserved orthologs of human γδ T cell modulators.

Butyrophilins (BTN), specifically BTN3A, play a central role in the modulation of γδ T cells, which are mainly present in gut and mucosal tissues. BTN3A1 is known, for example, to activate Vγ9Vδ2 T cells by means of a phosphoantigen interaction. In the extended HLA region, three genes are located, designated BTN3A1, BTN3A2 and BTN3A3, which were also defined in rhesus macaques. In contrast to humans, rhesus monkeys have an additional gene, BTN3A3Like, which has the features of a pseudogene. cDNA analysis of 32 Indian rhesus and 16 cynomolgus macaques originating from multiple-generation families revealed that all three genes are oligomorphic, and the deduced amino acids display limited variation. The macaque BTN3A alleles segregated together with MHC alleles, proving their location in the extended (Major Histocompatibility Complex) MHC. BTN3A nearly full-length transcripts of macaques and humans cluster tightly together in the phylogenetic tree, suggesting that the genes represent true orthologs of each other. Despite the limited level of polymorphism, 15 Mamu- and 14 Mafa-BTN3A haplotypes were defined, and, as in humans, all three BTN3A genes are transcribed in PBMCs and colon tissues. In addition to regular full-length transcripts, a high number of various alternative splicing (AS) products were observed for all BTN3A alleles, which may result in different isoforms. The comparable function of certain subsets of γδ T cells in human and non-human primates in concert with high levels of sequence conservation observed for the BTN3A transcripts presents the opportunity to study these not yet well understood molecules in macaques as a model species.

Determining Mhc-DRB profiles in wild populations of three congeneric true lemur species by noninvasive methods.

The major histocompatibility complex (MHC) is a highly polymorphic and polygenic genomic region that plays a crucial role in immune-related diseases. Given the need for comparative studies on the variability of immunologically important genes among wild populations and species, we investigated the allelic variation of MHC class II DRB among three congeneric true lemur species: the red-fronted lemur (Eulemur rufifrons), red-bellied lemur (Eulemur rubriventer), and black lemur (Eulemur macaco). We noninvasively collected hair and faecal samples from these species across different regions in Madagascar. We assessed DRB exon 2 polymorphism with a newly developed primer set, amplifying nearly all non-synonymous codons of the antigen-binding sites. We defined 26 DRB alleles from 45 individuals (17 alleles from E. rufifrons (N = 18); 5 from E. rubriventer (N = 7); and 4 from E. macaco (N = 20). All detected alleles are novel and show high levels of nucleotide (26.8%) and non-synonymous codon polymorphism (39.4%). In these lemur species, we found neither evidence of a duplication of DRB genes nor a sharing of alleles among sympatric groups or allopatric populations of the same species. The non-sharing of alleles may be the result of a geographical separation over a long time span and/or different pathogen selection pressures. We found dN/dS rates > 1 in the functionally important antigen recognition sites, providing evidence for balancing selection. Especially for small and isolated populations, quantifying and monitoring DRB variation are recommended to establish successful conservation plans that mitigate the possible loss of immunogenetic diversity in lemurs.

A Specialist Macaque MHC Class I Molecule with HLA-B*27-like Peptide-Binding Characteristics.

In different macaque species, the MHC A2*05 gene is present in abundance, and its gene products are characterized by low cell-surface expression and a highly conserved peptide-binding cleft. We have characterized the peptide-binding motif of Mamu-A2*05:01, and elucidated the binding capacity for virus-derived peptides. The macaque A2*05 allotype prefers the basic amino acid arginine at the second position of the peptide, and hydrophobic and polar amino acids at the C-terminal end. These preferences are shared with HLA-B*27 and Mamu-B*008, molecules shown to be involved in elite control in human HIV type 1 and macaque SIV infections, respectively. In contrast, however, Mamu-A2*05 preferentially binds 8-mer peptides. Retention in the endoplasmic reticulum seems to be the cause of the lower cell-surface expression. Subsequent peptide-binding studies have illustrated that Mamu-A2*05:01 is able to bind SIV-epitopes known to evoke a strong CD8+ T cell response in the context of the Mamu-B*008 allotype in SIV-infected rhesus macaques. Thus, the macaque A2*05 gene encodes a specialized MHC class I molecule, and is most likely transported to the cell surface only when suitable peptides become available.

Co-evolution of the MHC class I and KIR gene families in rhesus macaques: ancestry and plasticity.

Researchers dealing with the human leukocyte antigen (HLA) class I and killer immunoglobulin receptor (KIR) multi-gene families in humans are often wary of the complex and seemingly different situation that is encountered regarding these gene families in Old World monkeys. For the sake of comparison, the well-defined and thoroughly studied situation in humans has been taken as a reference. In macaques, both the major histocompatibility complex class I and KIR gene families are plastic entities that have experienced various rounds of expansion, contraction, and subsequent recombination processes. As a consequence, haplotypes in macaques display substantial diversity with regard to gene copy number variation. Additionally, for both multi-gene families, differential levels of polymorphism (allelic variation), and expression are observed as well. A comparative genetic approach has allowed us to answer questions related to ancestry, to shed light on unique adaptations of the species' immune system, and to provide insights into the genetic events and selective pressures that have shaped the range of these gene families.

MHC class I diversity of olive baboons (Papio anubis) unravelled by next-generation sequencing.

The olive baboon represents an important model system to study various aspects of human biology and health, including the origin and diversity of the major histocompatibility complex. After screening of a group of related animals for polymorphisms associated with a well-defined microsatellite marker, subsequent MHC class I typing of a selected population of 24 animals was performed on two distinct next-generation sequencing (NGS) platforms. A substantial number of 21 A and 80 B transcripts were discovered, about half of which had not been previously reported. Per animal, from one to four highly transcribed A alleles (majors) were observed, in addition to ones characterised by low transcripion levels (minors), such as members of the A*14 lineage. Furthermore, in one animal, up to 13 B alleles with differential transcription level profiles may be present. Based on segregation profiles, 16 Paan-AB haplotypes were defined. A haplotype encodes in general one or two major A and three to seven B transcripts, respectively. A further peculiarity is the presence of at least one copy of a B*02 lineage on nearly every haplotype, which indicates that B*02 represents a separate locus with probably a specialistic function. Haplotypes appear to be generated by recombination-like events, and the breakpoints map not only between the A and B regions but also within the B region itself. Therefore, the genetic makeup of the olive baboon MHC class I region appears to have been subject to a similar or even more complex expansion process than the one documented for macaque species.

Complex MHC Class I Gene Transcription Profiles and Their Functional Impact in Orangutans.

MHC haplotypes of humans and the African great ape species have one copy of the MHC-A, -B, and -C genes. In contrast, MHC haplotypes of orangutans, the Asian great ape species, exhibit variation in the number of gene copies. An in-depth analysis of the MHC class I gene repertoire in the two orangutan species, Pongo abelii and Pongo pygmaeus, is presented in this article. This analysis involved Sanger and next-generation sequencing methodologies, revealing diverse and complicated transcription profiles for orangutan MHC-A, -B, and -C. Thirty-five previously unreported MHC class I alleles are described. The data demonstrate that each orangutan MHC haplotype has one copy of the MHC-A gene, and that the MHC-B region has been subject to duplication, giving rise to at least three MHC-B genes. The MHC-B*03 and -B*08 lineages of alleles each account for a separate MHC-B gene. All MHC-B*08 allotypes have the C1-epitope motif recognized by killer cell Ig-like receptor. At least one other MHC-B gene is present, pointing to MHC-B alleles that are not B*03 or B*08. The MHC-C gene is present only on some haplotypes, and each MHC-C allotype has the C1-epitope. The transcription profiles demonstrate that MHC-A alleles are highly transcribed, whereas MHC-C alleles, when present, are transcribed at very low levels. The MHC-B alleles are transcribed to a variable extent and over a wide range. For those orangutan MHC class I allotypes that are detected by human monoclonal anti-HLA class I Abs, the level of cell-surface expression of proteins correlates with the level of transcription of the allele.

Genome typing of nonhuman primate models: implications for biomedical research.

The success of personalized medicine rests on understanding the genetic variation between individuals. Thus, as medical practice evolves and variation among individuals becomes a fundamental aspect of clinical medicine, a thorough consideration of the genetic and genomic information concerning the animals used as models in biomedical research also becomes critical. In particular, nonhuman primates (NHPs) offer great promise as models for many aspects of human health and disease. These are outbred species exhibiting substantial levels of genetic variation; however, understanding of the contribution of this variation to phenotypes is lagging behind in NHP species. Thus, there is a pivotal need to address this gap and define strategies for characterizing both genomic content and variability within primate models of human disease. Here, we discuss the current state of genomics of NHP models and offer guidelines for future work to ensure continued improvement and utility of this line of biomedical research.

Major histocompatibility complex haplotyping and long-amplicon allele discovery in cynomolgus macaques from Chinese breeding facilities.

Very little is currently known about the major histocompatibility complex (MHC) region of cynomolgus macaques (Macaca fascicularis; Mafa) from Chinese breeding centers. We performed comprehensive MHC class I haplotype analysis of 100 cynomolgus macaques from two different centers, with animals from different reported original geographic origins (Vietnamese, Cambodian, and Cambodian/Indonesian mixed-origin). Many of the samples were of known relation to each other (sire, dam, and progeny sets), making it possible to characterize lineage-level haplotypes in these animals. We identified 52 Mafa-A and 74 Mafa-B haplotypes in this cohort, many of which were restricted to specific sample origins. We also characterized full-length MHC class I transcripts using Pacific Biosciences (PacBio) RS II single-molecule real-time (SMRT) sequencing. This technology allows for complete read-through of unfragmented MHC class I transcripts (~1100 bp in length), so no assembly is required to unambiguously resolve novel full-length sequences. Overall, we identified 311 total full-length transcripts in a subset of 72 cynomolgus macaques from these Chinese breeding facilities; 130 of these sequences were novel and an additional 115 extended existing short database sequences to span the complete open reading frame. This significantly expands the number of Mafa-A, Mafa-B, and Mafa-I full-length alleles in the official cynomolgus macaque MHC class I database. The PacBio technique described here represents a general method for full-length allele discovery and genotyping that can be extended to other complex immune loci such as MHC class II, killer immunoglobulin-like receptors, and Fc gamma receptors.

The orthologs of HLA-DQ and -DP genes display abundant levels of variability in macaque species.

The human major histocompatibility complex (MHC) region encodes three types of class II molecules designated HLA-DR, -DQ, and -DP. Both the HLA-DQ and -DP gene region comprise a duplicated tandem of A and B genes, whereas in macaques, only one set of genes is present per region. A substantial sequencing project on the DQ and DP genes in various macaque populations resulted in the detection of previously 304 unreported full-length alleles. Phylogenetic studies showed that humans and macaques share trans-species lineages for the DQA1 and DQB1 genes, whereas the DPA1 and DPB1 lineages in macaques appear to be species-specific. Amino acid variability plot analyses revealed that each of the four genes displays more allelic variation in macaques than is encountered in humans. Moreover, the numbers of different amino acids at certain positions in the encoded proteins are higher than in humans. This phenomenon is remarkably prominent at the contact positions of the peptide-binding sites of the deduced macaque DPß-chains. These differences in the MHC class II DP regions of macaques and humans suggest separate evolutionary mechanisms in the generation of diversity.

Differential recombination dynamics within the MHC of macaque species.

A panel of 15 carefully selected microsatellites (short tandem repeats, STRs) has allowed us to study segregation and haplotype stability in various macaque species. The STRs span the major histocompatibility complex (MHC) region and map in more detail from the centromeric part of the Mhc-A to the DR region. Two large panels of Indian rhesus and Indonesian/Indochinese cynomolgus macaques have been subjected to pedigree analysis, allowing the definition of 161 and 36 different haplotypes and the physical mapping of 10 and 5 recombination sites, respectively. Although most recombination sites within the studied section of the Indian rhesus monkey MHC are situated between the Mhc-A and Mhc-B regions, the resulting recombination rate for this genomic segment is low and similar to that in humans. In contrast, in Indonesian/Indochinese macaques, two recombination sites, which appear to be absent in rhesus macaques, map between the class III and II regions. As a result, the mean recombination frequency of the core MHC, Mhc-A to class II, is higher in Indonesian/Indochinese cynomolgus than in Indian rhesus macaques, but as such is comparable to that in humans. The present communication demonstrates that the dynamics of recombination 'hot/cold spots' in the MHC, as well as their frequencies, may differ substantially between highly related macaque species.

Evolution of HLA-DRB genes.

The HLA region shows diversity concerning the number and content of DRB genes present per haplotype. Similar observations are made for the equivalent regions in other primate species. To elucidate the evolutionary history of the various HLA-DRB genes, a large panel of intron sequences obtained from humans, chimpanzees, rhesus macaques, and common marmosets has been subjected to phylogenetic analyses. Special attention was paid to the presence and absence of particular transposable elements and/or to their segments. The sharing of different parts of the same long interspersed nuclear element-2 (LINE2, L2) and various Alu insertions by the species studied demonstrates that one precursor gene must have been duplicated several times before the Old World monkey (OWM) and hominid (HOM) divergence. At least four ancestral DRB gene families appear to have been present before the radiation of OWM and HOM, and one of these even predates the speciation of Old and New World primates. Two of these families represent the pseudogenes DRB6/DRB2 and DRB7, which have been locked in the genomes of various primate species over long evolutionary time spans. Furthermore, all phylogenies of different intron segments show consistently that, apart from the pseudogenes, only DRB5 genes are shared by OWM and HOM, and they demonstrate the common history of certain DRB genes/lineages of humans and chimpanzees. In contrast, the evolutionary history of some other DRB loci is difficult to decipher, thus illustrating the complex history of the evolution of DRB genes due to a combination of mutations and recombination-like events. The selected approach allowed us to shed light on the ancestral DRB gene pool in primates and on the evolutionary relationship of the various HLA-DRB genes.

Gluten-sensitive enteropathy coincides with decreased capability of intestinal T cells to secrete IL-17 and IL-22 in a macaque model for celiac disease.

Celiac disease (CD) is an autoimmune disorder caused by intolerance to dietary gluten. The interleukin (IL)-17 and IL-22 function as innate regulators of mucosal integrity. Impaired but not well-understood kinetics of the IL-17/22 secretion was described in celiac patients. Here, the IL-17 and IL-22-producing intestinal cells were studied upon their in vitro stimulation with mitogens in class II major histocompatibility complex-defined, gluten-sensitive rhesus macaques. Pediatric biopsies were collected from distal duodenum during the stages of disease remission and relapse. Regardless of dietary gluten content, IL-17 and IL-22-producing cells consisted of CD4+ and CD8+ T lymphocytes as well as of lineage-negative (Lin-) cells. Upon introduction of dietary gluten, capability of intestinal T cells to secrete IL-17/22 started to decline (p<0.05), which was paralleled with gradual disruption of epithelial integrity. These data indicate that IL-17/22-producing cells play an important role in maintenance of intestinal mucosa in gluten-sensitive primates.

The repertoire of MHC class I genes in the common marmoset: evidence for functional plasticity.

In humans, the classical antigen presentation function of major histocompatibility complex (MHC) class I molecules is controlled by the human leukocyte antigen HLA -A, HLA-B and HLA-C loci. A similar observation has been made for great apes and Old World monkey species. In contrast, a New World monkey species such as the cotton-top tamarin (Saguinus oedipus) appears to employ the G locus for its classical antigen presentation function. At present, little is known about the classical MHC class I repertoire of the common marmoset (Callithrix jacchus), another New World monkey that is widely used in biomedical research. In the present population study, no evidence has been found for abundant transcription of classical I class genes. However, in each common marmoset, four to seven different G-like alleles were detected, suggesting that the ancestral locus has been subject to expansion. Segregation studies provided evidence for at least two G-like genes present per haplotype, which are transcribed by a variety of cell types. The alleles of these Caja-G genes cluster in separate lineages, suggesting that the loci diversified considerably after duplication. Phylogenetic analyses of the introns confirm that the Caja-G loci cluster in the vicinity of HLA-G, indicating that both genes shared an ancestor. In contrast to HLA-G, Caja-G shows considerable polymorphism at the peptide-binding sites. This observation, together with the lack of detectable transcripts of A and B-like genes, indicates that Caja-G genes have taken over the function of classical class I genes. These data highlight the extreme plasticity of the MHC class I gene system.

Unique peptide-binding motif for Mamu-B*037:01: an MHC class I allele common to Indian and Chinese rhesus macaques.

Indian and Chinese rhesus macaques are often used in biomedical research. Genetic analyses of the major histocompatibility class I region have revealed that these macaques display a substantial level of polymorphism at Mamu-A and Mamu-B loci, which have been subject to duplication. Only a few Mamu class I allotypes are characterised for their peptide-binding motifs, although more information of this nature would contribute to a better interpretation of T cell-mediated immune responses. Here, we present the results of the characterisation of the functional properties of Mamu-B*037:01, an allotype commonly encountered in rhesus macaques of Indian and Chinese origin. Mamu-B*037:01 is seen to have a strong preference for acidic amino acids at the third residue, and for arginine, lysine, and tyrosine at the carboxyl terminus. This peptide-binding motif is not described in the human population.

Haplotype diversity generated by ancient recombination-like events in the MHC of Indian rhesus macaques.

The Mamu-A, Mamu-B, and Mamu-DRB genes of the rhesus macaque show several levels of complexity such as allelic heterogeneity (polymorphism), copy number variation, differential segregation of genes/alleles present on a haplotype (diversity) and transcription level differences. A combination of techniques was implemented to screen a large panel of pedigreed Indian rhesus macaques (1,384 individuals representing the offspring of 137 founding animals) for haplotype diversity in an efficient and inexpensive manner. This approach allowed the definition of 140 haplotypes that display a relatively low degree of region variation as reflected by the presence of only 17 A, 18 B and 22 DRB types, respectively, exhibiting a global linkage disequilibrium comparable to that in humans. This finding contrasts with the situation observed in rhesus macaques from other geographic origins and in cynomolgus monkeys from Indonesia. In these latter populations, nearly every haplotype appears to be characterised by a unique A, B and DRB region. In the Indian population, however, a reshuffling of existing segments generated "new" haplotypes. Since the recombination frequency within the core MHC of the Indian rhesus macaques is relatively low, the various haplotypes were most probably produced by recombination events that accumulated over a long evolutionary time span. This idea is in accord with the notion that Indian rhesus macaques experienced a severe reduction in population during the Pleistocene due to a bottleneck caused by geographic changes. Thus, recombination-like processes appear to be a way to expand a diminished genetic repertoire in an isolated and relatively small founder population.

AIDS-protective HLA-B*27/B*57 and chimpanzee MHC class I molecules target analogous conserved areas of HIV-1/SIVcpz.

In the absence of treatment, most HIV-1-infected humans develop AIDS. However, a minority are long-term nonprogressors, and resistance is associated with the presence of particular HLA-B*27/B*57 molecules. In contrast, most HIV-1-infected chimpanzees do not contract AIDS. In comparison with humans, chimpanzees experienced an ancient selective sweep affecting the MHC class I repertoire. We have determined the peptide-binding properties of frequent chimpanzee MHC class I molecules, and show that, like HLA-B*27/B*57, they target similar conserved areas of HIV-1/SIV(cpz). In addition, many animals appear to possess multiple molecules targeting various conserved areas of the HIV-1/SIV(cpz) Gag protein, a quantitative aspect of the immune response that may further minimize the chance of viral escape. The functional characteristics of the contemporary chimpanzee MHC repertoire suggest that the selective sweep was caused by a lentiviral pandemic.

The extreme plasticity of killer cell Ig-like receptor (KIR) haplotypes differentiates rhesus macaques from humans.

NK cells are essential in shaping immune responses and play an important role during pregnancy and in controlling infections. Killer cell immunoglobulin-like receptors (KIRs) educate the NK cell and determine its state of activation. Our goal was to determine how the KIR repertoire of the rhesus macaque (Macaca mulatta) has been shaped during evolution. The presence or absence of 22 KIR gene groups was determined in 378 animals. Some unexpected observations were made in an outbred colony comprising animals of different origins. For instance, the KIR region appears to be highly plastic, and an unprecedented number of genotypes and haplotypes was observed. In contrast to humans, there is no distinction between group A and B haplotypes in the rhesus macaque, suggesting that different selective forces may be operative. Moreover, specific genes appear to be either present or absent in animals of different geographic origins. This extreme plasticity may have been propelled by co-evolution with the rhesus macaque MHC class I region, which shows signatures of expansion. The mosaic-like complexity of KIR genotypes as observed at the population level may represent an effective strategy for surviving epidemic infections.

DR haplotype diversity of the cynomolgus macaque as defined by its transcriptome.

The DR region of particular primate species may display allelic polymorphism and gene copy number variation (region configuration polymorphism). The sum of these distinct types of polymorphism is defined as complexity. To date, however, the DR region of cynomolgus macaques (Macaca fascicularis) has been poorly defined. Transcriptome analysis of a pedigreed colony, comprising animals from Indonesia and Indochina, revealed a total of 15 Mafa-DRA and 57 DRB alleles, specifying 28 different region configurations. The DRA alleles can be divided into two distinct lineages. One lineage is polymorphic, but the majority of the amino acid replacements map to the leader peptide. The second lineage is at best oligomorphic, and segregates with one specific Mafa-DRB allele. The number of Mafa-DRB genes ranges from two to five per haplotype. Due to the presence of pseudogenes, however, each haplotype encodes only one to three bona fide DRB transcripts. Depending on the region configuration in which the Mafa-DRB gene is embedded, identical alleles may display differential transcription levels. Region configurations appear to have been generated by recombination-like events. When genes or gene segments are relocated, it seems plausible that they may be placed in the context of distinct transcription control elements. As such, DRB region-related transcription level differences may add an extra layer of polymorphism to this section of the adaptive immune system.