Whole genome sequencing of SIV-infected macaques identifies candidate loci that may contribute to host control of virus replication

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.

CD8+ cells and small viral reservoirs facilitate post-ART control of SIV replication in M3+ Mauritian cynomolgus macaques initiated on ART two weeks post-infection

Sustainable HIV remission after antiretroviral therapy (ART) withdrawal, or post-treatment control (PTC), remains a top priority for HIV treatment. We observed surprising PTC in an MHC-haplomatched cohort of MHC-M3+ SIVmac239+ Mauritian cynomolgus macaques (MCMs) initiated on ART at two weeks post-infection (wpi). None of the MCMs possessed MHC haplotypes previously associated with SIV control. For six months after ART withdrawal, we observed undetectable or transient viremia in seven of the eight MCMs, despite detecting replication competent SIV using quantitative viral outgrowth assays. In vivo depletion of CD8α+ cells induced rebound in all animals, indicating the observed PTC was mediated, at least in part, by CD8α+ cells. With intact proviral DNA assays, we found that MCMs had significantly smaller viral reservoirs two wpi than a cohort of identically infected rhesus macaques, a population that rarely develops PTC. We found a similarly small viral reservoir among six additional SIV+ MCMs in which ART was initiated at eight wpi, some of whom exhibited viral rebound. These results suggest that an unusually small viral reservoir is a hallmark among SIV+ MCMs. By evaluating immunological differences between MCMs that did and did not rebound, we identified that PTC was associated with a reduced frequency of CD4+ and CD8+ lymphocyte subsets expressing exhaustion markers. Together, these results suggest a combination of small reservoirs and immune-mediated virus suppression contribute to PTC in MCMs. Further, defining the immunologic mechanisms that engender PTC in this model may identify therapeutic targets for inducing durable HIV remission in humans.

CD8+ cells and small viral reservoirs facilitate post-ART control of SIV in Mauritian cynomolgus macaques

Sustainable HIV remission after antiretroviral therapy (ART) withdrawal, or post-treatment control (PTC), remains a top priority for HIV treatment. We observed surprising PTC in an MHC-haplomatched cohort of MHC-M3+ SIVmac239+ Mauritian cynomolgus macaques (MCMs) initiated on ART at two weeks post-infection (wpi). For six months after ART withdrawal, we observed undetectable or transient viremia in seven of eight MCMs. In vivo depletion of CD8α+ cells induced rebound in all animals, indicating the PTC was mediated, at least in part, by CD8α+ cells. We found that MCMs had smaller acute viral reservoirs than a cohort of identically infected rhesus macaques, a population that rarely develops PTC. The mechanisms by which unusually small viral reservoirs and CD8α+ cell-mediated virus suppression enable PTC can be investigated using this MHC-haplomatched MCM model. Further, defining the immunologic mechanisms that engender PTC in this model may identify therapeutic targets for inducing durable HIV remission in humans.

Complete sequencing of a cynomolgus macaque major histocompatibility complex haplotype

Macaques provide the most widely used nonhuman primate models for studying the immunology and pathogenesis of human diseases. Although the macaque major histocompatibility complex (MHC) region shares most features with the human leukocyte antigen (HLA) region, macaques have an expanded repertoire of MHC class I genes. Although a chimera of two rhesus macaque MHC haplotypes was first published in 2004, the structural diversity of MHC genomic organization in macaques remains poorly understood owing to a lack of adequate genomic reference sequences. We used ultralong Oxford Nanopore and high-accuracy Pacific Biosciences (PacBio) HiFi sequences to fully assemble the ∼5.2-Mb M3 haplotype of an MHC-homozygous, Mauritian-origin cynomolgus macaque (Macaca fascicularis). The MHC homozygosity allowed us to assemble a single MHC haplotype unambiguously and avoid chimeric assemblies that hampered previous efforts to characterize this exceptionally complex genomic region in macaques. The high quality of this new assembly is exemplified by the identification of an extended cluster of six Mafa-AG genes that contains a recent duplication with a highly similar ∼48.5-kb block of sequence. The MHC class II region of this M3 haplotype is similar to the previously sequenced rhesus macaque haplotype and HLA class II haplotypes. The MHC class I region, in contrast, contains 13 MHC-B genes, four MHC-A genes, and three MHC-E genes (vs. 19 MHC-B, two MHC-A, and one MHC-E in the previously sequenced haplotype). These results provide an unambiguously assembled single contiguous cynomolgus macaque MHC haplotype with fully curated gene annotations that will inform infectious disease and transplantation research.

