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

Hypoimmune islets achieve insulin independence after allogeneic transplantation in a fully immunocompetent non-human primate

Allogeneic transplantation of pancreatic islets for patients with difficult-to-control diabetes mellitus is severely hampered by the requirement for continuous immunosuppression and its associated morbidity. We report that allogeneic transplantation of genetically engineered (B2M-/-, CIITA-/-, CD47+), primary, hypoimmune, pseudo-islets (p-islets) results in their engraftment into a fully immunocompetent, diabetic non-human primate wherein they provide stable endocrine function and enable insulin independence without inducing any detectable immune response in the absence of immunosuppression. Hypoimmune primary p-islets may provide a curative cell therapy for type 1 diabetes mellitus.

Long-read assembly of major histocompatibility complex and killer cell immunoglobulin-like receptor genome regions in cynomolgus macaque

BACKGROUND

The major histocompatibility complex (MHC) and the killer cell immunoglobulin-like receptors (KIR) are key regulators of immune responses. The cynomolgus macaque, an Old World monkey species, can be applied as an important preclinical model for studying human diseases, including coronavirus disease 2019 (COVID-19). Several MHC-KIR combinations have been associated with either a poor or good prognosis. Therefore, macaques with a well-characterized immunogenetic profile may improve drug evaluation and speed up vaccine development. At present, a complete overview of the MHC and KIR haplotype organizations in cynomolgus macaques is lacking, and characterization by conventional techniques is hampered by the extensive expansion of the macaque MHC-B region that complicates the discrimination between genes and alleles.

METHODS

We assembled complete MHC and KIR genomic regions of cynomolgus macaque using third-generation long-read sequencing approach. We identified functional Mafa-B loci at the transcriptome level using locus-specific amplification in a cohort of 33 Vietnamese cynomolgus macaques.

RESULTS

This is the first physical mapping of complete MHC and KIR gene regions in a Vietnamese cynomolgus macaque. Furthermore, we identified four functional Mafa-B loci (B2, B3, B5, and B6) and showed that alleles of the Mafa-I*01, -B*056, -B*034, and -B*001 functional lineages, respectively, are highly frequent in the Vietnamese cynomolgus macaque population.

CONCLUSION

The insights into the MHC and KIR haplotype organizations and the level of diversity may refine the selection of animals with specific genetic markers for future medical research.

The Cynomolgus Macaque MHC Polymorphism in Experimental Medicine

Among the non-human primates used in experimental medicine, cynomolgus macaques (Macaca fascicularis hereafter referred to as Mafa) are increasingly selected for the ease with which they are maintained and bred in captivity. Macaques belong to Old World monkeys and are phylogenetically much closer to humans than rodents, which are still the most frequently used animal model. Our understanding of the Mafa genome has progressed rapidly in recent years and has greatly benefited from the latest technical advances in molecular genetics. Cynomolgus macaques are widespread in Southeast Asia and numerous studies have shown a distinct genetic differentiation of continental and island populations. The major histocompatibility complex of cynomolgus macaque (Mafa MHC) is organized in the same way as that of human, but it differs from the latter by its high degree of classical class I gene duplication. Human polymorphic MHC regions play a pivotal role in allograft transplantation and have been associated with more than 100 diseases and/or phenotypes. The Mafa MHC polymorphism similarly plays a crucial role in experimental allografts of organs and stem cells. Experimental results show that the Mafa MHC class I and II regions influence the ability to mount an immune response against infectious pathogens and vaccines. MHC also affects cynomolgus macaque reproduction and impacts on numerous biological parameters. This review describes the Mafa MHC polymorphism and the methods currently used to characterize it. We discuss some of the major areas of experimental medicine where an effect induced by MHC polymorphism has been demonstrated.

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.

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.

Diversification of both KIR and NKG2 natural killer cell receptor genes in macaques - implications for highly complex MHC-dependent regulation of natural killer cells

