Bonobos Maintain Immune System Diversity with Three Functional Types of MHC-B

MHC-B Diversity and Signs of Respiratory Illness in Wild, East African Chimpanzees ( Pan troglodytes schweinfurthii )

Many traits, intrinsic and extrinsic to an organism, contribute to interindividual variation in immunity in wild habitats. The vertebrate Major Histocompatibility Complex (MHC) includes genes encoding antigen-presenting molecules that are highly variable, and that variation often predicts susceptibility/resistance to and recovery from pathogen infection. I compare MHC-B variation at two long-term chimpanzee research sites, Kibale National Park in Uganda and Gombe National Park in Tanzania. Using decades of respiratory health data available for these chimpanzees, I test hypotheses associated with maintenance of diversity at MHC loci, including heterozygote, divergent allele, and rare allele advantage hypotheses, and predictions for unique function of MHC-B in great apes. I found, despite confirmation of recent shared ancestry between Kibale and Gombe chimpanzees, including an overlapping MHC-B allele repertoire and similar MHC-B phenotype compositions, chimpanzees from the two research sites experienced differences in the occurrence of respiratory signs and had different associations of MHC-B diversity with signs of respiratory illness. Kibale chimpanzees with heterozygous genotypes and different peptide-binding supertypes were observed less often with respiratory signs than those homozygous or possessing the same supertypes, but this same association was not observed among Gombe chimpanzees. Gombe chimpanzees with specific MHC-B phenotypes that enable engagement of Natural Killer (NK) cells were observed more often with respiratory signs than chimpanzees with other phenotypes, but this was not observed at Kanyawara. This study emphasizes local adaptation in shaping genetic and phenotypic traits in different infectious disease contexts, even among close genetic relatives of the same subspecies, and highlights utility for continued and simultaneous tracking of host immune genes and specific pathogens in wild species.

Malaria-driven adaptation of MHC class I in wild bonobo populations

The malaria parasite Plasmodium falciparum causes substantial human mortality, primarily in equatorial Africa. Enriched in affected African populations, the B*53 variant of HLA-B, a cell surface protein that presents peptide antigens to cytotoxic lymphocytes, confers protection against severe malaria. Gorilla, chimpanzee, and bonobo are humans' closest living relatives. These African apes have HLA-B orthologs and are infected by parasites in the same subgenus (Laverania) as P. falciparum, but the consequences of these infections are unclear. Laverania parasites infect bonobos (Pan paniscus) at only one (TL2) of many sites sampled across their range. TL2 spans the Lomami River and has genetically divergent subpopulations of bonobos on each side. Papa-B, the bonobo ortholog of HLA-B, includes variants having a B*53-like (B07) peptide-binding supertype profile. Here we show that B07 Papa-B occur at high frequency in TL2 bonobos and that malaria appears to have independently selected for different B07 alleles in the two subpopulations.

Cell-Mediated Immune Ontogeny Is Affected by Sex but Not Environmental Context in a Long-Lived Primate Species

Ecoimmunology conceptualizes the role of immunity in shaping life history in a natural context. Within ecoimmunology, macroimmunology is a framework that explains the effects of habitat and spatial differences on variation in immune phenotypes across populations. Within these frameworks, immune ontogeny—the development of the immune system across an individual life span—has received little attention. Here, we investigated how immune ontogeny from birth until adulthood is affected by age, sex, and developmental environment in a long-lived primate species, the bonobo. We found a progressive, significant decline of urinary neopterin levels, a marker for the cell-mediated immune response, from birth until 5 years of age in both sexes. The overall pattern of age-related neopterin changes was sex-specific, with males having higher urinary neopterin levels than females in the first 3 years of life, and females having higher levels than males between 6 and 8 years. Environmental condition (zoo-housed vs. wild) did not influence neopterin levels, nor did age-related changes in neopterin levels differ between environments. Our data suggest that the post-natal development of cell-mediated immune ontogeny is sex-specific but does not show plasticity in response to environmental conditions in this long-lived primate species. This indicates that cell-mediated immune ontogeny in the bonobo follows a stereotypic and maybe a genetically determined pattern that is not affected by environmental differences in pathogen exposure and energy availability, but that sex is an important, yet often overlooked factor shaping patterns of immune ontogeny. Investigating the causes and consequences of variation in immunity throughout life is critical for our understanding of life-history evolution and strategies, mechanisms of sexual selection, and population dynamics with respect to pathogen susceptibility. A general description of sex-specific immune ontogeny as done here is a crucial step in this direction, particularly when it is considered in the context of a species’ ecology and evolutionary history.