Testing Sex-Biased Admixture Origin of Macaque Species Using Autosomal and X-Chromosomal Genomic Sequences

The role of sex-specific demography in hybridization and admixture of genetically diverged species and populations is essential to understand the origins of the genomic diversity of sexually reproducing organisms. In order to infer how sex-linked loci have been differentiated undergoing frequent hybridization and admixture, we examined 17 whole-genome sequences of seven species representing the genus Macaca, which shows frequent inter-specific hybridization and predominantly female philopatry. We found that hybridization and admixture were prevalent within these species. For three cases of suggested hybrid origin of species/subspecies, Macaca arctoides, Macaca fascicularis ssp. aurea, and Chinese Macaca mulatta, we examined the level of admixture of X chromosomes, which is less affected by male-biased migration than that of autosomes. In one case, we found that Macaca cyclopis and Macaca fuscata was genetically closer to Chinese M. mulatta than to the Indian M. mulatta, and the admixture level of Chinese M. mulatta and M. fuscata/cyclopis was more pronounced on the X chromosome than on autosomes. Since the mitochondrial genomes of Chinese M. mulatta, M. cyclopis, and M. fuscata were found to cluster together, and the mitochondrial genome of Indian M. mulatta is more distantly related, the observed pattern of genetic differentiation on X-chromosomal loci is consistent with the nuclear swamping hypothesis, in which strong, continuous male-biased introgression from the ancestral Chinese M. mulatta population to a population related to M. fuscata and M. cyclopis generated incongruencies between the genealogies of the mitochondrial and nuclear genomes.

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.

Unusual sequence characteristics of human chromosome 19 are conserved across 11 nonhuman primates

BACKGROUND

Human chromosome 19 has many unique characteristics including gene density more than double the genome-wide average and 20 large tandemly clustered gene families. It also has the highest GC content of any chromosome, especially outside gene clusters. The high GC content and concomitant high content of hypermutable CpG sites raises the possibility chromosome 19 exhibits higher levels of nucleotide diversity both within and between species, and may possess greater variation in DNA methylation that regulates gene expression.

RESULTS

We examined GC and CpG content of chromosome 19 orthologs across representatives of the primate order. In all 12 primate species with suitable genome assemblies, chromosome 19 orthologs have the highest GC content of any chromosome. CpG dinucleotides and CpG islands are also more prevalent in chromosome 19 orthologs than other chromosomes. GC and CpG content are generally higher outside the gene clusters. Intra-species variation based on SNPs in human common dbSNP, rhesus, crab eating macaque, baboon and marmoset datasets is most prevalent on chromosome 19 and its orthologs. Inter-species comparisons based on phyloP conservation show accelerated nucleotide evolution for chromosome 19 promoter flanking and enhancer regions. These same regulatory regions show the highest CpG density of any chromosome suggesting they possess considerable methylome regulatory potential.

CONCLUSIONS

The pattern of high GC and CpG content in chromosome 19 orthologs, particularly outside gene clusters, is present from human to mouse lemur representing 74 million years of primate evolution. Much CpG variation exists both within and between primate species with a portion of this variation occurring in regulatory regions.

Immunization of BLT Humanized Mice Redirects T Cell Responses to Gag and Reduces Acute HIV-1 Viremia

