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

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.

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 the major histocompatibility complex in humans and nonhuman primates

The major histocompatibility complex (MHC) is one of the most gene-dense regions of the mammalian genome. Multiple genes within the human MHC (HLA) show extensive polymorphism, and currently, more than 26,000 alleles divided over 39 different genes are known. Nonhuman primate (NHP) species are grouped into great and lesser apes and Old and New World monkeys, and their MHC is studied mostly because of their important role as animal models in preclinical research or in connection with conservation biology purposes. The evolutionary equivalents of many of the HLA genes are present in NHP species, and these genes may also show abundant levels of polymorphism. This review is intended to provide a comprehensive comparison relating to the organization and polymorphism of human and NHP MHC regions.

Characterization of 100 extended major histocompatibility complex haplotypes in Indonesian cynomolgus macaques

Many medical advancements—including improvements to anti-rejection therapies in transplantation and vaccine development—rely on preclinical studies conducted in cynomolgus macaques (Macaca fascicularis). Major histocompatibility complex (MHC) class I and class II genes of cynomolgus macaques are orthologous to human leukocyte antigen complex (HLA) class I and class II genes, respectively. Both encode cell-surface proteins involved in cell recognition and rejection of non-host tissues. MHC class I and class II genes are highly polymorphic, so comprehensive genotyping requires the development of complete databases of allelic variants. Our group used PacBio circular consensus sequencing of full-length cDNA amplicons to characterize MHC class I and class II transcript sequences for a cohort of 293 Indonesian cynomolgus macaques (ICM) in a large, pedigreed breeding colony. These studies allowed us to expand the existing database of Macaca fascicularis (Mafa) alleles by identifying an additional 141 MHC class I and 61 class II transcript sequences. In addition, we defined co-segregating combinations of allelic variants as regional haplotypes for 70 Mafa-A, 78 Mafa-B, and 45 Mafa-DRB gene clusters. Finally, we defined class I and class II transcripts that are associated with 100 extended MHC haplotypes in this breeding colony by combining our genotyping analyses with short tandem repeat (STR) patterns across the MHC region. Our sequencing analyses and haplotype definitions improve the utility of these ICM for transplantation studies as well as infectious disease and vaccine research.

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.

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.

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.

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.

Study of MHC class II region polymorphism in the Filipino cynomolgus macaque population

The cynomolgus macaque (Macaca fascicularis) is currently used as an animal model in various fields of immunology especially in the development of innovative vaccines for the prevention and treatment of infectious diseases. The polymorphism of the major histocompatibility complex (MHC) influences the development of adaptive immune responses and it is crucial to characterize the polymorphism of cynomolgus MHC genes. We present here a systematic study of the MHC class II haplotypes in the Filipino macaque population. By the study of a large sample of Filipino animals (N = 353), we have characterized 18 MHC class II haplotypes by means of genotyping seven microsatellites. The animals were DRB genotyped by means of PCR-SSO or DGGE-sequencing on genomic amplified fragments. We cloned and sequenced the complementary DNA (cDNA) of DQA, DQB, DPA, and DPB genes of 117 animals. Combining the microsatellite genotyping and cDNA characterized in the 117 animals, we defined genetic association between the cDNA and the microsatellites and characterized 18 MHC class II haplotypes. For 104 animals out of the 353 studied, the presence of a recombinant haplotype was highly probable. Thirty-four percent of recombination was located in 256 kb segment between D6S2876 and D6S2747 microsatellites, a region encompassing several hot spots of recombination in the human MHC.

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.