The killer immunoglobulin-like receptors (KIR) as well as their MHC class I ligands display enormous genetic diversity and polymorphism in macaque species. Signals resulting from interaction between KIR or CD94/NKG2 receptors and their cognate MHC class I proteins essentially regulate the activity of natural killer (NK) cells. Macaque and human KIR share many features, such as clonal expression patterns, gene copy number variations, specificity for particular MHC class I allotypes, or epistasis between KIR and MHC class I genes that influence susceptibility and resistance to immunodeficiency virus infection. In this review article we also annotated publicly available rhesus macaque BAC clone sequences and provide the first description of the CD94-NKG2 genomic region. Besides the presence of genes that are orthologous to human NKG2A and NKG2F, this region contains three NKG2C paralogues. Hence, the genome of rhesus macaques contains moderately expanded and diversified NKG2 genes in addition to highly diversified KIR genes. The presence of two diversified NK cell receptor families in one species has not been described before and is expected to require a complex MHC-dependent regulation of NK cells.

Comparative genomics of the human, macaque and mouse major histocompatibility complex

The MHC is a highly polymorphic genomic region that encodes the transplantation and immune regulatory molecules. It receives special attention for genetic investigation because of its important role in the regulation of innate and adaptive immune responses and its strong association with numerous infectious and/or autoimmune diseases. The MHC locus was first discovered in the mouse and for the past 50 years it has been studied most intensively in both mice and humans. However, in recent years the macaque species have emerged as some of the more important and advanced experimental animal models for biomedical research into MHC with important human immunodeficiency virus/simian immunodeficiency virus and transplantation studies undertaken in association with precise MHC genotyping and haplotyping methods using Sanger sequencing and next-generation sequencing. Here, in this special issue on 'Macaque Immunology' we provide a short review of the genomic similarities and differences among the human, macaque and mouse MHC class I and class II regions, with an emphasis on the association of the macaque class I region with MHC polymorphism, haplotype structure and function.

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.

Phylogenetic relationships of Malaysia's long-tailed macaques, Macaca fascicularis, based on cytochrome b sequences

Phylogenetic relationships among Malaysia's long-tailed macaques have yet to be established, despite abundant genetic studies of the species worldwide. The aims of this study are to examine the phylogenetic relationships of Macaca fascicularis in Malaysia and to test its classification as a morphological subspecies. A total of 25 genetic samples of M. fascicularis yielding 383 bp of Cytochrome b (Cyt b) sequences were used in phylogenetic analysis along with one sample each of M. nemestrina and M. arctoides used as outgroups. Sequence character analysis reveals that Cyt b locus is a highly conserved region with only 23% parsimony informative character detected among ingroups. Further analysis indicates a clear separation between populations originating from different regions; the Malay Peninsula versus Borneo Insular, the East Coast versus West Coast of the Malay Peninsula, and the island versus mainland Malay Peninsula populations. Phylogenetic trees (NJ, MP and Bayesian) portray a consistent clustering paradigm as Borneo's population was distinguished from Peninsula's population (99% and 100% bootstrap value in NJ and MP respectively and 1.00 posterior probability in Bayesian trees). The East coast population was separated from other Peninsula populations (64% in NJ, 66% in MP and 0.53 posterior probability in Bayesian). West coast populations were divided into 2 clades: the North-South (47%/54% in NJ, 26/26% in MP and 1.00/0.80 posterior probability in Bayesian) and Island-Mainland (93% in NJ, 90% in MP and 1.00 posterior probability in Bayesian). The results confirm the previous morphological assignment of 2 subspecies, M. f. fascicularis and M. f. argentimembris, in the Malay Peninsula. These populations should be treated as separate genetic entities in order to conserve the genetic diversity of Malaysia's M. fascicularis. These findings are crucial in aiding the conservation management and translocation process of M. fascicularis populations in Malaysia.

Identification of MHC class I sequences in four species of Macaca of China

Tibetan macaques (Macaca thibetana), stump-tailed macaques (M. arctoides), Assamese macaques (M. assamensis), and northern pig-tailed macaques (M. leonina) are four major species of Macaca in China. In order to effectively use these species in biomedical research, thorough investigations of their MHC immunogenetics are required. In this study, we identified MHC class I sequences using cDNA cloning and sequencing on a cohort of six M. thibetana, three M. arctoides, three M. assamensis, and three M. leonina derived from Sichuan and Yunnan provinces of China. Eighty new alleles were identified, including 26 MHC-A alleles, 46 MHC-B alleles, and 8 MHC-I alleles. Among them, Math-A1*126:01, Math-B*190:01, Math-B*191:01, Math-B*192:01, Maar-A1*127:01, Maar-A1*129:01, and Maas-A1*128:01 represent lineages that had not been reported earlier in Macaca. Phylogenetic analyses show that no obvious separation of lineages among these species of Macaca. This study provides important information about the MHC immunogenetics for the four major species of Chinese macaques and adds value to these species as model organisms in biomedical research.

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.

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.