BLT (bone marrow-liver-thymus) humanized mice, which reconstitute a functional human immune system, develop prototypic human virus-specific CD8⁺ T cell responses following infection with human immunodeficiency virus type 1 (HIV-1). We explored the utility of the BLT model for HIV-1 vaccine development by immunizing BLT mice against the conserved viral Gag protein, utilizing a rapid prime-boost protocol of poly(lactic-co-glycolic) acid microparticles and a replication-defective herpes simplex virus (HSV) recombinant vector. After HIV-1 challenge, the mice developed broad, proteome-wide gamma interferon-positive (IFN-γ⁺) T cell responses against HIV-1 that reached magnitudes equivalent to what is observed in HIV-1-infected individuals. The functionality of these responses was underscored by the consistent emergence of escape mutations in multiple CD8⁺ T cell epitopes during the course of infection. Although prechallenge vaccine-induced responses were largely undetectable, the Gag immunization increased both the magnitude and the kinetics of anamnestic Gag-specific T cell responses following HIV-1 infection, and the magnitude of these postchallenge Gag-specific responses was inversely correlated with acute HIV-1 viremia. Indeed, Gag immunization was associated with a modest but significant 0.5-log reduction in HIV-1 viral load when analyzed across four experimental groups of BLT mice. Notably, the HSV vector induced elevated plasma concentrations of polarizing cytokines and chemotactic factors, including interleukin-12p70 (IL-12p70) and MIP-1α, which were positively correlated with the magnitude of Gag-specific responses. Overall, these results support the ability of BLT mice to recapitulate human pathogen-specific T cell responses and to respond to immunization; however, additional improvements to the model are required to develop a robust system for testing HIV-1 vaccine efficacy.IMPORTANCE Advances in the development of humanized mice have raised the possibility of a small-animal model for preclinical testing of an HIV-1 vaccine. Here, we describe the capacity of BLT humanized mice to mount broadly directed HIV-1-specific human T cell responses that are functionally active, as indicated by the rapid emergence of viral escape mutations. Although immunization of BLT mice with the conserved viral Gag protein did not result in detectable prechallenge responses, it did increase the magnitude and kinetics of postchallenge Gag-specific T cell responses, which was associated with a modest but significant reduction in acute HIV-1 viremia. Additionally, the BLT model revealed immunization-associated increases in the plasma concentrations of immunomodulatory cytokines and chemokines that correlated with more robust T cell responses. These data support the potential utility of the BLT humanized mouse for HIV-1 vaccine development but suggest that additional improvements to the model are warranted.

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.

Cynomolgus macaque IL37 polymorphism and control of SIV infection

The association between gene polymorphisms and plasma virus load at the set point (SP-PVL) was investigated in Mauritian macaques inoculated with SIV. Among 44 macaques inoculated with 50 AID50, six individuals were selected: three with SP-PVL among the highest and three with SP-PVL among the lowest. The exons of 390 candidate genes of these six animals were sequenced. Twelve non-synonymous single nucleotide polymorphisms (NS-SNPs) lying in nine genes potentially associated with PVL were genotyped in 23 animals. Three NS-SNPs with probabilities of association with PVL less than 0.05 were genotyped in a total of 44 animals. One NS-SNP lying in exon 1 of the IL37 gene displayed a significant association (p = 3.33 × 10⁻⁴) and a strong odds ratio (19.52). Multiple linear regression modeling revealed three significant predictors of SP-PVL, including the IL37 exon 1 NS-SNP (p = 0.0004) and the MHC Class IB haplotypes M2 (p = 0.0007) and M6 (p = 0.0013). These three factors in conjunction explained 48% of the PVL variance (p = 4.8 × 10⁻⁶). The potential role of IL37 in the control of SIV infection is discussed.

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.

Preexisting Simian Immunodeficiency Virus Infection Increases Susceptibility to Tuberculosis in Mauritian Cynomolgus Macaques

Tuberculosis (TB), caused by Mycobacterium tuberculosis, is the leading cause of death among human immunodeficiency virus (HIV)-positive patients. The precise mechanisms by which HIV impairs host resistance to a subsequent M. tuberculosis infection are unknown. We modeled this coinfection in Mauritian cynomolgus macaques (MCM) using simian immunodeficiency virus (SIV) as an HIV surrogate. We infected seven MCM with SIVmac239 intrarectally and 6 months later coinfected them via bronchoscope with ∼10 CFU of M. tuberculosis Another eight MCM were infected with M. tuberculosis alone. TB progression was monitored by clinical parameters, by culturing bacilli in gastric and bronchoalveolar lavages, and by serial [¹⁸F]fluorodeoxyglucose (FDG) positron emission tomography/computed tomography (PET/CT) imaging. The eight MCM infected with M. tuberculosis alone displayed dichotomous susceptibility to TB, with four animals reaching humane endpoint within 13 weeks and four animals surviving >19 weeks after M. tuberculosis infection. In stark contrast, all seven SIV⁺ animals exhibited rapidly progressive TB following coinfection and all reached humane endpoint by 13 weeks. Serial PET/CT imaging confirmed dichotomous outcomes in MCM infected with M. tuberculosis alone and marked susceptibility to TB in all SIV⁺ MCM. Notably, imaging revealed a significant increase in TB granulomas between 4 and 8 weeks after M. tuberculosis infection in SIV⁺ but not in SIV-naive MCM and implies that SIV impairs the ability of animals to contain M. tuberculosis dissemination. At necropsy, animals with preexisting SIV infection had more overall pathology, increased bacterial loads, and a trend towards more extrapulmonary disease than animals infected with M. tuberculosis alone. We thus developed a tractable MCM model in which to study SIV-M. tuberculosis coinfection and demonstrate that preexisting SIV dramatically diminishes the ability to control M. tuberculosis coinfection.

