In-depth examination of the plasma proteomic response to infection with a wide variety of pathogens can assist in the development of new diagnostic paradigms, while providing insight into the interdependent pathogenic processes which encompass a host’s immunological and physiological responses. Ebola virus (EBOV) causes a highly lethal infection termed Ebola virus disease (EVD) in primates and humans. The Gram negative non-spore forming bacillus Burkholderia pseudomallei (Bp) causes melioidosis in primates and humans, characterized by severe pneumonia with high mortality. We sought to examine the host response to infection with these two bio-threat pathogens using established animal models to provide information on the feasibility of pre-symptomatic diagnosis, since the induction of host molecular signaling networks can occur before clinical presentation and pathogen detection. Herein we report the quantitative proteomic analysis of plasma collected at various times of disease progression from 10 EBOV-infected and 5 Bp-infected nonhuman primates (NHP). Our strategy employed high resolution LC–MS/MS and a peptide-tagging approach for relative protein quantitation. In each infection type, for all proteins with > 1.3 fold abundance change at any post-infection time point, a direct comparison was made with levels obtained from plasma collected daily from 5 naïve rhesus macaques, to determine the fold changes that were significant, and establish the natural variability of abundance for endogenous plasma proteins. A total of 41 plasma proteins displayed significant alterations in abundance during EBOV infection, and 28 proteins had altered levels during Bp infection, when compared to naïve NHPs. Many major acute phase proteins quantitated displayed similar fold-changes between the two infection types but exhibited different temporal dynamics. Proteins related to the clotting cascade, immune signaling and complement system exhibited significant differential abundance during infection with EBOV or Bp, indicating a specificity of the response. These results advance our understanding of the global plasma proteomic response to EBOV and Bp infection in relevant primate models for human disease and provide insight into potential innate immune response differences between viral and bacterial infections.

Whole-genome sequencing and comparative genomic analysis of immune-related genes of Cercocebus atys and Macaca mulatta identify candidate genes, such as ICAM2 and TLR4, that may explain the AIDS resistance of C. atys. Sooty mangabeys are a primate model for non-pathogenic simian immunodeficiency virus (SIV) infection. Guido Silvestri and colleagues sequenced the genome of a captive sooty mangabey from the Yerkes National Primate Research Center in Atlanta, Georgia. They compared this with the genomic sequences of other primates and suggest several candidate genes that may influence differences in SIV pathogenicity across nonhuman primates. The team also analysed genetic variation in the genomes of 10 sooty mangabeys from the same centre and the transcriptomes of 14 tissues. This work provides a resource for comparative genomics and studies on SIV pathogenicity. In contrast to infections with human immunodeficiency virus (HIV) in humans and simian immunodeficiency virus (SIV) in macaques, SIV infection of a natural host, sooty mangabeys (Cercocebus atys), is non-pathogenic despite high viraemia1. Here we sequenced and assembled the genome of a captive sooty mangabey. We conducted genome-wide comparative analyses of transcript assemblies from C. atys and AIDS-susceptible species, such as humans and macaques, to identify candidates for host genetic factors that influence susceptibility. We identified several immune-related genes in the genome of C. atys that show substantial sequence divergence from macaques or humans. One of these sequence divergences, a C-terminal frameshift in the toll-like receptor-4 (TLR4) gene of C. atys, is associated with a blunted in vitro response to TLR-4 ligands. In addition, we found a major structural change in exons 3–4 of the immune-regulatory protein intercellular adhesion molecule 2 (ICAM-2); expression of this variant leads to reduced cell surface expression of ICAM-2. These data provide a resource for comparative genomic studies of HIV and/or SIV pathogenesis and may help to elucidate the mechanisms by which SIV-infected sooty mangabeys avoid AIDS.