Expression profile of hepatic genes in cynomolgus macaques bred in Cambodia, China, and Indonesia: implications for cytochrome P450 genes

Cynomolgus macaques, frequently used in drug metabolism studies, are bred mainly in the countries of Asia; however, comparative studies of drug metabolism between cynomolgus macaques bred in these countries have not been conducted. In this study, hepatic gene expression profiles of cynomolgus macaques bred in Cambodia (mfCAM), China (mfCHN), and Indonesia (mfIDN) were analyzed. Microarray analysis revealed that expression of most hepatic genes, including drug-metabolizing enzyme genes, was not substantially different between mfCAM, mfCHN, and mfIDN; only 1.1% and 3.0% of all the gene probes detected differential expression (>2.5-fold) in mfCAM compared with mfCHN and mfIDN, respectively. Quantitative polymerase chain reaction showed that the expression levels of 14 cytochromes P450 (P450s) important for drug metabolism did not differ (>2.5-fold) in mfCAM, mfCHN, and mfIDN, validating the microarray data. In contrast, expression of CYP2B6 and CYP3A4 differed (>2.5-fold, p < 0.05) between cynomolgus (mfCAM, mfCHN, or mfIDN) and rhesus macaques, indicating greater differences in expression of P450 genes between the two lineages. Moreover, metabolic activities measured using 14 P450 substrates did not differ substantially (<1.5-fold) between mfCAM and mfCHN. These results suggest that gene expression profiles, including drug-metabolizing enzyme genes such as P450 genes, are similar in mfCAM, mfCHN, and mfIDN.

Diversity of MHC class I haplotypes in cynomolgus macaques

Cynomolgus macaques are widely used as a primate model for human diseases associated with an immunological process. Because there are individual differences in immune responsiveness, which are controlled by the polymorphic nature of the major histocompatibility (MHC) locus, it is important to reveal the diversity of MHC in the model animal. In this study, we analyzed 26 cynomolgus macaques from five families for MHC class I genes. We identified 32 Mafa-A, 46 Mafa-B, 6 Mafa-I, and 3 Mafa-AG alleles in which 14, 20, 3, and 3 alleles were novel. There were 23 MHC class I haplotypes and each haplotype was composed of one to three Mafa-A alleles and one to five Mafa-B alleles. Family studies revealed that there were two haplotypes which contained two Mafa-A1 alleles. These observations demonstrated further the complexity of MHC class I locus in the Old World monkey.

Characterisation of MHC haplotypes in a breeding colony of Indonesian cynomolgus macaques reveals a high level of diversity

Recent reports have revealed that cynomolgus macaques obtained from different geographic origins may be more or less suitable for particular studies depending on the specific question(s) being addressed, e.g. Mauritian cynomolgus macaques are particularly suitable for detailed immunological studies against a limited genetic background while less conserved populations may be more appropriate to predict breadth of vaccine coverage in the genetically diverse human population. We have characterised MHC haplotypes in 90 Indonesian cynomolgus macaques using microsatellite and reference strand conformational analysis. Thirty unique haplotypes were defined in the cohort, emphasising the high degree of diversity in this population of cynomolgus macaques. The majority of haplotypes were present at a frequency of ≤ 6%. Transcription profiles indicated that each haplotype was associated with two to eight transcribed class I alleles. The results corroborate previous reports of the extensive MHC diversity of Indonesian cynomolgus macaques and provide additional data to inform colony management decisions. Further, definition of the MHC diversity of the population satisfies one of the prerequisites to MHC association studies and detailed immunological investigations in this outbred non-human primate species.

Twenty-three novel major histocompatibility complex class I B alleles identified in cynomolgus macaques of Vietnamese origin

We report herein the identification of 23 novel Mafa-B alleles in cynomolgus macaques of Vietnamese origin.

Immunization with recombinant HLA classes I and II, HIV-1 gp140, and SIV p27 elicits protection against heterologous SHIV infection in rhesus macaques

