Improved full-length killer cell immunoglobulin-like receptor transcript discovery in Mauritian cynomolgus macaques

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.

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.

Identifying a Minor Histocompatibility Antigen in Mauritian Cynomolgus Macaques Encoded by APOBEC3C

Allogeneic hematopoietic stem cell transplants can lead to dramatic reductions in human immunodeficiency virus (HIV) reservoirs. This effect is partially mediated by donor T cells recognizing lymphocyte-expressed minor histocompatibility antigens (mHAgs). The potential to mark malignant and latently infected cells for destruction makes mHAgs attractive targets for cellular immunotherapies. However, testing such HIV reservoir reduction strategies will likely require preclinical studies in non-human primates (NHPs). In this study, we used a combination of alloimmunization, whole exome sequencing, and bioinformatics to identify an mHAg in Mauritian cynomolgus macaques (MCMs). We mapped the minimal optimal epitope to a 10-mer peptide (SW10) in apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3C (APOBEC3C) and determined the major histocompatibility complex class I restriction element as Mafa-A1∗063, which is expressed in almost 90% of MCMs. APOBEC3C SW10-specific CD8+ T cells recognized immortalized B cells but not fibroblasts from an mHAg-positive MCM. These results provide a framework for identifying mHAgs in a non-transplant setting and suggest that APOBEC3C SW10 could be used as a model antigen to test mHAg-targeted therapies in NHPs.

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.

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.

Nomenclature report for killer-cell immunoglobulin-like receptors (KIR) in macaque species: new genes/alleles, renaming recombinant entities and IPD-NHKIR updates

The Killer-cell Immunoglobulin-like Receptors (KIR) are encoded by a diverse group of genes, which are characterized by allelic polymorphism, gene duplications, and recombinations, which may generate recombinant entities. The number of reported macaque KIR sequences is steadily increasing, and these data illustrate a gene system that may match or exceed the complexity of the human KIR cluster. This report lists the names of quality controlled and annotated KIR genes/alleles with all the relevant references for two different macaque species: rhesus and cynomolgus macaques. Numerous recombinant KIR genes in these species necessitate a revision of some of the earlier-published nomenclature guidelines. In addition, this report summarizes the latest information on the Immuno Polymorphism Database (IPD)-NHKIR Database, which contains annotated KIR sequences from four non-human primate species.

Evaluation of an antibody to α4β7 in the control of SIVmac239-nef-stop infection

Treatment of SIV-infected rhesus macaques with short-term antiretroviral therapy (ART) and partially overlapping infusions of antibody to integrin α4β7 was reported to induce durable posttreatment viral suppression. In an attempt to replicate those observations, we treated macaques infected with the same virus and with the same ART and monoclonal antibody (mAb) regimens (anti-α4β7 versus control mAb). Sequencing demonstrated that the virus used was actually SIVmac239-nef-stop, not wild-type SIVmac239. A positive correlation was found at 2 weeks after infection between the frequency of repair of attenuated Nef-STOP virus to pathogenic Nef-OPEN and plasma SIV RNA levels. Levels of plasma viremia before the first antibody infusion and preinfection levels of α4β7hi CD4+ T cells, but not treatment with antibody to α4β7 , correlated with levels of viral replication upon discontinuation of all treatments. Follow-up plasma viremia, peripheral blood CD4+ T cell counts, and lymph node and rectal tissue viral load were not significantly different between anti-α4β7 and control mAb groups.

MHC genotyping from rhesus macaque exome sequences

Indian rhesus macaque major histocompatibility complex (MHC) variation can influence the outcomes of transplantation and infectious disease studies. Frequently, rhesus macaques are MHC genotyped to identify variants that could account for unexpected results. Since the MHC is only one region in the genome where variation could impact experimental outcomes, strategies for simultaneously profiling variation in the macaque MHC and the remainder of the protein coding genome would be useful. Here we determine MHC class I and class II genotypes using target-capture probes enriched for MHC sequences, a method we term macaque exome sequence (MES) genotyping. For a cohort of 27 Indian rhesus macaques, we describe two methods for obtaining MHC genotypes from MES data and demonstrate that the MHC class I and class II genotyping results obtained with these methods are 98.1% and 98.7% concordant, respectively, with expected MHC genotypes. In contrast, conventional MHC genotyping results obtained by deep sequencing of short multiplex PCR amplicons were only 92.6% concordant with expectations for this cohort.

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.