Flexible Mixture Model Approaches That Accommodate Footprint Size Variability for Robust Detection of Balancing Selection

Long-term balancing selection typically leaves narrow footprints of increased genetic diversity, and therefore most detection approaches only achieve optimal performances when sufficiently small genomic regions (i.e., windows) are examined. Such methods are sensitive to window sizes and suffer substantial losses in power when windows are large. Here, we employ mixture models to construct a set of five composite likelihood ratio test statistics, which we collectively term B statistics. These statistics are agnostic to window sizes and can operate on diverse forms of input data. Through simulations, we show that they exhibit comparable power to the best-performing current methods, and retain substantially high power regardless of window sizes. They also display considerable robustness to high mutation rates and uneven recombination landscapes, as well as an array of other common confounding scenarios. Moreover, we applied a specific version of the B statistics, termed B2, to a human population-genomic data set and recovered many top candidates from prior studies, including the then-uncharacterized STPG2 and CCDC169-SOHLH2, both of which are related to gamete functions. We further applied B2 on a bonobo population-genomic data set. In addition to the MHC-DQ genes, we uncovered several novel candidate genes, such as KLRD1, involved in viral defense, and SCN9A, associated with pain perception. Finally, we show that our methods can be extended to account for multiallelic balancing selection and integrated the set of statistics into open-source software named BalLeRMix for future applications by the scientific community.

The Genetic Mechanisms Driving Diversification of the KIR Gene Cluster in Primates

The activity and function of natural killer (NK) cells are modulated through the interactions of multiple receptor families, of which some recognize MHC class I molecules. The high level of MHC class I polymorphism requires their ligands either to interact with conserved epitopes, as is utilized by the NKG2A receptor family, or to co-evolve with the MHC class I allelic variation, which task is taken up by the killer cell immunoglobulin-like receptor (KIR) family. Multiple molecular mechanisms are responsible for the diversification of the KIR gene system, and include abundant chromosomal recombination, high mutation rates, alternative splicing, and variegated expression. The combination of these genetic mechanisms generates a compound array of diversity as is reflected by the contraction and expansion of KIR haplotypes, frequent birth of fusion genes, allelic polymorphism, structurally distinct isoforms, and variegated expression, which is in contrast to the mainly allelic nature of MHC class I polymorphism in humans. A comparison of the thoroughly studied human and macaque KIR gene repertoires demonstrates a similar evolutionarily conserved toolbox, through which selective forces drove and maintained the diversified nature of the KIR gene cluster. This hypothesis is further supported by the comparative genetics of KIR haplotypes and genes in other primate species. The complex nature of the KIR gene system has an impact upon the education, activity, and function of NK cells in coherence with an individual’s MHC class I repertoire and pathogenic encounters. Although selection operates on an individual, the continuous diversification of the KIR gene system in primates might protect populations against evolving pathogens.

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.

The HLA A03 Supertype and Several Pan Species Major Histocompatibility Complex Class I A Allotypes Share a Preference for Binding Positively Charged Residues in the F Pocket: Implications for Controlling Retroviral Infections