Whole Exome Sequencing of HIV-1 long-term non-progressors identifies rare variants in genes encoding innate immune sensors and signaling molecules

Common CCR5-∆32 and HLA alleles only explain a minority of the HIV long-term non-progressor (LTNP) and elite controller (EC) phenotypes. To identify rare genetic variants contributing to the slow disease progression phenotypes, we performed whole exome sequencing (WES) on seven LTNPs and four ECs. HLA and CCR5 allele status, total HIV DNA reservoir size, as well as variant-related functional differences between the ECs, LTNPs, and eleven age- and gender-matched HIV-infected non-controllers on antiretroviral therapy (NCARTs) were investigated. Several rare variants were identified in genes involved in innate immune sensing, CD4-dependent infectivity, HIV trafficking, and HIV transcription mainly within the LTNP group. ECs and LTNPs had a significantly lower HIV reservoir compared to NCARTs. Furthermore, three LTNPs with variants affecting HIV nuclear import showed integrated HIV DNA levels below detection limit after in vitro infection. HIV slow progressors with variants in the TLR and NOD2 pathways showed reduced pro-inflammatory responses compared to matched controls. Low-range plasma levels of fibronectin was observed in a LTNP harboring two FN1 variants. Taken together, this study identified rare variants in LTNPs as well as in one EC, which may contribute to understanding of HIV pathogenesis and these slow progressor phenotypes, especially in individuals without protecting CCR5-∆32 and HLA alleles.

Acute-Phase CD4+ T Cell Responses Targeting Invariant Viral Regions Are Associated with Control of Live Attenuated Simian Immunodeficiency Virus

We manipulated SIVmac239Δnef, a model of major histocompatibility complex (MHC)-independent viral control, to evaluate characteristics of effective cellular responses mounted by Mauritian cynomolgus macaques (MCMs) that express the M3 MHC haplotype, which has been associated with poor control of pathogenic simian immunodeficiency virus (SIV). We created SIVΔnef-8x to test the hypothesis that effective SIV-specific T cell responses targeting invariant viral regions can emerge in the absence of immunodominant CD8⁺ T cell responses targeting variable epitopes and that control is achievable in individuals lacking known "protective" MHC alleles. Full-proteome gamma interferon (IFN-γ) enzyme-linked immunospot (ELISPOT) assays identified six newly targeted immunogenic regions following SIVΔnef-8x infection of M3/M3 MCMs. We deep sequenced circulating virus and found that four of the six newly targeted regions rarely accumulated mutations. Six animals infected with SIVΔnef-8x had T cell responses that targeted at least one of the four invariant regions and had a lower set point viral load than two animals that did not have T cell responses that targeted any invariant regions. We found that MHC class II molecules restricted all four of the invariant peptide regions, while the two variable regions were restricted by MHC class I molecules. Therefore, in the absence of immunodominant CD8⁺ T cell responses that target variable regions during SIVmac239Δnef infection, individuals without protective MHC alleles developed predominantly CD4⁺ T cell responses specific for invariant regions that may improve control of virus replication. Our results provide some evidence that antiviral CD4⁺ T cells during acute SIV infection can contribute to effective viral control and should be considered in strategies to combat HIV infection.IMPORTANCE Studies defining effective cellular immune responses to human immunodeficiency virus (HIV) and SIV have largely focused on a rare population that express specific MHC class I alleles and control virus replication in the absence of antiretroviral treatment. This leaves in question whether similar effective immune responses can be achieved in the larger population. The majority of HIV-infected individuals mount CD8⁺ T cell responses that target variable viral regions that accumulate high-frequency escape mutations. Limiting T cell responses to these variable regions and targeting invariant viral regions, similar to observations in rare "elite controllers," may provide an ideal strategy for the development of effective T cell responses in individuals with diverse MHC genetics. Therefore, it is of paramount importance to determine whether T cell responses can be redirected toward invariant viral regions in individuals without protective MHC alleles and if these responses improve control of virus replication.