Due to the large geographical overlap of populations exposed to Zika virus (ZIKV) and human immunodeficiency virus (HIV), understanding the disease pathogenesis of co-infection is urgently needed. This warrants the development of an animal model for HIV-ZIKV co-infection. In this study, we used adult non-pregnant macaques that were chronically infected with simian immunodeficiency virus/chimeric simian human immunodeficiency virus (SIV/SHIV) and then inoculated with ZIKV. Plasma viral loads of both SIV/SHIV and ZIKV co-infected animals revealed no significant changes as compared to animals that were infected with ZIKV alone or as compared to SIV/SHIV infected animals prior to ZIKV inoculation. ZIKV tissue clearance of co-infected animals was similar to animals that were infected with ZIKV alone. Furthermore, in co-infected macaques, there was no statistically significant difference in plasma cytokines/chemokines levels as compared to prior to ZIKV inoculation. Collectively, these findings suggest that co-infection may not alter disease pathogenesis, thus warranting larger HIV-ZIKV epidemiological studies in order to validate these findings.

Zika Virus (ZIKV) is primarily transmitted through mosquito bites. It can also be transmitted during sexual intercourse and in utero from mother to fetus. To gain preliminary insight into ZIKV pathology and immune responses on route of transmission, rhesus macaques (RMs) were inoculated with ZIKV (PRVABC59) via intravaginal (IVAG) (n = 3) or subcutaneous (sub Q) (n = 2) routes. Systemic ZIKV infection was observed in all RMs, regardless of the route of inoculation. After 9 days postinfection (dpi), ZIKV was not detected in the plasma of IVAG- and sub-Q-inoculated RMs. Importantly, RMs harbored ZIKV up to 60 dpi in various anatomical locations. Of note, ZIKV was also present in several regions of the brain, including the caudate nucleus, parietal lobe, cortex, and amygdala. These observations appear to indicate that ZIKV infection may be systemic and persistent regardless of route of inoculation. In addition, we observed changes in key immune cell populations in response to ZIKV infection. Importantly, IVAG ZIKV infection of RMs is associated with increased depletion of CD11C hi myeloid cells, reduced PD-1 expression in NK cells, and elevated frequencies of Ki67(+) CD8(+) central memory cells as compared to sub Q ZIKV-infected RMs. These results need to interpreted with caution due to the small number of animals utilized in this study. Future studies involving large groups of animals that have been inoculated through both routes of transmission are needed to confirm our findings.

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.

Background Exome sequencing has been utilized to identify genetic variants associated with disease in humans. Identification of loss-of-function mutations with exome sequencing in rhesus macaques ( Macaca mulatta ) could lead to valuable animal models of genetic disease. Attempts have been made to identify variants in rhesus macaques by aligning exome data against the rheMac2 draft genome. However, such efforts have been impaired due to the incompleteness and annotation errors associated with rheMac2. We wished to determine whether aligning exome reads against our new, improved rhesus genome, MacaM, could be used to identify high impact, loss-of-function mutations in rhesus macaques that would be relevant to human disease. Results We compared alignments of exome reads from four rhesus macaques, the reference animal and three unrelated animals, against rheMac2 and MacaM. Substantially more reads aligned against MacaM than rheMac2. We followed the Broad Institute’s Best Practice guidelines for variant discovery which utilizes the Genome Analysis Toolkit to identify high impact mutations. When rheMac2 was used as the reference genome, a large number of apparent false positives were identified. When MacaM was used as the reference genome, the number of false positives was greatly reduced. After examining the variant analyses conducted with MacaM as reference genome, we identified two putative loss-of-function mutations, in the heterozygous state, in genes related to human health. Sanger sequencing confirmed the presence of these mutations. We followed the transmission of one of these mutations (in the butyrylthiocholine gene) through three generations of rhesus macaques. Further, we demonstrated a functional decrease in butyrylthiocholinesterase activity similar to that observed in human heterozygotes with loss-of-function mutations in the same gene. Conclusions The new MacaM genome can be effectively utilized to identify loss-of-function mutations in rhesus macaques without generating a high level of false positives. In some cases, heterozygotes may be immediately useful as models of human disease. For diseases where homozygous mutants are needed, directed breeding of loss-of-function heterozygous animals could be used to create rhesus macaque models of human genetic disease. The approach we describe here could be applied to other mammals, but only if their genomes have been improved beyond draft status.