Characterization of Mauritian Cynomolgus Macaque FcγR Alleles Using Long-Read Sequencing

The FcγRs are immune cell surface proteins that bind IgG and facilitate cytokine production, phagocytosis, and Ab-dependent, cell-mediated cytotoxicity. FcγRs play a critical role in immunity; variation in these genes is implicated in autoimmunity and other diseases. Cynomolgus macaques are an excellent animal model for many human diseases, and Mauritian cynomolgus macaques (MCMs) are particularly useful because of their restricted genetic diversity. Previous studies of MCM immune gene diversity have focused on the MHC and killer cell Ig-like receptor. In this study, we characterize FcγR diversity in 48 MCMs using PacBio long-read sequencing to identify novel alleles of each of the four expressed MCM FcγR genes. We also developed a high-throughput FcγR genotyping assay, which we used to determine allele frequencies and identify FcγR haplotypes in more than 500 additional MCMs. We found three alleles for FcγR1A, seven each for FcγR2A and FcγR2B, and four for FcγR3A; these segregate into eight haplotypes. We also assessed whether different FcγR alleles confer different Ab-binding affinities by surface plasmon resonance and found minimal difference in binding affinities across alleles for a panel of wild type and Fc-engineered human IgG. This work suggests that although MCMs may not fully represent the diversity of FcγR responses in humans, they may offer highly reproducible results for mAb therapy and toxicity studies.

Extensive Alternative Splicing of KIR Transcripts

The killer-cell Ig-like receptors (KIR) form a multigene entity involved in modulating immune responses through interactions with MHC class I molecules. The complexity of the KIR cluster is reflected by, for instance, abundant levels of allelic polymorphism, gene copy number variation, and stochastic expression profiles. The current transcriptome study involving human and macaque families demonstrates that KIR family members are also subjected to differential levels of alternative splicing, and this seems to be gene dependent. Alternative splicing may result in the partial or complete skipping of exons, or the partial inclusion of introns, as documented at the transcription level. This post-transcriptional process can generate multiple isoforms from a single KIR gene, which diversifies the characteristics of the encoded proteins. For example, alternative splicing could modify ligand interactions, cellular localization, signaling properties, and the number of extracellular domains of the receptor. In humans, we observed abundant splicing for KIR2DL4, and to a lesser extent in the lineage III KIR genes. All experimentally documented splice events are substantiated by in silico splicing strength predictions. To a similar extent, alternative splicing is observed in rhesus macaques, a species that shares a close evolutionary relationship with humans. Splicing profiles of Mamu-KIR1D and Mamu-KIR2DL04 displayed a great diversity, whereas Mamu-KIR3DL20 (lineage V) is consistently spliced to generate a homolog of human KIR2DL5 (lineage I). The latter case represents an example of convergent evolution. Although just a single KIR splice event is shared between humans and macaques, the splicing mechanisms are similar, and the predicted consequences are comparable. In conclusion, alternative splicing adds an additional layer of complexity to the KIR gene system in primates, and results in a wide structural and functional variety of KIR receptors and its isoforms, which may play a role in health and disease.

Novel full-length major histocompatibility complex class I allele discovery and haplotype definition in pig-tailed macaques

Pig-tailed macaques (Macaca nemestrina, Mane) are important models for human immunodeficiency virus (HIV) studies. Their infectability with minimally modified HIV makes them a uniquely valuable animal model to mimic human infection with HIV and progression to acquired immunodeficiency syndrome (AIDS). However, variation in the pig-tailed macaque major histocompatibility complex (MHC) and the impact of individual transcripts on the pathogenesis of HIV and other infectious diseases is understudied compared to that of rhesus and cynomolgus macaques. In this study, we used Pacific Biosciences single-molecule real-time circular consensus sequencing to describe full-length MHC class I (MHC-I) transcripts for 194 pig-tailed macaques from three breeding centers. We then used the full-length sequences to infer Mane-A and Mane-B haplotypes containing groups of MHC-I transcripts that co-segregate due to physical linkage. In total, we characterized full-length open reading frames (ORFs) for 313 Mane-A, Mane-B, and Mane-I sequences that defined 86 Mane-A and 106 Mane-B MHC-I haplotypes. Pacific Biosciences technology allows us to resolve these Mane-A and Mane-B haplotypes to the level of synonymous allelic variants. The newly defined haplotypes and transcript sequences containing full-length ORFs provide an important resource for infectious disease researchers as certain MHC haplotypes have been shown to provide exceptional control of simian immunodeficiency virus (SIV) replication and prevention of AIDS-like disease in nonhuman primates. The increased allelic resolution provided by Pacific Biosciences sequencing also benefits transplant research by allowing researchers to more specifically match haplotypes between donors and recipients to the level of nonsynonymous allelic variation, thus reducing the risk of graft-versus-host disease.