Vaccination with the Conserved Caveolin-1 Binding Motif in Human Immunodeficiency Virus Type 1 Glycoprotein gp41 Delays the Onset of Viral Infection and Provides Partial Protection in Simian/Human Immunodeficiency Virus-Challenged Cynomolgus Macaques

We have previously reported that the CBD1 peptide (SLEQIWNNMTWMQWDK), corresponding to the consensus caveolin-1 binding domain in human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein gp41, elicits peptide-specific antibodies. Here, we have investigated the cellular immune response and the protective efficacy against a simian/human immunodeficiency virus (SHIV_162P3) challenge. In addition to the CBD1 peptide, peptides overlapping the caveolin-binding-motif (CBM) (₆₂₂IWNNMTWMQW₆₃₁ or ₆₂₂IWNNMTW₆₂₈) were fused to a Gag-p24 T helper epitope for vaccination. All immunized cynomolgus macaques responded to a cocktail peptide immunization by inducing specific T cells and the production of high-titer CBD1/CBM peptide-specific antibodies. Six months after the fourth vaccine boost, six control and five vaccinated animals were challenged weekly by repeated exposure to SHIV_162P3 via the mucosal rectal route. All control animals were infected after 1 to 3 challenges with SHIV, while among the five vaccinated monkeys, three became infected after a delay compared to control; one was infected after the eighth viral challenge, and one remained uninfected even after the ninth SHIV challenge. Immunized animals maintained a CD4 T cell count, and their central memory CD4 T cells were less depleted than in the control group. Furthermore, SHIV challenge stimulates antigen-specific memory T cell response in vaccinated macaques. Our results indicate that peptides derived from the CBM region can be immunogenic and provide protection against SHIV infection in cynomolgus monkeys.IMPORTANCE In HIV-1-producing cells, gp41 exists in a complexed form with caveolin-1, an interaction most probably mediated by the caveolin-1 binding motif. This sequence is highly conserved in every single HIV-1 isolate, thus suggesting that there is constant selective pressure to preserve this sequence for a specific function in the HIV infectious cycle. Consequently, the CBM sequence may represent the "Achilles' heel" of HIV-1 in the development of an efficient vaccine. Our results demonstrate that macaques immunized with the CBM-based peptides displayed a delay in the onset of viral infection and CD4 depletion, as well as a significant induction of antigen-specific memory T cell response, which is essential for the control of HIV/SIV infections. Finally, as HIV-infected individuals lack anti-CBM immune responses, CBM-based vaccines could have applications as a therapeutic vaccine in AIDS patients.

Nonhuman primate models of human viral infections

Humans have a close phylogenetic relationship with nonhuman primates (NHPs) and share many physiological parallels, such as highly similar immune systems, with them. Importantly, NHPs can be infected with many human or related simian viruses. In many cases, viruses replicate in the same cell types as in humans, and infections are often associated with the same pathologies. In addition, many reagents that are used to study the human immune response cross-react with NHP molecules. As such, NHPs are often used as models to study viral vaccine efficacy and antiviral therapeutic safety and efficacy and to understand aspects of viral pathogenesis. With several emerging viral infections becoming epidemic, NHPs are proving to be a very beneficial benchmark for investigating human viral infections.

Whole genome sequencing in the search for genes associated with the control of SIV infection in the Mauritian macaque model

In the Mauritian macaque experimentally inoculated with SIV, gene polymorphisms potentially associated with the plasma virus load at a set point, approximately 100 days post inoculation, were investigated. Among the 42 animals inoculated with 50 AID50 of the same strain of SIV, none of which received any preventive or curative treatment, nine individuals were selected: three with a plasma virus load (PVL) among the lowest, three with intermediate PVL values and three among the highest PVL values. The complete genomes of these nine animals were then analyzed. Initially, attention was focused on variants with a potential functional impact on protein encoding genes (non-synonymous SNPs (NS-SNPs) and splicing variants). Thus, 424 NS-SNPs possibly associated with PVL were detected. The 424 candidates SNPs were genotyped in these 42 SIV experimentally infected animals (including the nine animals subjected to whole genome sequencing). The genes containing variants most probably associated with PVL at a set time point are analyzed herein.

