Dynamic evolution of Mhc haplotypes in cynomolgus macaques of different geographic origins

The major histocompatibility complex (MHC) plays a key role in immune defense, and the Mhc genes of cynomolgus macaque display a high degree of polymorphism. Based on their geographic distribution, different populations of cynomolgus macaques are recognized. Here we present the characterization of the Mhc class I and II repertoire of a large pedigreed group of cynomolgus macaques originating from the mainland north of the isthmus of Kra (N = 42). Segregation analyses resulted in the definition of 81 unreported Mafa-A/B/DRB/DQ/DP haplotypes, which include 32 previously unknown DRB regions. In addition, we report 13 newly defined Mafa-A/B/DRB/DQ/DP haplotypes in a group of cynomolgus macaques originating from the mainland south of the isthmus of Kra/Maritime Southeast Asia (N = 16). A relatively high level of sharing of Mafa-A (51%) and Mafa-B (40%) lineage groups is observed between the populations native to the north and the south of isthmus of Kra. At the allelic level, however, the Mafa-A/B haplotypes seem to be characteristic of a population. An overall comparison of all currently known data revealed that each geographic population has its own specific combinations of Mhc class I and II haplotypes. This illustrates the dynamic evolution of the cynomolgus macaque Mhc region, which was most likely generated by recombination and maintained by selection due to the differential pathogenic pressures encountered in different geographic areas.

Comparative genetics of KIR haplotype diversity in humans and rhesus macaques: the balancing act

The role of natural killer (NK) cells is tightly modulated by interactions of killer cell immunoglobulin-like receptors (KIR) with their ligands of the MHC class I family. Several characteristics of the KIR gene products are conserved in primate evolution, like the receptor structures and the variegated expression pattern. At the genomic level, however, the clusters encoding the KIR family display species-specific diversity, reflected by differential gene expansions and haplotype architecture. The human KIR cluster is extensively studied in large cohorts from various populations, which revealed two KIR haplotype groups, A and B, that represent more inhibitory and more activating functional profiles, respectively. So far, genomic KIR analyses in large outbred populations of non-human primate species are lacking. In this study, we roughly quadrupled the number of rhesus macaques studied for their KIR transcriptome (n = 298). Using segregation analysis, we defined 112 unique KIR region configurations, half of which display a more inhibitory profile, whereas the other half has a more activating potential. The frequencies and functional potential of these profiles might mirror the human KIR haplotype groups. However, whereas the human group A and B KIR haplotypes are confined to largely fixed organizations, the haplotypes in macaques feature highly variable gene content. Moreover, KIR homozygosity was hardly encountered in this panel of macaques. This study exhibits highly diverse haplotype architectures in humans and macaques, which nevertheless might have an equivalent effect on the modulation of NK cell activity.

Unparalleled Rapid Evolution of KIR Genes in Rhesus and Cynomolgus Macaque Populations

The killer cell Ig-like receptors (KIR) modulate immune responses through interactions with MHC class I molecules. The KIR region in large cohorts of rhesus and cynomolgus macaque populations were characterized, and the experimental design enabled the definition of a considerable number of alleles (n = 576) and haplotypes, which are highly variable with regard to architecture. Although high levels of polymorphism were recorded, only a few alleles are shared between species and populations. The rapid evolution of allelic polymorphism, accumulated by point mutations, was further confirmed by the emergence of a novel KIR allele in a rhesus macaque family. In addition to allelic variation, abundant orthologous and species-specific KIR genes were identified, the latter of which are frequently generated by fusion events. The concerted action of both genetic mechanisms, in combination with differential selective pressures at the population level, resulted in the unparalleled rapid evolution of the KIR gene region in two closely related macaque species. The variation of the KIR gene repertoire at the species and population level might have an impact on the outcome of preclinical studies with macaque models.

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.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

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.

Spontaneous endometriosis in rhesus macaques: evidence for a genetic association with specific Mamu-A1 alleles

Endometriosis is a poorly understood common debilitating women's reproductive disorder resulting from proliferative and ectopic endometrial tissue associated with variable clinical symptoms including dysmenorrhea (painful menstrual periods), dyspareunia (pain on intercourse), female infertility, and an increased risk of malignant transformation. The rhesus macaque (Macaca mulatta) develops a spontaneous endometriosis that is very similar to that seen in women. We hypothesized that specific major histocompatibility complex (MHC) alleles may contribute to the pathogenesis of endometriosis. As part of a collaboration between the Biomedical Primate Research Centre (BPRC) in the Netherlands and the New England Primate Research Center (NEPRC) in the United States, we analyzed DNA sequences of MHC class I (Macaca mulatta, Mamu-A1) and class II (Mamu-DRB) alleles from rhesus macaques with endometriosis and compared the allele frequencies with those of age-matched healthy macaques. We demonstrate that two MHC class I alleles are overrepresented in diseased macaques compared to controls: Mamu-A1*001, 33.3 % in BPRC animals with endometriosis vs. 11.6 % in healthy macaques (p= 0.007), and Mamu-A1*007, 21.9 % NEPRC rhesus macaques vs. 6.7 %, (p= 0.003). We provide evidence that select MHC class I alleles are associated with endometriosis in rhesus macaques and suggest that the disease pathogenesis contribution of MHC class I warrants further research.

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.