The major histocompatibility complex (MHC) class I region of humans, chimpanzees (Pan troglodytes), and bonobos (Pan paniscus) is highly similar, and orthologues of HLA-A, -B, and -C are present in both Pan species. Based on functional characteristics, the different HLA-A allotypes are classified into different supertypes. One of them, the HLA A03 supertype, is widely distributed among different human populations. All contemporary known chimpanzee and bonobo MHC class I A allotypes cluster genetically into one of the six HLA-A families, the HLA-A1/A3/A11/A30 family. We report here that the peptide-binding motif of the Patr-A*05:01 allotype, which is commonly present in a cohort of western African chimpanzees, has a strong preference for binding peptides with basic amino acids at the carboxyl terminus. This phenomenon is shared with the family members of the HLA A03 supertype. Based on the chemical similarities in the peptide-binding pocket, we inferred that the preference for binding peptides with basic amino acids at the carboxyl terminus is widely present among the human, chimpanzee, and bonobo MHC-A allotypes. Subsequent in silico peptide-binding predictions illustrated that these allotypes have the capacity to target conserved parts of the proteome of human immunodeficiency virus type 1 (HIV-1) and the simian immunodeficiency virus SIVcpz.IMPORTANCE Most experimentally infected chimpanzees seem to control an HIV-1 infection and are therefore considered to be relatively resistant to developing AIDS. Contemporary free-ranging chimpanzees may carry SIVcpz, and there is evidence for AIDS-like symptoms in these free-ranging animals, whereas SIV infections in bonobos appear to be absent. In humans, the natural control of an HIV-1 infection is strongly associated with the presence of particular HLA class I allotypes. The ancestor of the contemporary living chimpanzees and bonobos survived a selective sweep targeting the MHC class I repertoire. We have put forward a hypothesis that this may have been caused by an ancestral retroviral infection similar to SIVcpz. Characterization of the relevant MHC allotypes may contribute to understanding the shaping of their immune repertoire. The abundant presence of MHC-A allotypes that prefer peptides with basic amino acids at the C termini suggests that these molecules may contribute to the control of retroviral infections in humans, chimpanzees, and bonobos.

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.

Differences in MHC-B diversity and KIR epitopes in two populations of wild chimpanzees

The major histocompatibility complex (MHC) class I genes play a critical role within the immune system, both by the presentation of antigens from intracellular pathogens to immunocompetent cells and by the interaction with killer cell immunoglobulin-like receptors (KIR) on natural killer cells (NK cells). Genes of the MHC are highly diverse, and MHC variation can have effects on the immune functionality of individuals; hence, comparisons of MHC diversity among closely related phylogenetic taxa may give insight into the factors responsible for the shaping of its diversity. The four geographically separated chimpanzee subspecies differ in their overall genetic diversity, have different population histories, and are confronted with different pathogens in their natural habitat, all of which may affect MHC class I DNA sequence diversity. Here, we compare the MHC-B exon two DNA sequence diversity from 24 wild western and 46 wild eastern chimpanzees using necropsy and noninvasively collected fecal samples, respectively. We found a higher MHC-B exon two nucleotide diversity, in our western than eastern chimpanzees. The inclusion of previously published MHC-B exon two data from other western and eastern chimpanzees supported this finding. In addition, our results confirm and extend the finding of a very low C1 epitope frequency at eastern chimpanzee MHC-B molecules, which likely affects the ability of these molecules to interact with NK cells. While the understanding of the differing pathogen environments encountered by disparate populations of a species is a challenging endeavor, these findings highlight the potential for these pathogens to selectively shape immune system variation.

Nomenclature report 2019: major histocompatibility complex genes and alleles of Great and Small Ape and Old and New World monkey species

The major histocompatibility complex (MHC) is central to the innate and adaptive immune responses of jawed vertebrates. Characteristic of the MHC are high gene density, gene copy number variation, and allelic polymorphism. Because apes and monkeys are the closest living relatives of humans, the MHCs of these non-human primates (NHP) are studied in depth in the context of evolution, biomedicine, and conservation biology. The Immuno Polymorphism Database (IPD)-MHC NHP Database (IPD-MHC NHP), which curates MHC data of great and small apes, as well as Old and New World monkeys, has been upgraded. The curators of the database are responsible for providing official designations for newly discovered alleles. This nomenclature report updates the 2012 report, and summarizes important nomenclature issues and relevant novel features of the IPD-MHC NHP Database.

Two to Tango: Co-evolution of Hominid Natural Killer Cell Receptors and MHC