Mauritian cynomolgus macaques with M3M4 MHC genotype control SIVmac251 infection

BACKGROUND

Understanding natural HIV control may lead to new preventative or therapeutic strategies. Several protective major histocompatibility complex (MHC) genotypes were found in humans and rhesus macaques. Here, we report a simian immunodeficiency virus (SIV) controller MHC genotype in Mauritian cynomolgus macaques (MCMs).

METHODS

Twelve MHC-genotyped MCMs were infected with SIVmac251 and monitored for viral loads and CD4+ T-cell counts.

RESULTS

Two macaques with M3M4 genotype exhibited the lowest peak viral loads (log plasma SIV RNA copies/mL), nearly 3 logs lower than those in most macaques with other MHC haplotype combinations, and set point viral loads below the level of detection limit by RT-qPCR (<2 log RNA copies/mL). They maintained healthy CD4+ T-cell counts of >500 cells/μL blood, while CD4 counts in the vast majority of other macaques were below this level.

CONCLUSIONS

The M3M4 MHC genotype may confer enhanced control of SIV replication in MCMs.

Development of genotyping method for functionally relevant variants of cytochromes P450 in cynomolgus macaques

In cynomolgus macaques (Macaca fascicularis), widely used in drug metabolism studies, CYP2C9, CYP2C76, CYP2D6, CYP3A4, and CYP3A5, important drug-metabolizing enzymes, are abundantly expressed in liver and metabolize cytochrome P450 substrates. CYP2C9 (c.334A>C), CYP2C76 (c.449TG>A), CYP2D6 (c.891A>G), CYP3A4 (IVS3 + 1G>del), and CYP3A5 (c.625A>T) substantially influence metabolic activity of enzymes, and thus are important variants in drug metabolism studies. In this study, a real-time PCR method was developed for genotyping these variants. The validity of the methods was verified by genotyping two wild type, two heterozygous, and two homozygous DNAs and was used to genotype 41 cynomolgus macaques (from Cambodia, Indonesia, the Philippines, or Vietnam) for the five variants, along with another important variant CYP2C19 (c.308C>T). The CYP2C9 and CYP2C19 variants were found only in Cambodian and Vietnamese animals, while the CYP2C76 and CYP2D6 variants were found only in Indonesian and Philippine animals. The CYP3A4 and CYP3A5 variants were not found in any of the animals analyzed. Mauritian animals, genotyped using next-generation sequencing data for comparison, possessed the CYP2C19 and CYP2D6 variants, but not the other variants. These results indicated differences in prevalence of these important variants among animal groups. Therefore, the genotyping tool developed is useful for drug metabolism studies using cynomolgus macaques.

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

Killer cell immunoglobulin-like receptors (KIRs) modulate disease progression of pathogens including HIV, malaria, and hepatitis C. Cynomolgus and rhesus macaques are widely used as nonhuman primate models to study human pathogens, and so, considerable effort has been put into characterizing their KIR genetics. However, previous studies have relied on cDNA cloning and Sanger sequencing that lack the throughput of current sequencing platforms. In this study, we present a high throughput, full-length allele discovery method utilizing Pacific Biosciences circular consensus sequencing (CCS). We also describe a new approach to Macaque Exome Sequencing (MES) and the development of the Rhexome1.0, an adapted target capture reagent that includes macaque-specific capture probe sets. By using sequence reads generated by whole genome sequencing (WGS) and MES to inform primer design, we were able to increase the sensitivity of KIR allele discovery. We demonstrate this increased sensitivity by defining nine novel alleles within a cohort of Mauritian cynomolgus macaques (MCM), a geographically isolated population with restricted KIR genetics that was thought to be completely characterized. Finally, we describe an approach to genotyping KIRs directly from sequence reads generated using WGS/MES reads. The findings presented here expand our understanding of KIR genetics in MCM by associating new genes with all eight KIR haplotypes and demonstrating the existence of at least one KIR3DS gene associated with every haplotype.

Diversity and regulatory impact of copy number variation in the primate Macaca fascicularis

BACKGROUND