Acetylcholinesterase (AChE) is the physiologically important target for organophosphorus toxicants (OP) including nerve agents and pesticides. Butyrylcholinesterase (BChE) in blood serves as a bioscavenger that protects AChE in nerve synapses from inhibition by OP. Mass spectrometry methods can detect exposure to OP by measuring adducts on the active site serine of plasma BChE. Genetic variants of human AChE and BChE do exist, but loss of function mutations have been identified only in the BCHE gene. The most common AChE variant, His353Asn (H322N), also known as the Yt blood group antigen, has normal AChE activity. The most common BChE variant, Ala567Thr (A539T) or the K-variant in honor of Werner Kalow, has 33% reduced plasma BChE activity. The genetic variant most frequently associated with prolonged response to muscle relaxants, Asp98Gly (D70G) or atypical BChE, has reduced activity and reduced enzyme concentration. Early studies in young, healthy males, performed at a time when it was legal to test nerve agents in humans, showed that individuals responded differently to the same low dose of sarin with toxic symptoms ranging in severity from minimal to moderate. Additionally, animal studies indicated that BChE protects from toxicants that have a higher reactivity with AChE than with BChE (e.g., nerve agents) but not from toxicants that have a higher reactivity with BChE than with AChE (e.g., OP pesticides). As a corollary, we hypothesize that individuals with genetic variants of BChE may be at increased risk of toxicity from nerve agents but not from OP pesticides.

The rhesus macaque (Macaca mulatta) is widely used in molecular evolutionary analyses, particularly to identify genes under adaptive or unique evolution in the human lineage. For such studies, it is necessary to align nucleotide sequences of homologous protein-coding genes among multiple species. The validity of these analyses is dependent on high quality genomic data. However, for most mammalian species (other than humans and mice), only draft genomes are available. There has been concern that some results obtained from evolutionary analyses using draft genomes may not be correct. The rhesus macaque provides a unique opportunity to determine whether an improved genome (MacaM) yields better results than a draft genome (rheMac2) for evolutionary studies. We compared protein-coding genes annotated in the rheMac2 and MacaM genomes with their human orthologs. We found many genes annotated in rheMac2 had apparently spurious sequences not present in genes derived from MacaM. The rheMac2 annotations also appeared to inflate a frequently used evolutionary index, ω (the ratio of nonsynonymous to synonymous substitution rates). Genes with these spurious sequences must be filtered out from evolutionary analyses to obtain correct results. With the MacaM genome, improved sequence information means many more genes can be examined for indications of selection. These results indicate how upgrading genomes from draft status to a higher level of quality can improve interpretation of evolutionary patterns.

Background Common chimpanzees ( Pan troglodytes ) and bonobos ( Pan paniscus ) are the species most closely related to humans. For this reason, it is especially important to have complete and accurate chimpanzee nucleotide and protein sequences to understand how humans evolved their unique capabilities. We provide transcriptome data from four untransformed cell types derived from the reference Pan troglodytes , “Clint”, to better annotate the chimpanzee genome and provide empirical validation for proposed gene models of this important species. Findings RNA was extracted from primary cells cultured from four tissues: skin, adipose stroma, vascular smooth muscle and skeletal muscle. These four RNA samples were sequenced on the Illumina HiSeq 2000 platform. Sequences were deposited in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA). Transcripts were assembled, annotated and deposited in the NCBI Transcriptome Shotgun Assembly (TSA) database. Conclusions We have provided a high quality annotation of 44,275 transcripts with full-length coding sequence (CDS). This set represented a total of 10,110 unique genes, thus providing empirical support for their existence. This dataset can be used to improve the annotation of the Pan troglodytes genome.

