Genomic plasticity of the immune-related Mhc class I B region in macaque species

Complete Genomic Assembly of Mauritian Cynomolgus Macaque Killer Ig-like Receptor and Natural Killer Group 2 Haplotypes

Mauritian-origin cynomolgus macaques (MCMs) serve as a powerful nonhuman primate model in biomedical research due to their unique genetic homogeneity, which simplifies experimental designs. Despite their extensive use, a comprehensive understanding of crucial immune-regulating gene families, particularly killer Ig-like receptors (KIR) and NK group 2 (NKG2), has been hindered by the lack of detailed genomic reference assemblies. In this study, we employ advanced long-read sequencing techniques to completely assemble eight KIR and seven NKG2 genomic haplotypes, providing an extensive insight into the structural and allelic diversity of these immunoregulatory gene clusters. Leveraging these genomic resources, we prototype a strategy for genotyping KIR and NKG2 using short-read, whole-exome capture data, illustrating the potential for cost-effective multilocus genotyping at colony scale. These results mark a significant enhancement for biomedical research in MCMs and underscore the feasibility of broad-scale genetic investigations.

Can extreme MHC class I diversity be a feature of a wide geographic range? The example of Seba's short-tailed bat (Carollia perspicillata)

The major histocompatibility complex (MHC) is one of the most diverse genetic regions under pathogen-driven selection because of its central role in antigen binding and immunity. The highest MHC variability, both in terms of the number of individual alleles and gene copies, has so far been found in passerine birds; this is probably attributable to passerine adaptation to both a wide geographic range and a diverse array of habitats. If extraordinary high MHC variation and duplication rates are adaptive features under selection during the evolution of ecologically and taxonomically diverse species, then similarly diverse MHC architectures should be found in bats. Bats are an extremely species-rich mammalian group that is globally widely distributed. Many bat species roost in multitudinous groups and have high contact rates with pathogens, conspecifics, and allospecifics. We have characterized the MHC class I diversity in 116 Panamanian Seba's short-tailed bats (Carollia perspicillata), a widely distributed, generalist, neotropical species. We have detected a remarkable individual and population-level diversity of MHC class I genes, with between seven and 22 alleles and a unique genotype in each individual. This diversity is comparable with that reported in passerine birds and, in both taxonomic groups, further variability has evolved through length polymorphisms. Our findings support the hypothesis that, for species with a geographically broader range, high MHC class I variability is particularly adaptive. Investigation of the details of the underlying adaptive processes and the role of the high MHC diversity in pathogen resistance are important next steps for a better understanding of the role of bats in viral evolution and as carriers of several deadly zoonotic viruses.

Determination of major histocompatibility class I and class II genetic composition of the Caribbean Primate Center specific pathogen-free rhesus macaque (Macaca mulatta) colony based on massively parallel sequencing

BACKGROUND

Knowledge of major histocompatibility complex (MHC) composition and distribution in rhesus macaque colonies is critical for management strategies that maximize the utility of this model for biomedical research.

METHODS

Variation within the Mamu-A and Mamu-B (class I) and DRB, DQA/B, and DPA/B (class II) regions of 379 animals from the Caribbean Primate Research Center's (CPRC) specific pathogen free (SPF) colony was examined using massively parallel sequencing.

RESULTS

Analyses of the 7 MHC loci revealed a background of Indian origin with high levels of variation despite past genetic bottlenecks. All loci exhibited mutual linkage disequilibria while conforming to Hardy-Weinberg expectations suggesting the achievement of mutation-selection balance.

CONCLUSION

The CPRC's SPF colony is a significant resource for research on AIDS and other infectious agents. Characterizing colony-wide MHC variability facilitates the breeding and selection of animals bearing desired haplotypes and increases the investigator's ability to understand the immune responses mounted by these animals.

Genetics of Natural Killer Cells in Human Health, Disease, and Survival

Natural killer (NK) cells have vital functions in human immunity and reproduction. In the innate and adaptive immune responses to infection, particularly by viruses, NK cells respond by secreting inflammatory cytokines and killing infected cells. In reproduction, NK cells are critical for genesis of the placenta, the organ that controls the supply of oxygen and nutrients to the growing fetus. Controlling NK cell functions are interactions of HLA class I with inhibitory NK cell receptors. First evolved was the conserved interaction of HLA-E with CD94:NKG2A; later established were diverse interactions of HLA-A, -B, and -C with killer cell immunoglobulin-like receptors. Characterizing the latter interactions is rapid evolution, which distinguishes human populations and all species of higher primate. Driving this evolution are the different and competing selections imposed by pathogens on NK cell-mediated immunity and by the constraints of human reproduction on NK cell-mediated placentation. Promoting rapid evolution is independent segregation of polymorphic receptors and ligands throughout human populations.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

MHC and KIR Polymorphisms in Rhesus Macaque SIV Infection

