ABSTRACT
Allogeneic hematopoietic stem cell transplantation (alloHSCT) from donors lacking C-C chemokine receptor 5 (CCR5Δ32/Δ32) can cure HIV, yet mechanisms remain speculative. To define how alloHSCT mediates HIV cure, we performed MHC-matched alloHSCT in SIV+, anti-retroviral therapy (ART)-suppressed Mauritian cynomolgus macaques (MCMs) and demonstrated that allogeneic immunity was the major driver of reservoir clearance, occurring first in peripheral blood, then peripheral lymph nodes, and finally in mesenteric lymph nodes draining the gastrointestinal tract. While allogeneic immunity could extirpate the latent viral reservoir and did so in two alloHSCT-recipient MCMs that remained aviremic >2.5 years after stopping ART, in other cases, it was insufficient without protection of engrafting cells afforded by CCR5-deficiency, as CCR5-tropic virus spread to donor CD4+ T cells despite full ART suppression. These data demonstrate the individual contributions of allogeneic immunity and CCR5 deficiency to HIV cure and support defining targets of alloimmunity for curative strategies independent of HSCT.
Subject(s)
HIV Infections , Hematopoietic Stem Cell Transplantation , Simian Acquired Immunodeficiency Syndrome , Simian Immunodeficiency Virus , Animals , Macaca fascicularis , Viral LoadABSTRACT
BACKGROUND: Naturally occurring colorectal cancers (CRC) in rhesus macaques share many features with their human counterparts and are useful models for cancer immunotherapy; but mechanistic data are lacking regarding the comparative molecular pathogenesis of these cancers. METHODS: We conducted state-of-the-art imaging including CT and PET, clinical assessments, and pathological review of 24 rhesus macaques with naturally occurring CRC. Additionally, we molecularly characterized these tumors utilizing immunohistochemistry (IHC), microsatellite instability assays, DNAseq, transcriptomics, and developed a DNA methylation-specific qPCR assay for MLH1, CACNA1G, CDKN2A, CRABP1, and NEUROG1, human markers for CpG island methylator phenotype (CIMP). We furthermore employed Monte-Carlo simulations to in-silico model alterations in DNA topology in transcription-factor binding site-rich promoter regions upon experimentally demonstrated DNA methylation. RESULTS: Similar cancer histology, progression patterns, and co-morbidities could be observed in rhesus as reported for human CRC patients. IHC identified loss of MLH1 and PMS2 in all cases, with functional microsatellite instability. DNA sequencing revealed the close genetic relatedness to human CRCs, including a similar mutational signature, chromosomal instability, and functionally-relevant mutations affecting KRAS (G12D), TP53 (R175H, R273*), APC, AMER1, ALK, and ARID1A. Interestingly, MLH1 mutations were rarely identified on a somatic or germline level. Transcriptomics not only corroborated the similarities of rhesus and human CRCs, but also demonstrated the significant downregulation of MLH1 but not MSH2, MSH6, or PMS2 in rhesus CRCs. Methylation-specific qPCR suggested CIMP-positivity in 9/16 rhesus CRCs, but all 16/16 exhibited significant MLH1 promoter hypermethylation. DNA hypermethylation was modelled to affect DNA topology, particularly propeller twist and roll profiles. Modelling the DNA topology of a transcription factor binding motif (TFAP2A) in the MLH1 promoter that overlapped with a methylation-specific probe, we observed significant differences in DNA topology upon experimentally shown DNA methylation. This suggests a role of transcription factor binding interference in epigenetic silencing of MLH1 in rhesus CRCs. CONCLUSIONS: These data indicate that epigenetic silencing suppresses MLH1 transcription, induces the loss of MLH1 protein, abrogates mismatch repair, and drives genomic instability in naturally occurring CRC in rhesus macaques. We consider this spontaneous, uninduced CRC in immunocompetent, treatment-naïve rhesus macaques to be a uniquely informative model for human CRC.
Subject(s)
Brain Neoplasms , Colorectal Neoplasms , Microsatellite Instability , Neoplastic Syndromes, Hereditary , Humans , Animals , Macaca mulatta/genetics , Macaca mulatta/metabolism , MutL Protein Homolog 1/genetics , Mismatch Repair Endonuclease PMS2/genetics , Mismatch Repair Endonuclease PMS2/metabolism , Colorectal Neoplasms/pathology , DNA Methylation/genetics , Epigenesis, Genetic , Transcription Factors/genetics , Transcription Factors/metabolism , DNA/metabolism , DNA Mismatch Repair/geneticsABSTRACT
The CCR5-specific antibody Leronlimab is being investigated as a novel immunotherapy that can suppress HIV replication with minimal side effects. Here we studied the virological and immunological consequences of Leronlimab in chronically CCR5-tropic HIV-1 infected humans (n = 5) on suppressive antiretroviral therapy (ART) and in ART-naïve acutely CCR5-tropic SHIV infected rhesus macaques (n = 4). All five human participants transitioned from daily combination ART to self-administered weekly subcutaneous (SC) injections of 350 mg or 700 mg Leronlimab and to date all participants have sustained virologic suppression for over seven years. In all participants, Leronlimab fully occupied CCR5 receptors on peripheral blood CD4+ T cells and monocytes. In ART-naïve rhesus macaques acutely infected with CCR5-tropic SHIV, weekly SC injections of 50 mg/kg Leronlimab fully suppressed plasma viremia in half of the macaques. CCR5 receptor occupancy by Leronlimab occurred concomitant with rebound of CD4+ CCR5+ T-cells in peripheral blood, and full CCR5 receptor occupancy was found in multiple anatomical compartments. Our results demonstrate that weekly, self-administered Leronlimab was safe, well-tolerated, and efficacious for long-term virologic suppression and should be included in the arsenal of safe, easily administered, longer-acting antiretroviral treatments for people living with HIV-1. Trial Registration: ClinicalTrials.gov Identifiers: NCT02175680 and NCT02355184.
Subject(s)
Simian Immunodeficiency Virus , Animals , Antibodies, Monoclonal, Humanized/pharmacology , HIV Antibodies , Humans , Macaca mulatta , Receptors, CCR5ABSTRACT
MOTIVATION: Single-cell sequencing methods provide previously impossible resolution into the transcriptome of individual cells. Cell hashing reduces single-cell sequencing costs by increasing capacity on droplet-based platforms. Cell hashing methods rely on demultiplexing algorithms to accurately classify droplets; however, assumptions underlying these algorithms limit accuracy of demultiplexing, ultimately impacting the quality of single-cell sequencing analyses. RESULTS: We present Bimodal Flexible Fitting (BFF) demultiplexing algorithms BFFcluster and BFFraw, a novel class of algorithms that rely on the single inviolable assumption that barcode count distributions are bimodal. We integrated these and other algorithms into cellhashR, a new R package that provides integrated QC and a single command to execute and compare multiple demultiplexing algorithms. We demonstrate that BFFcluster demultiplexing is both tunable and insensitive to issues with poorly behaved data that can confound other algorithms. Using two well-characterized reference datasets, we demonstrate that demultiplexing with BFF algorithms is accurate and consistent for both well-behaved and poorly behaved input data. AVAILABILITY AND IMPLEMENTATION: cellhashR is available as an R package at https://github.com/BimberLab/cellhashR. cellhashR version 1.0.3 was used for the analyses in this manuscript and is archived on Zenodo at https://www.doi.org/10.5281/zenodo.6402477. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.
Subject(s)
Algorithms , Software , Electronic Data Processing , Sequence Analysis , Single-Cell AnalysisABSTRACT
CD8+ T cells are key mediators of antiviral and antitumor immunity. The isolation and study of Ag-specific CD8+ T cells, as well as mapping of their MHC restriction, has practical importance to the study of disease and the development of therapeutics. Unfortunately, most experimental approaches are cumbersome, owing to the highly variable and donor-specific nature of MHC-bound peptide/TCR interactions. Here we present a novel system for rapid identification and characterization of Ag-specific CD8+ T cells, particularly well suited for samples with limited primary cells. Cells are stimulated ex vivo with Ag of interest, followed by live cell sorting based on surface-trapped TNF-α. We take advantage of major advances in single-cell sequencing to generate full-length sequence data from the paired TCR α- and ß-chains from these Ag-specific cells. The paired TCR chains are cloned into retroviral vectors and used to transduce donor CD8+ T cells. These TCR transductants provide a virtually unlimited experimental reagent, which can be used for further characterization, such as minimal epitope mapping or identification of MHC restriction, without depleting primary cells. We validated this system using CMV-specific CD8+ T cells from rhesus macaques, characterizing an immunodominant Mamu-A1*002:01-restricted epitope. We further demonstrated the utility of this system by mapping a novel HLA-A*68:02-restricted HIV Gag epitope from an HIV-infected donor. Collectively, these data validate a new strategy to rapidly identify novel Ags and characterize Ag-specific CD8+ T cells, with applications ranging from the study of infectious disease to immunotherapeutics and precision medicine.
Subject(s)
CD8-Positive T-Lymphocytes , HIV Infections , Animals , Epitopes , Epitopes, T-Lymphocyte , Macaca mulatta , Receptors, Antigen, T-Cell , Tumor Necrosis Factor-alphaABSTRACT
CRISPR/Cas systems are some of the most promising tools for therapeutic genome editing. The use of these systems is contingent on the optimal designs of guides and homology-directed repair (HDR) templates. While this design can be achieved in silico, validation and further optimization are usually performed with the help of reporter systems. Here, we describe a novel reporter system, termed BETLE, that allows for the fast, sensitive, and cell-specific detection of genome editing and template-specific HDR by encoding multiple reporter proteins in different open-reading frames. Out-of-frame non-homologous end joining (NHEJ) leads to the expression of either secretable NanoLuc luciferase, enabling a highly sensitive and low-cost analysis of editing, or fluorescent mTagBFP2, allowing for the enumeration and tissue-specific localization of genome-edited cells. BETLE includes a site to validate CRISPR/Cas systems for a sequence-of-interest, making it broadly adaptable. We evaluated BETLE using a defective moxGFP with a 39-base-pair deletion and showed spCas9, saCas9, and asCas12a editing as well as sequence-specific HDR and the repair of moxGFP in cell lines with single and multiple reporter integrants. Taken together, these data show that BETLE allows for the rapid detection and optimization of CRISPR/Cas genome editing and HDR in vitro and represents a state-of the art tool for future applications in vivo.
Subject(s)
CRISPR-Cas Systems , DNA Breaks, Double-Stranded , CRISPR-Cas Systems/genetics , Gene Editing , DNA End-Joining Repair , GenomeABSTRACT
Cytomegaloviruses (CMVs) are highly adapted to their host species resulting in strict species specificity. Hence, in vivo examination of all aspects of CMV biology employs animal models using host-specific CMVs. Infection of rhesus macaques (RM) with rhesus CMV (RhCMV) has been established as a representative model for infection of humans with HCMV due to the close evolutionary relationships of both host and virus. However, the only available RhCMV clone that permits genetic modifications is based on the 68-1 strain which has been passaged in fibroblasts for decades resulting in multiple genomic changes due to tissue culture adaptations. As a result, 68-1 displays reduced viremia in RhCMV-naïve animals and limited shedding compared to non-clonal, low passage isolates. To overcome this limitation, we used sequence information from primary RhCMV isolates to construct a full-length (FL) RhCMV by repairing all mutations affecting open reading frames (ORFs) in the 68-1 bacterial artificial chromosome (BAC). Inoculation of adult, immunocompetent, RhCMV-naïve RM with the reconstituted virus resulted in significant viremia in the blood similar to primary isolates of RhCMV and furthermore led to high viral genome copy numbers in many tissues at day 14 post infection. In contrast, viral dissemination was greatly reduced upon deletion of genes also lacking in 68-1. Transcriptome analysis of infected tissues further revealed that chemokine-like genes deleted in 68-1 are among the most highly expressed viral transcripts both in vitro and in vivo consistent with an important immunomodulatory function of the respective proteins. We conclude that FL-RhCMV displays in vitro and in vivo characteristics of a wildtype virus while being amenable to genetic modifications through BAC recombineering techniques.
Subject(s)
Cytomegalovirus Infections/virology , Cytomegalovirus/genetics , Genome, Viral/genetics , Viremia , Animals , Cell Line , Chromosomes, Artificial, Bacterial , Cytomegalovirus/pathogenicity , DNA, Recombinant , Disease Models, Animal , Female , Fibroblasts/virology , Humans , Macaca mulatta , Male , Mutation , Open Reading Frames/genetics , Phylogeny , Species SpecificityABSTRACT
SUMMARY: Large scale genomic studies produce millions of sequence variants, generating datasets far too massive for manual inspection. To ensure variant and genotype data are consistent and accurate, it is necessary to evaluate variants prior to downstream analysis using quality control (QC) reports. Variant call format (VCF) files are the standard format for representing variant data; however, generating summary statistics from these files is not always straightforward. While tools to summarize variant data exist, they generally produce simple text file tables, which still require additional processing and interpretation. VariantQC fills this gap as a user friendly, interactive visual QC report that generates and concisely summarizes statistics from VCF files. The report aggregates and summarizes variants by dataset, chromosome, sample and filter type. The VariantQC report is useful for high-level dataset summary, quality control and helps flag outliers. Furthermore, VariantQC operates on VCF files, so it can be easily integrated into many existing variant pipelines. AVAILABILITY AND IMPLEMENTATION: DISCVRSeq's VariantQC tool is freely available as a Java program, with the compiled JAR and source code available from https://github.com/BimberLab/DISCVRSeq/. Documentation and example reports are available at https://bimberlab.github.io/DISCVRSeq/.
Subject(s)
Software , Genetic Variation , Genomics , Genotype , Quality ControlABSTRACT
Allogeneic hematopoietic stem cell transplantation (HSCT) and xenotransplantation are accompanied by viral reactivations and virus-associated complications resulting from immune deficiency. Here, in a Mauritian cynomolgus macaque model of fully MHC-matched allogeneic HSCT, we report reactivations of cynomolgus polyomavirus, lymphocryptovirus, and cytomegalovirus, macaque viruses analogous to HSCT-associated human counterparts BK virus, Epstein-Barr virus, and human cytomegalovirus. Viral replication in recipient macaques resulted in characteristic disease manifestations observed in HSCT patients, such as polyomavirus-associated hemorrhagic cystitis and tubulointerstitial nephritis or lymphocryptovirus-associated post-transplant lymphoproliferative disorder. However, in most cases, the reconstituted immune system, alone or in combination with short-term pharmacological intervention, exerted control over viral replication, suggesting engraftment of functional donor-derived immunity. Indeed, the donor-derived reconstituted immune systems of two long-term engrafted HSCT recipient macaques responded to live attenuated yellow fever 17D vaccine (YFV 17D) indistinguishably from untransplanted controls, mounting 17D-targeted neutralizing antibody responses and clearing YFV 17D within 14 days. Together, these data demonstrate that this macaque model of allogeneic HSCT recapitulates clinical situations of opportunistic viral infections in transplant patients and provides a pre-clinical model to test novel prophylactic and therapeutic modalities.
Subject(s)
Disease Models, Animal , Hematopoietic Stem Cell Transplantation , Opportunistic Infections , Virus Diseases , Allografts , Animals , Hematopoietic Stem Cell Transplantation/adverse effects , Humans , Macaca fascicularis , Opportunistic Infections/virologyABSTRACT
MHC-E is a highly conserved nonclassical MHC class Ib molecule that predominantly binds and presents MHC class Ia leader sequence-derived peptides for NK cell regulation. However, MHC-E also binds pathogen-derived peptide Ags for presentation to CD8+ T cells. Given this role in adaptive immunity and its highly monomorphic nature in the human population, HLA-E is an attractive target for novel vaccine and immunotherapeutic modalities. Development of HLA-E-targeted therapies will require a physiologically relevant animal model that recapitulates HLA-E-restricted T cell biology. In this study, we investigated MHC-E immunobiology in two common nonhuman primate species, Indian-origin rhesus macaques (RM) and Mauritian-origin cynomolgus macaques (MCM). Compared to humans and MCM, RM expressed a greater number of MHC-E alleles at both the population and individual level. Despite this difference, human, RM, and MCM MHC-E molecules were expressed at similar levels across immune cell subsets, equivalently upregulated by viral pathogens, and bound and presented identical peptides to CD8+ T cells. Indeed, SIV-specific, Mamu-E-restricted CD8+ T cells from RM recognized antigenic peptides presented by all MHC-E molecules tested, including cross-species recognition of human and MCM SIV-infected CD4+ T cells. Thus, MHC-E is functionally conserved among humans, RM, and MCM, and both RM and MCM represent physiologically relevant animal models of HLA-E-restricted T cell immunobiology.
Subject(s)
CD4-Positive T-Lymphocytes/immunology , CD8-Positive T-Lymphocytes/immunology , Histocompatibility Antigens Class I/metabolism , Histocompatibility Antigens/metabolism , Killer Cells, Natural/immunology , Simian Acquired Immunodeficiency Syndrome/immunology , Simian Immunodeficiency Virus/immunology , Animals , Antigen Presentation , Antigens, Viral/immunology , Antigens, Viral/metabolism , Cells, Cultured , Conserved Sequence/genetics , Histocompatibility Antigens/genetics , Histocompatibility Antigens/immunology , Histocompatibility Antigens Class I/genetics , Histocompatibility Antigens Class I/immunology , Humans , Macaca fascicularis , Macaca mulatta , Models, Animal , Peptides/immunology , Peptides/metabolism , HLA-E AntigensABSTRACT
BACKGROUND: Non-human primates (NHPs), particularly macaques, serve as critical and highly relevant pre-clinical models of human disease. The similarity in human and macaque natural disease susceptibility, along with parallel genetic risk alleles, underscores the value of macaques in the development of effective treatment strategies. Nonetheless, there are limited genomic resources available to support the exploration and discovery of macaque models of inherited disease. Notably, there are few public databases tailored to searching NHP sequence variants, and no other database making use of centralized variant calling, or providing genotype-level data and predicted pathogenic effects for each variant. RESULTS: The macaque Genotype And Phenotype (mGAP) resource is the first public website providing searchable, annotated macaque variant data. The mGAP resource includes a catalog of high confidence variants, derived from whole genome sequence (WGS). The current mGAP release at time of publication (1.7) contains 17,087,212 variants based on the sequence analysis of 293 rhesus macaques. A custom pipeline was developed to enable annotation of the macaque variants, leveraging human data sources that include regulatory elements (ENCODE, RegulomeDB), known disease- or phenotype-associated variants (GRASP), predicted impact (SIFT, PolyPhen2), and sequence conservation (Phylop, PhastCons). Currently mGAP includes 2767 variants that are identical to alleles listed in the human ClinVar database, of which 276 variants, spanning 258 genes, are identified as pathogenic. An additional 12,472 variants are predicted as high impact (SnpEff) and 13,129 are predicted as damaging (PolyPhen2). In total, these variants are predicted to be associated with more than 2000 human disease or phenotype entries reported in OMIM (Online Mendelian Inheritance in Man). Importantly, mGAP also provides genotype-level data for all subjects, allowing identification of specific individuals harboring alleles of interest. CONCLUSIONS: The mGAP resource provides variant and genotype data from hundreds of rhesus macaques, processed in a consistent manner across all subjects ( https://mgap.ohsu.edu ). Together with the extensive variant annotations, mGAP presents unprecedented opportunity to investigate potential genetic associations with currently characterized disease models, and to uncover new macaque models based on parallels with human risk alleles.
Subject(s)
Computational Biology/methods , Genetic Variation , Genotype , Phenotype , Animals , Disease Models, Animal , Humans , Information Storage and Retrieval , Internet , Macaca mulattaABSTRACT
Natural killer (NK) cells play a central role in immune responses through direct cytotoxicity and the release of cytokines that prime adaptive immunity. In simian primates, NK cell responses are regulated by interactions between two highly polymorphic sets of molecules: the killer-cell immunoglobulin-like receptors (KIRs) and their major histocompatibility complex (MHC) class I ligands. KIR-MHC class I interactions in humans have been implicated in the outcome of a number viral diseases and cancers. However, studies to address the role of KIRs in animal models have been limited by the complex immunogenetics and lack of defined ligands for KIRs in non-human primates. Due to the rapid evolution of KIRs, there is little conservation among the KIR genes of different primate species and it is not possible to predict the specificity of KIRs from known KIR-MHC class I interactions in humans. Hence, the MHC class I ligands for KIRs in species other than humans are poorly defined. Here, we review the KIR genes of the rhesus macaque, an important animal model for human immunodeficiency virus infection and other infectious diseases, and the MHC class I ligands that have been identified for KIRs in this species.
Subject(s)
Histocompatibility Antigens Class I/immunology , Killer Cells, Natural/immunology , Macaca/immunology , Receptors, KIR/immunology , Animals , Evolution, Molecular , Genetic Variation/genetics , Genetic Variation/immunology , Haplotypes , Histocompatibility Antigens Class I/genetics , Histocompatibility Antigens Class I/metabolism , Humans , Killer Cells, Natural/metabolism , Macaca/genetics , Macaca/metabolism , Protein Binding/immunology , Receptors, KIR/genetics , Receptors, KIR/metabolismABSTRACT
The increasing number of Killer Immunoglobulin-like Receptor (KIR) sequences available for non-human primate species and cattle has prompted development of a centralized database, guidelines for a standardized nomenclature, and minimum requirements for database submission. The guidelines and nomenclature are based on those used for human KIR and incorporate modifications made for inclusion of non-human species in the companion IPD-NHKIR database. Included in this first release are the rhesus macaque (Macaca mulatta), chimpanzee (Pan troglodytes), orangutan (Pongo abelii and Pongo pygmaeus), and cattle (Bos taurus).
Subject(s)
Receptors, KIR , Terminology as Topic , Animals , Cattle , Humans , Macaca mulatta/genetics , Pan troglodytes/genetics , Pongo pygmaeus/geneticsABSTRACT
Cytomegaloviruses (CMV) are highly species-specific due to millennia of co-evolution and adaptation to their host, with no successful experimental cross-species infection in primates reported to date. Accordingly, full genome phylogenetic analysis of multiple new CMV field isolates derived from two closely related nonhuman primate species, Indian-origin rhesus macaques (RM) and Mauritian-origin cynomolgus macaques (MCM), revealed distinct and tight lineage clustering according to the species of origin, with MCM CMV isolates mirroring the limited genetic diversity of their primate host that underwent a population bottleneck 400 years ago. Despite the ability of Rhesus CMV (RhCMV) laboratory strain 68-1 to replicate efficiently in MCM fibroblasts and potently inhibit antigen presentation to MCM T cells in vitro, RhCMV 68-1 failed to productively infect MCM in vivo, even in the absence of host CD8+ T and NK cells. In contrast, RhCMV clone 68-1.2, genetically repaired to express the homologues of the HCMV anti-apoptosis gene UL36 and epithelial cell tropism genes UL128 and UL130 absent in 68-1, efficiently infected MCM as evidenced by the induction of transgene-specific T cells and virus shedding. Recombinant variants of RhCMV 68-1 and 68-1.2 revealed that expression of either UL36 or UL128 together with UL130 enabled productive MCM infection, indicating that multiple layers of cross-species restriction operate even between closely related hosts. Cumulatively, these results implicate cell tropism and evasion of apoptosis as critical determinants of CMV transmission across primate species barriers, and extend the macaque model of human CMV infection and immunology to MCM, a nonhuman primate species with uniquely simplified host immunogenetics.
Subject(s)
Cytomegalovirus Infections/transmission , Cytomegalovirus/genetics , Disease Models, Animal , Macaca fascicularis/virology , Macaca mulatta/virology , Animals , Cytomegalovirus Infections/genetics , DNA, Viral/analysis , DNA, Viral/genetics , High-Throughput Nucleotide Sequencing , Phylogeny , Species SpecificityABSTRACT
Rhesus macaques are an important pre-clinical model of human disease. To advance our understanding of genomic variation that may influence disease, we surveyed genome-wide variation in 21 rhesus macaques. We employed best-practice variant calling, validated with Mendelian inheritance. Next, we used alignment data from our cohort to detect genomic regions likely to produce inaccurate genotypes, potentially due to either gene duplication or structural variation between individuals. We generated a final dataset of >16 million high confidence variants, including 13 million in Chinese-origin rhesus macaques, an increasingly important disease model. We detected an average of 131 mutations predicted to severely alter protein coding per animal, and identified 45 such variants that coincide with known pathogenic human variants. These data suggest that expanded screening of existing breeding colonies will identify novel models of human disease, and that increased genomic characterization can help inform research studies in macaques.
Subject(s)
Disease Models, Animal , Genetic Diseases, Inborn/genetics , Macaca mulatta/genetics , Mutation , Polymorphism, Genetic , Animals , Genomics , Sequence Analysis, DNAABSTRACT
BACKGROUND: Rhesus macaques are widely used in biomedical research, but the application of genomic information in this species to better understand human disease is still in its infancy. Whole-genome sequence (WGS) data in large pedigreed macaque colonies could provide substantial experimental power for genetic discovery, but the collection of WGS data in large cohorts remains a formidable expense. Here, we describe a cost-effective approach that selects the most informative macaques in a pedigree for 30X WGS, followed by low-cost genotyping-by-sequencing (GBS) at 30X on the remaining macaques in order to generate sparse genotype data at high accuracy. Dense variants from the selected macaques with WGS data are then imputed into macaques having only sparse GBS data, resulting in dense genome-wide genotypes throughout the pedigree. RESULTS: We developed GBS for the macaque genome using a digestion with PstI, followed by sequencing of size-selected fragments at 30X coverage. From GBS sequence data collected on all individuals in a 16-member pedigree, we characterized high-confidence genotypes at 22,455 single nucleotide variant (SNV) sites that were suitable for guiding imputation of dense sequence data from WGS. To characterize dense markers for imputation, we performed WGS at 30X coverage on nine of the 16 individuals, yielding 10,193,425 high-confidence SNVs. To validate the use of GBS data for facilitating imputation, we initially focused on chromosome 19 as a test case, using an optimized panel of 833 sparse, evenly-spaced markers from GBS and 5,010 dense markers from WGS. Using the method of "Genotype Imputation Given Inheritance" (GIGI), we evaluated the effects on imputation accuracy of 3 different strategies for selecting individuals for WGS, including 1) using "GIGI-Pick" to select the most informative individuals, 2) using the most recent generation, or 3) using founders only. We also evaluated the effects on imputation accuracy of using a range of from 1 to 9 WGS individuals for imputation. We found that the GIGI-Pick algorithm for selection of WGS individuals outperformed common heuristic approaches, and that genotype numbers and accuracy improved very little when using >5 WGS individuals for imputation. Informed by our findings, we used 4 macaques with WGS data to impute variants at up to 7,655,491 sites spanning all 20 autosomes in the 12 remaining macaques, based on their GBS genotypes at only 17,158 loci. Using a strict confidence threshold, we imputed an average of 3,680,238 variants per individual at >99 % accuracy, or an average 4,458,883 variants per individual at a more relaxed threshold, yielding >97 % accuracy. CONCLUSIONS: We conclude that an optimal tradeoff between genotype accuracy, number of imputed genotypes, and overall cost exists at the ratio of one individual selected for WGS using the GIGI-Pick algorithm, per 3-5 relatives selected for GBS. This approach makes feasible the collection of accurate, dense genome-wide sequence data in large pedigreed macaque cohorts without the need for more expensive WGS data on all individuals.
Subject(s)
Genotyping Techniques/methods , Macaca mulatta/genetics , Sequence Analysis, DNA/methods , Algorithms , Animals , Chromosomes/genetics , Computational Biology/economics , Computational Biology/methods , Genotyping Techniques/economics , Polymorphism, Single Nucleotide , Sequence Analysis, DNA/economicsABSTRACT
Compensatory mutations offset fitness defects resulting from CD8(+) T lymphocyte (CD8(TL))-mediated escape, but their impact on viral evolution following transmission to naive hosts remains unclear. Here, we investigated the reversion kinetics of Gag(181-189)CM9 CD8(TL) escape-associated compensatory mutations in simian immunodeficiency virus (SIV)-infected macaques. Preexisting compensatory mutations did not result in acute-phase escape of the SIVmac239 CD8(TL) epitope Gag(181-189)CM9 and instead required a tertiary mutation for stabilization in the absence of Gag(181-189)CM9 escape mutations. Therefore, transmitted compensatory mutations do not necessarily predict rapid CD8(TL) escape.
Subject(s)
CD8-Positive T-Lymphocytes/immunology , Gene Products, gag/genetics , Gene Products, gag/immunology , Mutation , Simian Acquired Immunodeficiency Syndrome/immunology , Simian Immunodeficiency Virus/immunology , Amino Acid Sequence , Animals , CD8-Positive T-Lymphocytes/virology , Epitopes, T-Lymphocyte/genetics , Epitopes, T-Lymphocyte/immunology , Macaca mulatta , Molecular Sequence Data , Simian Acquired Immunodeficiency Syndrome/transmission , Simian Acquired Immunodeficiency Syndrome/virology , Simian Immunodeficiency Virus/geneticsABSTRACT
Simian immunodeficiency virus (SIV) vaccines based upon 68-1 Rhesus Cytomegalovirus (RhCMV) vectors show remarkable protection against pathogenic SIVmac239 challenge. Across multiple independent rhesus macaque (RM) challenge studies, nearly 60% of vaccinated RM show early, complete arrest of SIVmac239 replication after effective challenge, whereas the remainder show progressive infection similar to controls. Here, we performed viral sequencing to determine whether the failure to control viral replication in non-protected RMs is associated with the acquisition of viral escape mutations. While low level viral mutations accumulated in all animals by 28 days-post-challenge, which is after the establishment of viral control in protected animals, the dominant circulating virus in virtually all unprotected RMs was nearly identical to the challenge stock, and there was no difference in mutation patterns between this cohort and unvaccinated controls. These data definitively demonstrate that viral mutation does not explain lack of viral control in RMs not protected by RhCMV/SIV vaccination. We further demonstrate that during chronic infection RhCMV/SIV vaccinated RMs do not acquire escape mutation in epitopes targeted by RhCMV/SIV, but instead display mutation in canonical MHC-Ia epitopes similar to unvaccinated RMs. This suggests that after the initial failure of viral control, unconventional T cell responses induced by 68-1 RhCMV/SIV vaccination do not exert strong selective pressure on systemically replicating SIV.
Subject(s)
Macaca mulatta , Mutation , SAIDS Vaccines , Simian Acquired Immunodeficiency Syndrome , Simian Immunodeficiency Virus , Animals , Simian Immunodeficiency Virus/immunology , Simian Immunodeficiency Virus/genetics , Simian Acquired Immunodeficiency Syndrome/immunology , Simian Acquired Immunodeficiency Syndrome/virology , Simian Acquired Immunodeficiency Syndrome/prevention & control , SAIDS Vaccines/immunology , SAIDS Vaccines/genetics , Cytomegalovirus/immunology , Cytomegalovirus/genetics , Virus Replication/immunology , Vaccination , Immune Evasion/geneticsABSTRACT
Molecular interactions between killer immunoglobulin-like receptors (KIRs) and their MHC class I ligands play a central role in the regulation of natural killer (NK) cell responses to viral pathogens and tumors. Here we identify Mamu-A1*00201 (Mamu-A*02), a common MHC class I molecule in the rhesus macaque with a canonical Bw6 motif, as a ligand for Mamu-KIR3DL05. Mamu-A1*00201 tetramers folded with certain SIV peptides, but not others, directly stained primary NK cells and Jurkat cells expressing multiple allotypes of Mamu-KIR3DL05. Differences in binding avidity were associated with polymorphisms in the D0 and D1 domains of Mamu-KIR3DL05, whereas differences in peptide-selectivity mapped to the D1 domain. The reciprocal exchange of the third predicted MHC class I-contact loop of the D1 domain switched the specificity of two Mamu-KIR3DL05 allotypes for different Mamu-A1*00201-peptide complexes. Consistent with the function of an inhibitory KIR, incubation of lymphocytes from Mamu-KIR3DL05(+) macaques with target cells expressing Mamu-A1*00201 suppressed the degranulation of tetramer-positive NK cells. These observations reveal a previously unappreciated role for D1 polymorphisms in determining the selectivity of KIRs for MHC class I-bound peptides, and identify the first functional KIR-MHC class I interaction in the rhesus macaque. The modulation of KIR-MHC class I interactions by viral peptides has important implications to pathogenesis, since it suggests that the immunodeficiency viruses, and potentially other types of viruses and tumors, may acquire changes in epitopes that increase the affinity of certain MHC class I ligands for inhibitory KIRs to prevent the activation of specific NK cell subsets.
Subject(s)
HLA-B Antigens/genetics , Peptides/metabolism , Polymorphism, Single Nucleotide , Potassium Channels, Inwardly Rectifying/genetics , Animals , Cloning, Molecular , Epitopes, T-Lymphocyte/metabolism , Histocompatibility Antigens Class I/metabolism , Host-Pathogen Interactions , Humans , Jurkat Cells , Macaca mulatta , Potassium Channels, Inwardly Rectifying/metabolism , Protein Binding , Protein Conformation , Simian Acquired Immunodeficiency Syndrome/genetics , Simian Acquired Immunodeficiency Syndrome/metabolism , Simian Acquired Immunodeficiency Syndrome/virology , T-Lymphocytes/metabolism , TransfectionABSTRACT
We present JBrowse 2, a general-purpose genome annotation browser offering enhanced visualization of complex structural variation and evolutionary relationships. It retains core features of JBrowse while adding new views for synteny, dotplots, breakpoints, gene fusions, and whole-genome overviews. It allows users to share sessions, open multiple genomes, and navigate between views. It can be embedded in a web page, used as a standalone application, or run from Jupyter notebooks or R sessions. These improvements are enabled by a ground-up redesign using modern web technology. We describe application functionality, use cases, performance benchmarks, and implementation notes for web administrators and developers.