Gorilla MHC class I gene and sequence variation in a comparative context

Using PacBio SMRT data for identification of class I MHC alleles in a wildlife species, Zalophus californianus (California sea lion)

High allelic polymorphism and association with disease susceptibility has made the genes encoding major histocompatibility complex (MHC) antigen presentation molecules in humans, domesticated animals, and wildlife species of wide interest to ecologists, evolutionary biologists, and health specialists. The often multifaceted polygenism and extreme polymorphism of this immunogenetic system have made it especially difficult to characterize in non-model species. Here we compare and contrast the workflows of traditional Sanger sequencing of plasmid-cloned amplicons to Pacific Biosciences SMRT circular consensus sequencing (CCS) in their ability to capture alleles of MHC class I in a wildlife species where characterization of these genes was absent. We assessed two California sea lions (Zalophus californianus), a species suffering from a high prevalence of an aggressive cancer associated with a sexually transmitted gamma herpesvirus. In this pilot study, SMRT CCS proved superior in identifying more alleles from each animal than the more laborious plasmid cloning/Sanger workflow (12:7, 10:7), and no alleles were identified with the cloning/Sanger approach that were not identified by SMRT CCS. We discuss the advantages and disadvantages of each approach including cost, allele rarefaction, and sequence fidelity.

Enteroviruses from Humans and Great Apes in the Republic of Congo: Recombination within Enterovirus C Serotypes

Enteroviruses (EVs) are viruses of the family Picornaviridae that cause mild to severe infections in humans and in several animal species, including non-human primates (NHPs). We conducted a survey and characterization of enteroviruses circulating between humans and great apes in the Congo. Fecal samples (N = 24) of gorillas and chimpanzees living close to or distant from humans in three Congolese parks were collected, as well as from healthy humans (N = 38) living around and within these parks. Enteroviruses were detected in 29.4% of gorilla and 13.15% of human feces, including wild and human-habituated gorillas, local humans and eco-guards. Two identical strains were isolated from two humans coming from two remote regions. Their genomes were similar and all genes showed their close similarity to coxsackieviruses, except for the 3C, 3D and 5′-UTR regions, where they were most similar to poliovirus 1 and 2, suggesting recombination. Recombination events were found between these strains, poliovirus 1 and 2 and EV-C99. It is possible that the same EV-C species circulated in both humans and apes in different regions in the Congo, which must be confirmed in other investigations. In addition, other studies are needed to further investigate the circulation and genetic diversity of enteroviruses in the great ape population, to draw a definitive conclusion on the different species and types of enteroviruses circulating in the Republic of Congo.

Similar patterns of genetic diversity and linkage disequilibrium in Western chimpanzees (Pan troglodytes verus) and humans indicate highly conserved mechanisms of MHC molecular evolution

BACKGROUND

Many species are threatened with extinction as their population sizes decrease with changing environments or face novel pathogenic threats. A reduction of genetic diversity at major histocompatibility complex (MHC) genes may have dramatic effects on populations' survival, as these genes play a key role in adaptive immunity. This might be the case for chimpanzees, the MHC genes of which reveal signatures of an ancient selective sweep likely due to a viral epidemic that reduced their population size a few million years ago. To better assess how this past event affected MHC variation in chimpanzees compared to humans, we analysed several indexes of genetic diversity and linkage disequilibrium across seven MHC genes on four cohorts of chimpanzees and we compared them to those estimated at orthologous HLA genes in a large set of human populations.

RESULTS

Interestingly, the analyses uncovered similar patterns of both molecular diversity and linkage disequilibrium across the seven MHC genes in chimpanzees and humans. Indeed, in both species the greatest allelic richness and heterozygosity were found at loci A, B, C and DRB1, the greatest nucleotide diversity at loci DRB1, DQA1 and DQB1, and both significant global linkage disequilibrium and the greatest proportions of haplotypes in linkage disequilibrium were observed at pairs DQA1 ~ DQB1, DQA1 ~ DRB1, DQB1 ~ DRB1 and B ~ C. Our results also showed that, despite some differences among loci, the levels of genetic diversity and linkage disequilibrium observed in contemporary chimpanzees were globally similar to those estimated in small isolated human populations, in contrast to significant differences compared to large populations.