Natural killer lymphocytes are essentially involved as the first line of defense against agents such as viruses and malignant cells. The activity of these cells is regulated via interaction of specific and diverse killer cell immunoglobulin-like receptors (KIR) with the highly polymorphic cognate MHC class I proteins on target cells. Genetic variability of both KIR and MHC-I ligands has been shown to be associated with resistance to many diseases, including infection with the immunodeficiency virus. Disease course and progression to AIDS after infection with human immunodeficiency virus-1 (HIV-1) is essentially influenced by the presence of the stimulatory KIR3DS1 receptor in combination with HLA-Bw4. Knowledge of such genetic interactions that contribute to not only disease resistance but also susceptibility are just as important. Such combined genetic factors were recently reported in the rhesus macaque AIDS model. Here, we review the rhesus macaque MHC class I and KIR gene systems and the role of their polymorphisms in the SIV infection model.

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.

Mamu-B genes and their allelic repertoires in different populations of Chinese-origin rhesus macaques

Since rhesus monkeys of Chinese origin have gained greater utilization in recent years, it is urgent to investigate the major histocompatibility complex (MHC) immunogenetics of Chinese rhesus macaques. In this study, we identified 81 Mamu-B sequences using complementary DNA cloning and sequencing on a cohort of 58 rhesus monkeys derived from three local populations of China. Twenty of these Mamu-B alleles are novel and four of them represent new lineages. Although more alleles are shared among different populations than Mamu-A locus, the Mamu-B allelic repertoires found in these three populations of Chinese macaques are largely independent, which underscores the MHC polymorphism among different populations of Chinese rhesus macaques. Our results are an important addition to the limited MHC immunogenetic information available for rhesus macaques of Chinese origin.

Nonhuman primate transplant models finally evolve: detailed immunogenetic analysis creates new models and strengthens the old

Nonhuman primate (NHP) models play a critical role in the translation of novel therapies for transplantation to the clinic. However, although MHC disparity significantly affects the outcome of transplantation, until recently, experiments using NHP models were performed without the ability to rigorously control the degree of MHC disparity in transplant cohorts. In this review, we discuss several key technical breakthroughs in the field, which have finally enabled detailed immunogenetic data to be incorporated into NHP transplantation studies. These advances have created a new gold-standard for NHP transplantation research, which incorporates detailed information regarding the degree of relatedness and the degree of MHC haplotype disparity between transplant pairs and the precise MHC alleles that both donors and recipients express. The adoption of this new standard promises to increase the rigor of NHP transplantation studies and to ensure that these experiments are optimally translatable to patient care.

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.

Genomic plasticity of the MHC class I A region in rhesus macaques: extensive haplotype diversity at the population level as revealed by microsatellites

The Mamu-A genes of the rhesus macaque show different degrees of polymorphism, transcription level variation, and differential haplotype distribution. Per haplotype, usually one “major” transcribed gene is present, A1 (A7), in various combinations with “minor” genes, A2 to A6. In silico analysis of the physical map of a heterozygous animal revealed the presence of similar Mamu-A regions consisting of four duplication units, but with dissimilar positions of the A1 genes on both haplotypes, and in combination with different minor genes. Two microsatellites, D6S2854 and D6S2859, have been selected as potential tools to characterize this complex region. Subsequent analysis of a large breeding colony resulted in the description of highly discriminative patterns, displaying copy number variation in concert with microsatellite repeat length differences. Sequencing and segregation analyses revealed that these patterns are unique for each Mamu-A haplotype. In animals of Indian, Burmese, and Chinese origin, 19, 15, or 9 haplotypes, respectively, could be defined, illustrating the occurrence of differential block duplications and subsequent rearrangements by recombination. The haplotypes can be assigned to 12 unique combinations of genes (region configurations). Although most configurations harbor two transcribed A genes, one or three genes per haplotype are also present. Additionally, haplotypes lacking an A1 gene or with an A1 duplication appear to exist. The presence of different transcribed A genes/alleles in monkeys from various origins may have an impact on differential disease susceptibilities. The high-throughput microsatellite technique will be a valuable tool in animal selection for diverse biomedical research projects.

Erratum to: rhesus macaque MHC class I molecules show differential subcellular localizations

The MHC class I gene family of rhesus macaques is characterised by considerable gene duplications. While a HLA-C-orthologous gene is absent, the Mamu-A and in particular the Mamu-B genes have expanded, giving rise to plastic haplotypes with differential gene content. Although some of the rhesus macaque MHC class I genes are known to be associated with susceptibility/resistance to infectious diseases, the functional significance of duplicated Mamu-A and Mamu-B genes and the expression pattern of their encoded proteins are largely unknown. Here, we present data of the subcellular localization of AcGFP-tagged Mamu-A and Mamu-B molecules. We found strong cell surface and low intracellular expression for Mamu-A1, Mamu-A2 and Mamu-A3-encoded molecules as well as for Mamu-B*01704, Mamu-B*02101, Mamu-B*04801, Mamu-B*06002 and Mamu-B*13401. In contrast, weak cell surface and strong intracellular expression was seen for Mamu-A4*1403, Mamu-B*01202, Mamu-B*02804, Mamu-B*03002, Mamu-B*05704, Mamu-I*010201 and Mamu-I*0121. The different expression patterns were assigned to the antigen-binding alpha1 and alpha2 domains, suggesting failure of peptide binding is responsible for retaining 'intracellular' Mamu class I molecules in the endoplasmic reticulum. These findings indicate a diverse functional role of the duplicated rhesus macaque MHC class I genes.