Major histocompatibility complex (MHC) molecules expressed on the surface of human immunodeficiency virus (HIV) are potential targets for neutralizing antibodies. Since MHC molecules are polymorphic, nonself MHC can also be immunogenic. We have used combinations of novel recombinant HLA class I and II and HIV/simian immunodeficiency virus (SIV) antigens, all linked to dextran, to investigate whether they can elicit protective immunity against heterologous simian/human immunodeficiency virus (SHIV) challenge in rhesus macaques. Three groups of animals were immunized with HLA (group 1, n = 8), trimeric YU2 HIV type 1 (HIV-1) gp140 and SIV p27 (HIV/SIV antigens; group 2, n = 8), or HLA plus HIV/SIV antigens (group 3, n = 8), all with Hsp70 and TiterMax Gold adjuvant. Another group (group 4, n = 6) received the same vaccine as group 3 without TiterMax Gold. Two of eight macaques in group 3 were completely protected against intravenous challenge with 18 50% animal infective doses (AID(50)) of SHIV-SF162P4/C grown in human cells expressing HLA class I and II lineages represented in the vaccine, while the remaining six macaques showed decreased viral loads compared to those in unimmunized animals. Complement-dependent neutralizing activity in serum and high levels of anti-HLA antibodies were elicited in groups 1 and 3, and both were inversely correlated with the plasma viral load at 2 weeks postchallenge. Antibody-mediated protection was strongly supported by the fact that transfer of pooled serum from the two challenged but uninfected animals protected two naïve animals against repeated low-dose challenge with the same SHIV stock. This study demonstrates that immunization with recombinant HLA in combination with HIV-1 antigens might be developed into an alternative strategy for a future AIDS vaccine.

Primate-specific regulation of natural killer cells

Natural killer (NK) cells are circulating lymphocytes that function in innate immunity and placental reproduction. Regulating both development and function of NK cells is an array of variable and conserved receptors that interact with major histocompatibility complex (MHC) class I molecules. Families of lectin-like and immunoglobulin-like receptors are determined by genes in the natural killer complex (NKC) and leukocyte receptor complex (LRC), respectively. As a consequence of the strong, varying pressures on the immune and reproductive systems, NK cell receptors and their MHC class I ligands evolve rapidly, are highly diverse and exhibit dramatic species-specific differences. The variable, polymorphic family of killer cell immunoglobulin-like receptors (KIR) that regulate human NK cell development and function arose recently, from a single-copy gene during the evolution of simian primates. Our studies of KIR and MHC class I genes in representative species show how these two unlinked but functionally intertwined genetic complexes have co-evolved. In humans, combinations of KIR and HLA class I factors are associated with infectious diseases, including HIV/AIDS, autoimmunity, reproductive success and the outcome of therapeutic transplantation. The extraordinary, and unanticipated, divergence of human NK cell receptors and MHC class I ligands from their mouse counterparts can in part explain the difficulties experienced in finding informative mouse models for human diseases. Non-human primate models have far greater potential, but to realize their promise will first require more complete definition of the genetics and function of KIR and MHC variation in non-human primate species, at a level comparable to that achieved for the human species.

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.

Diversity of MHC class I genes in Burmese-origin rhesus macaques

Rhesus macaques (Macaca mulatta) are widely used in developing a strategy for vaccination against human immunodeficiency virus by using simian immunodeficiency virus infection as a model system. Because the genome diversity of major histocompatibility complex (MHC) is well known to control the immune responsiveness to foreign antigens, MHC loci in Indian- and Chinese-origin macaques used in the experiments have been characterized, and it was revealed that the diversity of MHC in macaques was larger than the human MHC. To further characterize the diversity of Mamu-A and Mamu-B loci, we investigated a total of 73 different sequences of Mamu-A, 83 sequences of Mamu-B, and 15 sequences of Mamu-I cDNAs isolated from Burmese-origin macaques. It was found that there were one to five expressing genes in each locus. Among the Mamu-A, Mamu-B, and Mamu-I sequences, 44 (60.2%), 45 (54.2%), and 8 (53.3%), respectively, were novel, and most of the other known alleles were identical to those reported from Chinese- or Indian-origin macaques, demonstrating a genetic mixture between the geographically distinct populations of present day China and India. In addition, it was found that a Mamu haplotype contained at least two highly transcribed Mamu-A genes, because multiple Mamu-A1 cDNAs were obtained from one haplotype. These findings further revealed the diversity and complexity of MHC locus in the rhesus macaques.