Natural killer (NK) cells have diverse roles in hominid immunity and reproduction. Modulating these functions are the interactions between major histocompatibility complex (MHC) class I molecules that are ligands for two NK cell surface receptor types. Diverse killer cell immunoglobulin-like receptors (KIR) bind specific motifs encoded within the polymorphic MHC class I cell surface glycoproteins, while, in more conserved interactions, CD94:NKG2A receptors recognize MHC-E with bound peptides derived from MHC class I leader sequences. The hominid lineage presents a choreographed co-evolution of KIR with their MHC class I ligands. MHC-A, -B, and -C are present in all great apes with species-specific haplotypic variation in gene content. The Bw4 epitope recognized by lineage II KIR is restricted to MHC-B but also present on some gorilla and human MHC-A. Common to great apes, but rare in humans, are MHC-B possessing a C1 epitope recognized by lineage III KIR. MHC-C arose from duplication of MHC-B and is fixed in all great apes except orangutan, where it exists on approximately 50% of haplotypes and all allotypes are C1-bearing. Recent study showed that gorillas possess yet another intermediate MHC organization compared to humans. Like orangutans, but unlike the Pan-Homo species, duplication of MHC-B occurred. However, MHC-C is fixed, and the MHC-C C2 epitope (absent in orangutans) emerges. The evolution of MHC-C drove expansion of its cognate lineage III KIR. Recently, position −21 of the MHC-B leader sequence has been shown to be critical in determining NK cell educational outcome. In humans, methionine (−21M) results in CD94:NKG2A-focused education whereas threonine (−21T) produces KIR-focused education. This is another dynamic position among hominids. Orangutans have exclusively −21M, consistent with their intermediate stage in lineage III KIR-focused evolution. Gorillas have both −21M and −21T, like humans, but they are unequally encoded by their duplicated B genes. Chimpanzees have near-fixed −21T, indicative of KIR-focused NK education. Harmonious with this observation, chimpanzee KIR exhibit strong binding and, compared to humans, smaller differences between binding levels of activating and inhibitory KIR. Consistent between these MHC-NK cell receptor systems over the course of hominid evolution is the evolution of polymorphism favoring the more novel and dynamic KIR system.

Reduced bonobo MHC class I diversity predicts a reduced viral peptide binding ability compared to chimpanzees

BACKGROUND

The highly polymorphic genes of the major histocompatibility complex (MHC) class I are involved in defense against viruses and other intracellular pathogens. Although several studies found reduced MHC class I diversity in bonobos in comparison to the closely related chimpanzee, it is unclear if this lower diversity also influences the functional ability of MHC class I molecules in bonobos. Here, we use a bioinformatic approach to analyze the viral peptide binding ability of all published bonobo MHC class I molecules (n = 58) in comparison to all published chimpanzee MHC class I molecules (n = 161) for the class I loci A, B, C and A-like.

RESULTS

We examined the peptide binding ability of all 219 different MHC class I molecules to 5,788,712 peptides derived from 1432 different primate viruses and analyzed the percentage of bound peptides and the overlap of the peptide binding repertoires of the two species. We conducted multiple levels of analysis on the "species"-, "population"- and "individual"-level to account for the characterization of MHC variation in a larger number of chimpanzees and their broader geographic distribution. We found a lower percentage of bound peptides in bonobos at the B locus in the "population"-level comparison and at the B and C loci in the "individual"-level comparison. Furthermore, we found evidence of a limited peptide binding repertoire in bonobos by tree-based visualization of functional clustering of MHC molecules, as well as an analysis of peptides bound by both species.

CONCLUSION

Our results suggest a reduced MHC class I viral peptide binding ability at the B and C loci in bonobos compared to chimpanzees. The effects of this finding on the immune defense against viruses in wild living bonobos are unclear. However, special caution is needed to prevent introduction and spread of new viruses to bonobos, as their defensive ability to cope with new viruses could be limited compared to chimpanzees.

Killer-cell immunoglobulin-like receptors on the cusp of modern immunogenetics