Copy number variations (CNVs) are a significant source of genetic diversity and commonly found in mammalian genomes. We have generated a genome-wide CNV map for Cynomolgus monkeys (Macaca fascicularis). This crab-eating macaque is the closest animal model to humans that is used in biomedical research.

RESULTS

We show that Cynomolgus monkey CNVs are in general much smaller in size than gene loci and are specific to the population of origin. Genome-wide expression data from five vitally important organs demonstrates that CNVs in close proximity to transcription start sites associate strongly with expression changes. Among these eQTL genes we find an overrepresentation of genes involved in metabolism, receptor activity, and transcription.

CONCLUSION

These results provide evidence that CNVs shape tissue transcriptomes in monkey populations, potentially offering an adaptive advantage. We suggest that this genetic diversity should be taken into account when using Cynomolgus macaques as models.

Discovery of novel MHC-class I alleles and haplotypes in Filipino cynomolgus macaques (Macaca fascicularis) by pyrosequencing and Sanger sequencing: Mafa-class I polymorphism

Although the low polymorphism of the major histocompatibility complex (MHC) transplantation genes in the Filipino cynomolgus macaque (Macaca fascicularis) is expected to have important implications in the selection and breeding of animals for medical research, detailed polymorphism information is still lacking for many of the duplicated class I genes. To better elucidate the degree and types of MHC polymorphisms and haplotypes in the Filipino macaque population, we genotyped 127 unrelated animals by the Sanger sequencing method and high-resolution pyrosequencing and identified 112 different alleles, 28 at cynomolgus macaque MHC (Mafa)-A, 54 at Mafa-B, 12 at Mafa-I, 11 at Mafa-E, and seven at Mafa-F alleles, of which 56 were newly described. Of them, the newly discovered Mafa-A8*01:01 lineage allele had low nucleotide similarities (<86%) with primate MHC class I genes, and it was also conserved in the Vietnamese and Indonesian populations. In addition, haplotype estimations revealed 17 Mafa-A, 23 Mafa-B, and 12 Mafa-E haplotypes integrated with 84 Mafa-class I haplotypes and Mafa-F alleles. Of these, the two Mafa-class I haplotypes, F/A/E/B-Hp1 and F/A/E/B-Hp2, had the highest haplotype frequencies at 10.6 and 10.2%, respectively. This suggests that large scale genetic screening of the Filipino macaque population would identify these and other high-frequency Mafa-class I haplotypes that could be used as MHC control animals for the benefit of biomedical research.

Whole-genome sequencing of six Mauritian Cynomolgus macaques (Macaca fascicularis) reveals a genome-wide pattern of polymorphisms under extreme population bottleneck

Cynomolgus macaques (Macaca fascicularis) were introduced to the island of Mauritius by humans around the 16th century. The unique demographic history of the Mauritian cynomolgus macaques provides the opportunity to not only examine the genetic background of well-established nonhuman primates for biomedical research but also understand the effect of an extreme population bottleneck on the pattern of polymorphisms in genomes. We sequenced the whole genomes of six Mauritian cynomolgus macaques and obtained an average of 20-fold coverage of the genome sequences for each individual. The overall level of nucleotide diversity was 23% smaller than that of the Malaysian cynomolgus macaques, and a reduction of low-frequency polymorphisms was observed. In addition, we also confirmed that the Mauritian cynomolgus macaques were genetically closer to a representative of the Malaysian population than to a representative of the Indochinese population. Excess of nonsynonymous polymorphisms in low frequency, which has been observed in many other species, was not very strong in the Mauritian samples, and the proportion of heterozygous nonsynonymous polymorphisms relative to synonymous polymorphisms is higher within individuals in Mauritian than Malaysian cynomolgus macaques. Those patterns indicate that the extreme population bottleneck made purifying selection overwhelmed by the power of genetic drift in the population. Finally, we estimated the number of founding individuals by using the genome-wide site frequency spectrum of the six samples. Assuming a simple demographic scenario with a single bottleneck followed by exponential growth, the estimated number of founders (∼20 individuals) is largely consistent with previous estimates.

Bringing non-human primate research into the post-genomic era: how monkeys are teaching us about elite controllers of HIV/AIDS

Whole-genome sequencing of Mauritian cynomolgus macaques reveals novel candidate loci for controlling simian immunodeficiency virus replication.

See related Research, http://genomebiology.com/2014/15/11/478