Multilocus definition of MHC haplotypes in pedigreed cynomolgus macaques (Macaca fascicularis)

Cynomolgus macaques (Macaca fascicularis) are used widely in biomedical research, and the genetics of their MHC (Mhc-Mafa) has become the focus of considerable attention in recent years. The cohort of Indonesian pedigreed macaques that we present here was typed for Mafa-A, -B, and -DR, by sequencing, as described in earlier studies. Additionally, the DRB region of these animals was characterised by microsatellite analyses. In this study, full-length sequencing of Mafa-DPA/B and -DQA/B in these animals was performed. A total of 75 different alleles were observed; 22 of which have not previously been reported, plus 18 extended exon 2 alleles that were already known. Furthermore, two microsatellites, D6S2854 and D6S2859, were used to characterise the complex Mafa-A region. Sequencing and segregation analyses revealed that the length patterns of these microsatellites are unique for each Mafa-A haplotype. In this work, we present a pedigreed colony of approximately 120 cynomolgus macaques; all of which are typed for the most significant polymorphic MHC class I and class II markers. Offspring of these pedigreed animals are easily characterised for their MHC by microsatellite analyses on the Mafa-A and -DRB regions, which makes the cumbersome sequencing analyses redundant.

Extensive DRB region diversity in cynomolgus macaques: recombination as a driving force

The DR region of primate species is generally complex and displays diversity concerning the number and combination of distinct types of DRB genes present per region configuration. A highly variable short tandem repeat (STR) present in intron 2 of nearly all primate DRB genes can be utilized as a quick and accurate high through-put typing procedure. This approach resulted previously in the description of unique and haplotype-specific DRB-STR length patterns in humans, chimpanzees, and rhesus macaques. For the present study, a cohort of 230 cynomolgus monkeys, including self-sustaining breeding groups, has been examined. MtDNA analysis showed that most animals originated from the Indonesian islands, but some are derived from the mainland, south and north of the Isthmus of Kra. Haplotyping and subsequent sequencing resulted in the detection of 118 alleles, including 28 unreported ones. A total of 49 Mafa-DRB region configurations were detected, of which 28 have not yet been described. Humans and chimpanzees possess a low number of different DRB region configurations in concert with a high degree of allelic variation. In contrast, however, allelic heterogeneity within a given Mafa-DRB configuration is even less frequently observed than in rhesus macaques. Several of these region configurations appear to have been generated by recombination-like events, most probably propagated by a retroviral element mapping within DRB6 pseudogenes, which are present on the majority of haplotypes. This undocumented high level of DRB region configuration-associated diversity most likely represents a species-specific strategy to cope with various pathogens.

Reshuffling of ancient peptide binding motifs between HLA-DRB multigene family members: old wine served in new skins

In most primate species, the class II region of the Major Histocompatibility Complex (MHC) displays diversity with regard to gene copy number and combination of DRB genes present per region configuration. Some of these loci exhibit extremely high levels of allelic variability, whereas others display only moderate levels of polymorphism. To understand the evolutionary history of the various HLA-DR region genes, a large number of full-length sequences of rhesus macaques, chimpanzees and humans were determined. The exon-intron organisation of the DRA gene, displaying only low levels of polymorphism, appears to have been highly conserved during primate evolution. The physical length of various DRB genes/alleles, however, fluctuates significantly in primates due to the presence of indels (insertions/deletions), mainly mapping to intron 1. Phylogenetic evidence supports the notion that the generation of new DRB genes is a dynamic and steadily ongoing process. Indeed, most of the primate DRB alleles investigated represent relatively young entities, possessing species-unique sequences. This seems to contradict the current view that the highly similar peptide binding motifs of many HLA-, Patr- and Mamu-DR molecules, encoded by exon 2 of the DRB gene, represent old entities, which predate primate speciation. As no evidence was found for convergent evolution, the combination of these two observations indicates that ancient peptide binding motifs are frequently reshuffled among duplicated members of the HLA-DRB multigene family.

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.

A highly divergent microsatellite facilitating fast and accurate DRB haplotyping in humans and rhesus macaques

The DRB region of the MHC in primate species is known to display abundant region configuration polymorphism with regard to the number and content of genes present per haplotype. Furthermore, depending on the species studied, the different DRB genes themselves may display varying degrees of allelic polymorphism. Because of this combination of diversity (differential gene number) and polymorphism (allelic variation), molecular typing methods for the primate DRB region are cumbersome. All intact DRB genes present in humans and rhesus macaques appear to possess, however, a complex and highly divergent microsatellite. Microsatellite analysis of a sizeable panel of outbred rhesus macaques, covering most of the known Mamu-DRB haplotypes, resulted in the definition of unique genotyping patterns that appear to be specific for a given haplotype. Subsequent examination of a representative panel of human cells illustrated that this approach also facilitates high-resolution HLA-DRB typing in an easy, quick, and reproducible fashion. The genetic composition of this complex microsatellite is shown to be in concordance with the phylogenetic relationships of various HLA-DRB and Mamu-DRB exon 2 gene/lineage sequences. Moreover, its length variability segregates with allelic variation of the respective gene. This simple protocol may find application in a variety of research avenues such as transplantation biology, disease association studies, molecular ecology, paternity testing, and forensic medicine.