Killer-cell immunoglobulin-like receptors (KIR) and their human leukocyte antigen (HLA) ligands play a central role in immunity and human health. These molecules are encoded by gene families with copy number variation, extreme levels of sequence diversity and complex expression patterns. The rapid evolution of KIR and HLA genes and their associations with infectious diseases, pregnancy disorders, immunopathologies and outcome of cell transplantation have generated considerable interest from immunologists, geneticists and clinicians. Until recently, however, analyses have been stuck at low-level resolution, focusing primarily on presence or absence of KIR genes. This is changing with the advent of modern high throughput sequencing, cell phenotyping and bioinformatics. These developments allow high-resolution analysis and much deeper understanding of KIR evolution and KIR function. The impending deluge of high dimensional data brings inevitably new challenges in analysis, interpretation and communication of results, but the benefits are already tangible. The diversity of KIR across worldwide human populations is being catalogued at the allele level. Structures of KIR molecules and their interactions with HLA-peptide complexes are being determined. How KIR modulate natural killer cell education is being defined. Ligands for activating KIR, elusive for many years, are being discovered. KIR gene complexes and their related receptor gene families are being characterized in animal models and livestock breeds. These advances are helping to generate a more complete picture of the impact of KIR variation in health and disease and offer new opportunities for immunotherapy, as highlighted in a recent meeting (The Tenth KIR Workshop, April 2017 Cambridge, UK).

AIDS in chimpanzees: the role of MHC genes

The ancestral progenitor of common chimpanzees and bonobos experienced a selective sweep that ravaged its major histocompatibility complex (MHC) class I repertoire. The causative agent was probably an ancestral retrovirus, highly related to the contemporary HIV-1 strain, which initiated the acquired immunodeficiency syndrome pandemic in the human population. As a direct result, MHC class I allotypes with the capability of targeting conserved retroviral elements were enriched in the ancestral progenitor. Even today, the impact can be traced back by studying the functional capacities of the contemporary MHC class I allotypes of common chimpanzees. Viruses, however, have developed several strategies to manipulate the cell-surface expression of MHC class I genes. Monitoring the presence and absence of the MHC class I allotypes on the cell surface is conducted, for instance, by the hosts' gene products of the killer cell immunoglobulin-like receptor (KIR) complex. Hence, one may wonder whether-in the future-any clues with regard to the signature of the MHC class I selective sweep might be unearthed for the KIR genes as well.

Limited MHC class I intron 2 repertoire variation in bonobos

Common chimpanzees (Pan troglodytes) experienced a selective sweep, probably caused by a SIV-like virus, which targeted their MHC class I repertoire. Based on MHC class I intron 2 data analyses, this selective sweep took place about 2-3 million years ago. As a consequence, common chimpanzees have a skewed MHC class I repertoire that is enriched for allotypes that are able to recognise conserved regions of the SIV proteome. The bonobo (Pan paniscus) shared an ancestor with common chimpanzees approximately 1.5 to 2 million years ago. To investigate whether the signature of this selective sweep is also detectable in bonobos, the MHC class I gene repertoire of two bonobo panels comprising in total 29 animals was investigated by Sanger sequencing. We identified 14 Papa-A, 20 Papa-B and 11 Papa-C alleles, of which eight, five and eight alleles, respectively, have not been reported previously. Within this pool of MHC class I variation, we recovered only 2 Papa-A, 3 Papa-B and 6 Papa-C intron 2 sequences. As compared to humans, bonobos appear to have an even more diminished MHC class I intron 2 lineage repertoire than common chimpanzees. This supports the notion that the selective sweep may have predated the speciation of common chimpanzees and bonobos. The further reduction of the MHC class I intron 2 lineage repertoire observed in bonobos as compared to the common chimpanzee may be explained by a founding effect or other subsequent selective processes.

MHC class I diversity in chimpanzees and bonobos

Major histocompatibility complex (MHC) class I genes are critically involved in the defense against intracellular pathogens. MHC diversity comparisons among samples of closely related taxa may reveal traces of past or ongoing selective processes. The bonobo and chimpanzee are the closest living evolutionary relatives of humans and last shared a common ancestor some 1 mya. However, little is known concerning MHC class I diversity in bonobos or in central chimpanzees, the most numerous and genetically diverse chimpanzee subspecies. Here, we used a long-read sequencing technology (PacBio) to sequence the classical MHC class I genes A, B, C, and A-like in 20 and 30 wild-born bonobos and chimpanzees, respectively, with a main focus on central chimpanzees to assess and compare diversity in those two species. We describe in total 21 and 42 novel coding region sequences for the two species, respectively. In addition, we found evidence for a reduced MHC class I diversity in bonobos as compared to central chimpanzees as well as to western chimpanzees and humans. The reduced bonobo MHC class I diversity may be the result of a selective process in their evolutionary past since their split from chimpanzees.