The timing of type I interferon signalling determines the disease course of SIV infection. Type I interferon (IFN-I) is shown here to have dual effects in rhesus macaques exposed to simian immunodeficiency virus (SIV): it is beneficial at the onset of infection but as infection progresses it becomes detrimental. IFN signaling was manipulated in two ways. IFN-I receptor blockade results in increased plasma viraemia, accelerated CD4 T cell loss and progression to AIDS. In contrast, IFN-α2a administration prior to high-dose intrarectal SIV challenge increases resistance to systemic infection. However, continued IFN-α2a treatment induces IFN-I desensitization and facilitates SIV infection. Overall, the benefits of early antiviral activity appear to outweigh the detrimental effects of immune activation during acute SIV infection. Inflammation in HIV infection is predictive of non-AIDS morbidity and death1, higher set point plasma virus load2 and virus acquisition3; thus, therapeutic agents are in development to reduce its causes and consequences. However, inflammation may simultaneously confer both detrimental and beneficial effects. This dichotomy is particularly applicable to type I interferons (IFN-I) which, while contributing to innate control of infection4,5,6,7,8,9,10, also provide target cells for the virus during acute infection, impair CD4 T-cell recovery, and are associated with disease progression6,7,11,12,13,14,15,16,17,18,19. Here we manipulated IFN-I signalling in rhesus macaques (Macaca mulatta) during simian immunodeficiency virus (SIV) transmission and acute infection with two complementary in vivo interventions. We show that blockade of the IFN-I receptor caused reduced antiviral gene expression, increased SIV reservoir size and accelerated CD4 T-cell depletion with progression to AIDS despite decreased T-cell activation. In contrast, IFN-α2a administration initially upregulated expression of antiviral genes and prevented systemic infection. However, continued IFN-α2a treatment induced IFN-I desensitization and decreased antiviral gene expression, enabling infection with increased SIV reservoir size and accelerated CD4 T-cell loss. Thus, the timing of IFN-induced innate responses in acute SIV infection profoundly affects overall disease course and outweighs the detrimental consequences of increased immune activation. Yet, the clinical consequences of manipulation of IFN signalling are difficult to predict in vivo and therapeutic interventions in human studies should be approached with caution.

Background: The rhesus macaque (Macaca mulatta) is a key species for advancing biomedical research. Like all draft mammalian genomes, the draft rhesus assembly (rheMac2) has gaps, sequencing errors and misassemblies that have prevented automated annotation pipelines from functioning correctly. Another rhesus macaque assembly, CR_1.0, is also available but is substantially more fragmented than rheMac2 with smaller contigs and scaffolds. Annotations for these two assemblies are limited in completeness and accuracy. High quality assembly and annotation files are required for a wide range of studies including expression, genetic and evolutionary analyses.Results: We report a new de novo assembly of the rhesus macaque genome (MacaM) that incorporates both the original Sanger sequences used to assemble rheMac2 and new Illumina sequences from the same animal. MacaM has a weighted average (N50) contig size of 64 kilobases, more than twice the size of the rheMac2 assembly and almost five times the size of the CR_1.0 assembly. The MacaM chromosome assembly incorporates information from previously unutilized mapping data and preliminary annotation of scaffolds. Independent assessment of the assemblies using Ion Torrent read alignments indicates that MacaM is more complete and accurate than rheMac2 and CR_1.0. We assembled messenger RNA sequences from several rhesus tissues into transcripts which allowed us to identify a total of 11,712 complete proteins representing 9,524 distinct genes. Using a combination of our assembled rhesus macaque transcripts and human transcripts, we annotated 18,757 transcripts and 16,050 genes with complete coding sequences in the MacaM assembly. Further, we demonstrate that the new annotations provide greatly improved accuracy as compared to the current annotations of rheMac2. Finally, we show that the MacaM genome provides an accurate resource for alignment of reads produced by RNA sequence expression studies.Conclusions: The MacaM assembly and annotation files provide a substantially more complete and accurate representation of the rhesus macaque genome than rheMac2 or CR_1.0 and will serve as an important resource for investigators conducting next-generation sequencing studies with nonhuman primates.Reviewers: This article was reviewed by Dr. Lutz Walter, Dr. Soojin Yi and Dr. Kateryna Makova.