Microsatellite typing of the rhesus macaque MHC region

To improve the results gained by serotyping rhesus macaque major histocompatibility complex (MHC) antigens, molecular typing techniques have been established for class I and II genes. Like the rhesus macaque Mamu-DRB loci, the Mamu-A and -B are not only polymorphic but also polygenic. As a consequence, sequence-based typing of these genes is time-consuming. Therefore, eight MHC-linked microsatellites, or short tandem repeats (STRs), were evaluated for their use in haplotype characterization. Polymorphism analyses in rhesus macaques of Indian and Chinese origin showed high STR allelic diversity in both populations but different patterns of allele frequency distribution between the groups. Pedigree data for class I and II loci and the eight STRs allowed us to determine extended MHC haplotypes in rhesus macaque breeding groups. STR sequencing and comparisons with the complete rhesus macaque MHC genomic map allowed the exact positioning of the markers. Strong linkage disequilibria were observed between Mamu-DR and -DQ loci and adjacent STRs. Microsatellite typing provides an efficient, robust, and quick method of genotyping and deriving MHC haplotypes for rhesus macaques regardless of their geographical origin. The incorporation of MHC-linked STRs into routine genetic tests will contribute to efforts to improve the genetic characterization of the rhesus macaque for biomedical research and can provide comparative information about the evolution of the MHC region.

Genetic makeup of the DR region in rhesus macaques: gene content, transcripts, and pseudogenes

In the human population, five major HLA-DRB haplotypes have been identified, whereas the situation in rhesus macaques (Macaca mulatta) is radically different. At least 30 Mamu-DRB region configurations, displaying polymorphism with regard to number and combination of DRB loci present per haplotype, have been characterized. Until now, Mamu-DRB region genes have been studied mainly by genomic sequencing of polymorphic exon 2 segments. However, relatively little is known about the expression status of these genes. To understand which exon 2 segments may represent functional genes, full-length cDNA analyses of -DRA and -DRB were initiated. In the course of the study, 11 cDRA alleles were identified, representing four distinct gene products. Amino acid replacements are confined to the leader peptide and cytoplasmatic tail, whereas residues of the alpha1 domain involved in peptide binding, are conserved between humans, chimpanzees, and rhesus macaques. Furthermore, from the 11 Mamu-DRB region configurations present in this panel, 28 cDRB alleles were isolated, constituting 12 distinct cDRA/cDRB configurations. Evidence is presented that a single configuration expresses maximally up to three -DRB genes. For some exon 2 DRB sequences, the corresponding transcripts could not be detected, rendering such alleles as probable pseudogenes. The full-length cDRA and cDRB sequences are necessary to construct Mhc class II tetramers, as well as transfectant cell lines. As the rhesus macaque is an important animal model in AIDS vaccine studies, the information provided in this communication is essential to define restriction elements and to monitor immune responses in SIV/simian human immunodeficiency virus-infected rhesus macaques.

Evolutionary stability of MHC class II haplotypes in diverse rhesus macaque populations

A thoroughly characterized breeding colony of 172 pedigreed rhesus macaques was used to analyze exon 2 of the polymorphic Mamu- DPB1, -DQA1, -DQB1, and - DRB loci. Most of the monkeys or their ancestors originated in India, though the panel also included animals from Burma and China, as well as some of unknown origin and mixed breeds. In these animals, mtDNA appears to correlate with the aforementioned geographic origin, and a large number of Mamu class II alleles were observed. The different Mamu- DPB1 alleles were largely shared between monkeys of different origin, whereas in humans particular alleles appear to be unique for ethnic populations. In contrast to Mamu-DPB1, the highly polymorphic - DQA1/DQB1 alleles form tightly linked pairs that appear to be about two-thirds population specific. For most of the DQA1/DQB1 pairs, Mamu- DRB region configurations present on the same chromosome have been ascertained, resulting in 41 different -DQ/DRB haplotypes. These distinct DQ/DRB haplotypes seem to be specific for monkeys of a determined origin. Thus, in evolutionary terms, the Mamu-DP, -DQ, and -DR regions show increasing instability with regard to allelic polymorphism, such as for -DP/DQ, or gene content and allelic polymorphism, such as for -DR, resulting in population-specific class II haplotypes. Furthermore, novel haplotypes are generated by recombination-like events. The results imply that mtDNA analysis in combination with Mhc typing is a helpful tool for selecting animals for biomedical experiments.

Poly (U, s 4-U) as a synthetic messenger RNA for polyphenylalanine synthesis in a cell-free system derivem derived from rat liver

Two different U and S4-U containing polymers with U:s4U ratio of 3:1 and 5:1 were synthesized. Their activity as messenger RNA was studied in an amino acid incorporation cell-free system from rat liver. It was shown that both copolymers can stimulate the incorporation of phenylalanine into oligophynylalamyl-tRNA.