Rhesus macaque MHC class I molecules show differential subcellular localizations

The MHC class I gene family of rhesus macaques is characterised by considerable gene duplications. While a HLA-C-orthologous gene is absent, the Mamu-A and in particular the Mamu-B genes have expanded, giving rise to plastic haplotypes with differential gene content. Although some of the rhesus macaque MHC class I genes are known to be associated with susceptibility/resistance to infectious diseases, the functional significance of duplicated Mamu-A and Mamu-B genes and the expression pattern of their encoded proteins are largely unknown. Here, we present data of the subcellular localization of AcGFP-tagged Mamu-A and Mamu-B molecules. We found strong cell surface and low intracellular expression for Mamu-A1, Mamu-A2 and Mamu-A3-encoded molecules as well as for Mamu-B*01704, Mamu-B*02101, Mamu-B*04801, Mamu-B*06002 and Mamu-B*13401. In contrast, weak cell surface and strong intracellular expression was seen for Mamu-A4*1403, Mamu-B*01202, Mamu-B*02804, Mamu-B*03002, Mamu-B*05704, Mamu-I*010201 and Mamu-I*0121. The different expression patterns were assigned to the antigen-binding α1 and α2 domains, suggesting failure of peptide binding is responsible for retaining ‘intracellular’ Mamu class I molecules in the endoplasmic reticulum. These findings indicate a diverse functional role of the duplicated rhesus macaque MHC class I genes.

Definition of Mafa-A and -B haplotypes in pedigreed cynomolgus macaques (Macaca fascicularis)

The major histocompatibility complex (MHC) class I B gene/allelic repertoire was investigated in a pedigreed population of cynomolgus macaques of mixed Indonesian/Malaysian origin. The Mafa-B alleles detected in this cohort are mostly specific for a given geographic area, and only a small number of alleles appears to be shared with other populations. This suggests the fast evolution of Mafa-B alleles due to adaptation to new environments. In contrast to humans, the B locus in Old World monkeys displays extensive copy number variation. The Mafa-B and previously defined -A gene combinations segregate in families and thus allowed the definition of extended haplotypes. In many cases it was possible to assign a particular Mafa-I allele to one of these Mafa-A/B haplotypes as well. The presence of a large number of stable haplotypes in this cohort of animals, which was pedigreed for up to eight generations, looks promising for developing discriminative MHC typing tools that are less cumbersome. Furthermore, the discovery of 53 unreported Mafa-B sequences expands the lexicon of alleles significantly, and may help in understanding the complex organisation of the macaque B region.

Compound evolutionary history of the rhesus macaque MHC class I B region revealed by microsatellite analysis and localization of retroviral sequences

In humans, the single polymorphic B locus of the major histocompatibility complex is linked to the microsatellite MIB. In rhesus macaques, however, haplotypes are characterized by the presence of unique combinations of multiple B genes, which may display different levels of polymorphism. The aim of the study was to shed light on the evolutionary history of this highly complex region. First, the robustness of the microsatellite MIB-linked to almost half of the B genes in rhesus macaques (Mamu-B)–for accurate B haplotyping was studied. Based on the physical map of an established haplotype comprising 7 MIB loci, each located next to a certain Mamu-B gene, two MIB loci, MIB1 and MIB6, were investigated in a panel of MHC homozygous monkeys. MIB1 revealed a complex genotyping pattern, whereas MIB6 analysis resulted in the detection of one or no amplicon. Both patterns are specific for a given B haplotype, show Mendelian segregation, and even allow a more precise haplotype definition than do traditional typing methods. Second, a search was performed for retroelements that may have played a role in duplication processes as observed in the macaque B region. This resulted in the description of two types of duplicons. One basic unit comprises an expressed Mamu-B gene, adjacent to an HERV16 copy closely linked to MIB. The second type of duplicon comprises a Mamu-B (pseudo)gene, linked to a truncated HERV16 structure lacking its MIB segment. Such truncation seems to coincide with the loss of B gene transcription. Subsequent to the duplication processes, recombination between MIB and Mamu-B loci appears to have occurred, resulting in a hyperplastic B region. Thus, analysis of MIB in addition to B loci allows deciphering of the compound evolutionary history of the class I B region in Old World monkeys.