The study of nonhuman primates (NHP) is key to understanding human evolution, in addition to being an important model for biomedical research. NHPs are especially important for translational medicine. There are now exciting opportunities to greatly increase the utility of these models by incorporating Next Generation (NextGen) sequencing into study design. Unfortunately, the draft status of nonhuman genomes greatly constrains what can currently be accomplished with available technology. Although all genomes contain errors, draft assemblies and annotations contain so many mistakes that they make currently available nonhuman primate genomes misleading to investigators conducting evolutionary studies; and these genomes are of insufficient quality to serve as references for NextGen studies. Fortunately, NextGen sequencing can be used in the production of greatly improved genomes. Existing Sanger sequences can be supplemented with NextGen whole genome, and exomic genomic sequences to create new, more complete and correct assemblies. Additional physical mapping, and an incorporation of information about gene structure, can be used to improve assignment of scaffolds to chromosomes. In addition, mRNA-sequence data can be used to economically acquire transcriptome information, which can be used for annotation. Some highly polymorphic and complex regions, for example MHC class I and immunoglobulin loci, will require extra effort to properly assemble and annotate. However, for the vast majority of genes, a modest investment in money, and a somewhat greater investment in time, can greatly improve assemblies and annotations sufficient to produce true, reference grade nonhuman primate genomes. Such resources can reasonably be expected to transform nonhuman primate research.

Finished genome sequences and assemblies are available for only a few vertebrates. Thus, investigators studying many species must rely on draft genomes. Using the rhesus macaque as an example, we document the effects of sequencing errors, gaps in sequence and misassemblies on one automated gene model pipeline, Gnomon. The combination of draft genome with automated gene finding software can result in spurious sequences. We estimate that approximately 50% of the rhesus gene models are missing, incomplete or incorrect. The problems identified in this work likely apply to all draft vertebrate genomes annotated with any automated gene model pipeline and thus represent a pervasive challenge to the analysis of draft genomes.

BACKGROUND:Affymetrix GeneChips utilize 25-mer oligonucleotides probes linked to a silica surface to detect targets in solution. Mismatches due to single nucleotide polymorphisms (SNPs) can affect the hybridization between probes and targets. Previous research has indicated that binding between probes and targets strongly depends on the positions of these mismatches. However, there has been substantial variability in the effect of mismatch type across studies.METHODS:By taking advantage of naturally occurring mismatches between rhesus macaque transcripts and human probes from the Affymetrix U133 Plus 2 GeneChip, we collected the largest 25-mer probes dataset with single-base mismatches at each of the 25 positions on the probe ever used in this type of analysis.RESULTS:A mismatch at the center of a probe led to a greater loss in signal intensity than a mismatch at the ends of the probe, regardless of the mismatch type. There was a slight asymmetry between the ends of a probe: effects of mismatches at the 3' end of a probe were greater than those at the 5' end. A cross study comparison of the effect of mismatch types revealed that results were not in good agreement among different reports. However, if the mismatch types were consolidated to purine or pyrimidine mismatches, cross study conclusions could be generated.CONCLUSION:The comprehensive assessment of the effects of single-base mismatches on microarrays provided in this report can be useful for improving future versions of microarray platform design and the corresponding data analysis algorithms.

Background: Gene targeting in nonhuman primates has the potential to produce critical animal models for translational studies related to human diseases. Successful gene targeting in fibroblasts followed by somatic cell nuclear transfer (SCNT) has been achieved in several species of large mammals but not yet in primates. Our goal was to establish the protocols necessary to achieve gene targeting in primary culture of adult rhesus macaque fibroblasts as a first step in creating nonhuman primate models of genetic disease using nuclear transfer technology.Results: A primary culture of adult male fibroblasts was transfected with hTERT to overcome senescence and allow long term in vitro manipulations. Successful gene targeting of the HPRT locus in rhesus macaques was achieved by electroporating S-phase synchronized cells with a construct containing a SV40 enhancer.Conclusion: The cell lines reported here could be used for the production of null mutant rhesus macaque models of human genetic disease using SCNT technology. In addition, given the close evolutionary relationship and biological similarity between rhesus macaques and humans, the protocols described here may prove useful in the genetic engineering of human somatic cells.