CONCLUSIONS

We conclude, first, that highly conserved mechanisms shaped the diversity of orthologous MHC genes in chimpanzees and humans. Furthermore, our findings support the hypothesis that an ancient demographic decline affecting the chimpanzee populations - like that ascribed to a viral epidemic - exerted a substantial effect on the molecular diversity of their MHC genes, albeit not more pronounced than that experienced by HLA genes in human populations that underwent rapid genetic drift during humans' peopling history. We thus propose a model where chimpanzees' MHC genes regenerated molecular variation through recombination/gene conversion and/or balancing selection after the selective sweep.

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.

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.

HLAIb worldwide genetic diversity: New HLA-H alleles and haplotype structure description

The classical HLA class I genes (HLA Ia) were extensively studied because of their implication in clinical fields and anthropology. Less is known about worldwide genetic diversity and linkage disequilibrium for non-classical HLA class I genes (HLA Ib) and HLA pseudogenes. Notably, HLA-H, which is deleted in a fraction of the population, remains scarcely explored. The aims of this study were 1/ to get further insight into HLA-H genetic diversity and into how this variability potentially affects its expression and 2/ to define HLA Ib worldwide allelic diversity and linkage. Exome sequence data from the 1000 Genomes Project were used to define second field HLA-A, -E, -F, -G and -H typing using PolyPheMe software. Allelic and two-loci haplotype frequencies were estimated using Gene[Rate] software both at worldwide and continental levels. Eleven novel HLA-H alleles identified in exome data were validated by NGS performed on 25 genomic DNA samples from the same cohort. Phylogenetic analysis and frequency distribution of HLA-H alleles revealed three clades, each predominantly represented in Admixed American, European and East Asian populations, African populations and South Asian populations. Among these eleven novel alleles, two potentially encode complete transmembrane HLA proteins. We confirm the high LD between HLA-H and -A, and between HLA-H and -G, and show the three genes have distinct worldwide allelic distribution. Conversely, HLA-E and HLA-F both showed little LD, displayed restricted allelic diversity and practically no difference in their distribution across the planet. Our work thus reveals an unexpectedly high HLA-H genetic diversity, with alleles highly represented in Asia possibly encoding a functional HLA protein. Functional implication of these results remains to be explored, both in physiological and pathological contexts.

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.

Discovery of gorilla MHC-C expressing C1 ligand for KIR

In comparison to humans and chimpanzees, gorillas show low diversity at MHC class I genes (Gogo), as reflected by an overall reduced level of allelic variation as well as the absence of a functionally important sequence motif that interacts with killer cell immunoglobulin-like receptors (KIR). Here, we use recently generated large-scale genomic sequence data for a reassessment of allelic diversity at Gogo-C, the gorilla orthologue of HLA-C. Through the combination of long-range amplifications and long-read sequencing technology, we obtained, among the 35 gorillas reanalyzed, three novel full-length genomic sequences including a coding region sequence that has not been previously described. The newly identified Gogo-C*03:01 allele has a divergent recombinant structure that sets it apart from other Gogo-C alleles. Domain-by-domain phylogenetic analysis shows that Gogo-C*03:01 has segments in common with Gogo-B*07, the additional B-like gene that is present on some gorilla MHC haplotypes. Identified in ~ 50% of the gorillas analyzed, the Gogo-C*03:01 allele exclusively encodes the C1 epitope among Gogo-C allotypes, indicating its important function in controlling natural killer cell (NK cell) responses via KIR. We further explored the hypothesis whether gorillas experienced a selective sweep which may have resulted in a general reduction of the gorilla MHC class I repertoire. Our results provide little support for a selective sweep but rather suggest that the overall low Gogo class I diversity can be best explained by drastic demographic changes gorillas experienced in the ancient and recent past.

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.

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.