Restricted MHC class I A locus diversity in olive and hybrid olive/yellow baboons from the Southwest National Primate Research Center

Baboons are valuable models for complex human diseases due to their genetic and physiologic similarities to humans. Deep sequencing methods to characterize full-length major histocompatibility complex (MHC) class I (MHC-I) alleles in different nonhuman primate populations were used to identify novel MHC-I alleles in baboons. We combined data from Illumina MiSeq sequencing and Roche/454 sequencing to characterize novel full-length MHC-I transcripts in a cohort of olive and hybrid olive/yellow baboons from the Southwest National Primate Research Center (SNPRC). We characterized 57 novel full-length alleles from 24 baboons and found limited genetic diversity at the MHC-I A locus, with significant sharing of two MHC-I A lineages between 22 out of the 24 animals characterized. These shared alleles provide the basis for development of tools such as MHC:peptide tetramers for studying cellular immune responses in this important animal model.

Current and Future Methods for mRNA Analysis: A Drive Toward Single Molecule Sequencing

The transcriptome encompasses a range of species including messenger RNA, and other noncoding RNA such as rRNA, tRNA, and short and long noncoding RNAs. Due to the huge role played by mRNA in development and disease, several methods have been developed to sequence and characterize mRNA, with RNA sequencing (RNA-Seq) emerging as the current method of choice particularly for large high-throughput studies. Short-read RNA-Seq which involves sequencing of short cDNA fragments and computationally assembling them to reconstruct the transcriptome, or aligning them to a reference is the most widely used approach. However, due to inherent limitations of this approach in de novo transcriptome assembly and isoform quantification, long-read RNA-Seq approaches, which also happen to be single molecule sequencing approaches, are increasingly becoming the standard for de novo transcriptome assembly and isoform quantification. In this chapter, we review the technical aspects of the current methods of RNA-Seq, both short and long-read approaches, and data analysis methods available. We discuss recent advances in single-cell RNA-Seq and direct RNA-Seq approaches, which perhaps will dominate the future of RNA-Seq.

KIR3DL01 upregulation on gut natural killer cells in response to SIV infection of KIR- and MHC class I-defined rhesus macaques

Natural killer cells provide an important early defense against viral pathogens and are regulated in part by interactions between highly polymorphic killer-cell immunoglobulin-like receptors (KIRs) on NK cells and their MHC class I ligands on target cells. We previously identified MHC class I ligands for two rhesus macaque KIRs: KIR3DL01 recognizes Mamu-Bw4 molecules and KIR3DL05 recognizes Mamu-A1*002. To determine how these interactions influence NK cell responses, we infected KIR3DL01+ and KIR3DL05+ macaques with and without defined ligands for these receptors with SIVmac239, and monitored NK cell responses in peripheral blood and lymphoid tissues. NK cell responses in blood were broadly stimulated, as indicated by rapid increases in the CD16+ population during acute infection and sustained increases in the CD16+ and CD16-CD56- populations during chronic infection. Markers of proliferation (Ki-67), activation (CD69 & HLA-DR) and antiviral activity (CD107a & TNFα) were also widely expressed, but began to diverge during chronic infection, as reflected by sustained CD107a and TNFα upregulation by KIR3DL01+, but not by KIR3DL05+ NK cells. Significant increases in the frequency of KIR3DL01+ (but not KIR3DL05+) NK cells were also observed in tissues, particularly in the gut-associated lymphoid tissues, where this receptor was preferentially upregulated on CD56+ and CD16-CD56- subsets. These results reveal broad NK cell activation and dynamic changes in the phenotypic properties of NK cells in response to SIV infection, including the enrichment of KIR3DL01+ NK cells in tissues that support high levels of virus replication.