BACKGROUND:Somatic cell nuclear transfer (SCNT) requires cytoplast-mediated reprogramming of the donor nucleus. Cytoplast factors such as maturation promoting factor are implicated based on their involvement in nuclear envelope breakdown (NEBD) and premature chromosome condensation (PCC). Given prior difficulties in SCNT in primates using conventional protocols, we hypothesized that the ability of cytoplasts to induce nuclear remodeling was instrumental in efficient reprogramming. METHODS:NEBD and PCC in monkey (Macaca mulatta) SCNT embryos were monitored by lamin A/C immunolabeling. RESULTS:Initially, a persistent lamin A/C signal from donor cell nuclei after fusion with cytoplasts was observed indicative of incomplete NEBD following SCNT and predictive of developmental arrest. We then identified fluorochrome-assisted enucleation and donor cell electrofusion as likely candidates for inducing premature cytoplast activation and a consequent lack of nuclear remodeling. Modified protocols designed to prevent premature cytoplast activation during SCNT showed robust NEBD and PCC. Coincidently, over 20% of SCNT embryos reconstructed with fetal fibroblasts progressed to blastocysts. Similar results were obtained with other somatic cells. Reconstructed blastocysts displayed patterns of Oct-4 expression similar to fertilized embryos reflecting successful reprogramming. CONCLUSIONS:Our results represent a significant breakthrough in elucidating the role of nuclear remodeling events in reprogramming following SCNT.

BACKGROUND:Rhesus macaques serve a critical role in the study of human biomedical research. While both Indian and Chinese rhesus macaques are commonly used, genetic differences between these two subspecies affect aspects of their behavior and physiology, including response to simian immunodeficiency virus (SIV) infection. Single nucleotide polymorphisms (SNPs) can play an important role in both establishing ancestry and in identifying genes involved in complex diseases. We sequenced the 3' end of rhesus macaque genes in an effort to identify gene-based SNPs that could distinguish between Indian and Chinese rhesus macaques and aid in association analysis.RESULTS:We surveyed the 3' end of 94 genes in 20 rhesus macaque animals. The study included 10 animals each of Indian and Chinese ancestry. We identified a total of 661 SNPs, 457 of which appeared exclusively in one or the other population. Seventy-nine additional animals were genotyped at 44 of the population-exclusive SNPs. Of those, 38 SNPs were confirmed as being population-specific.CONCLUSION:This study demonstrates that the 3' end of genes is rich in sequence polymorphisms and is suitable for the efficient discovery of gene-linked SNPs. In addition, the results show that the genomic sequences of Indian and Chinese rhesus macaque are remarkably divergent, and include numerous population-specific SNPs. These ancestral SNPs could be used for the rapid scanning of rhesus macaques, both to establish animal ancestry and to identify gene alleles that may contribute to the phenotypic differences observed in these populations.

BACKGROUND:The non-human primate (NHP) research community has been intensely interested in obtaining whole-genome expression arrays for their work. Recently, novel approaches were used to generate the DNA sequence information for a rhesus GeneChip. To test the reliability of the rhesus GeneChip across different centers, RNA was isolated from five sources: cerebral cortex, pancreas, thymus, testis, and an immortalized fibroblast cell line. Aliquots of this RNA were sent to each of three centers: Yerkes National Primate Research Center, Oregon National Primate Research Center and the University of Nebraska Medical Center. Each center labeled the samples and hybridized them with two rhesus macaque GeneChips. In addition, rhesus samples were hybridzed with human GeneChips to compare with samples hybridized with the rhesus GeneChip.RESULTS:The results indicate that center effects were minimal and the rhesus GeneChip appears highly reliable. To test the validity of the rhesus GeneChip, five of the most differentially expressed genes among tissues identified in the reliability experiments were chosen for analysis with Quantitative PCR. For all 5 genes, the qPCR and GeneChip results were in agreement with regard to differential expression between tissues. Significantly more probesets were called present when rhesus samples were hybridized with the rhesus GeneChip than when these same samples were hybridized with a human GeneChip.CONCLUSION:The rhesus GeneChip is both a reliable and a valid tool for examining gene expression and represents a significant improvement over the use of the human GeneChip for rhesus macaque gene expression studies.

Non-human primates (NHPs) serve as excellent models for transplantation research. Gene expression profiling has been shown to be an important tool in understanding many aspects of transplantation. The many advantages of using expression arrays, and the unavailability of monkey specific arrays (until recently), has prompted some investigators to use human expression arrays with samples from rhesus macaques (Macaca mulatta). This review examines attempts to use human expression arrays with monkey samples and the recent development of monkey-specific expression arrays.