INTERSPECIES HYBRIDIZATION AND THE STRATIFICATION OF NUCLEAR GENETIC VARIATION OF RHESUS (MACACA MULATTA) AND LONG-TAILED MACAQUES (MACACA FASCICULARIS)

Prevalence and genetic diversity of simian malaria in wild macaque populations across Thailand: Implications for human health

Over the past year, P. falciparum infections have declined in Thailand, yet nonhuman primate malaria infections have correspondingly increased, including Plasmodium knowlesi and P. cynomolgi. Nevertheless, little is known about simian malaria in its natural macaque hosts, Macaca mulatta and Macaca fascicularis. This study aims to address several research questions, including the prevalence and distribution of simian malaria in these two Thai wild macaque species, variations in infection between different macaque species and between M. fascicularis subspecies, and the genetic composition of these pathogens. Blood samples were collected from 82 M. mulatta and 690 M. fascicularis across 15 locations in Thailand, as well as two locations in Vietnam and Myanmar. We employed quantitative real-time PCR targeting the Plasmodium genus-specific 18S ribosomal RNA (rRNA) gene to detect malaria infection, with a limit of detection set at 1,215.98 parasites per mL. We genotyped eight microsatellite markers, and the P. cynomolgi dihydrofolate reductase gene (DHFR) was sequenced (N = 29). In total, 100 of 772 samples (13 %) tested positive for malaria, including 45 (13 %) for P. cynomolgi, 37 (13 %) for P. inui, 16 (5 %) for P. coatneyi, and 2 (0.25 %) for Hepatocystis sp. in Saraburi, central and Ranong, southern Thailand. Notably, simian malaria infection was observed exclusively in M. fascicularis and not in M. mulatta (P = 0.0002). Particularly, P. cynomolgi was detected in 21.7 % (45/207) of M. f. fascicularis living in Wat Tham Phrapothisat, Saraburi Province. The infection with simian malaria was statistically different between M. fascicularis and M. mulatta (P = 0.0002) but not within M. fascicularis subspecies (P = 0.78). A haplotype network analysis revealed that P. cynomolgi shares a lineage with reference strains obtained from macaques. No mutation in the predicted binding pocket of PcyDHFR to pyrimethamine was observed. This study reveals a significant prevalence of simian malaria infection in M. fascicularis. The clonal genotypes of P. cynomolgi suggest in-reservoir breeding. These findings raise concerns about the potential spread of nonhuman primate malaria to humans and underscore the need for preventive measures.

Gene flow from rhesus (Macaca mulatta) to cynomolgus macaques (M. fascicularis) and effects of introgressive hybridization on reproduction in two biomedically relevant non-human primate species

BACKGROUND

We compared the reproductive patterns of wild Indochinese and Sundaic cynomolgus macaques (Mf) exhibiting different levels of genetic admixture with rhesus macaques (Mm).

METHODS

Ten adult females from each Indochinese (WHM) and Sundaic (KN/KTK) Mf populations, which exhibited 50% and 15% of Mm autosomal SNPs, were selected as focal animals. Animals were observed for 12 months, and the frequencies of sexual proceptivity, attractivity and receptivity, number of newborns, and changes in sex skin were recorded.

RESULTS

Both populations showed all three sexual behaviors throughout the year, but they were classified as moderately seasonal breeders because their 3-month birth counts were as high as ~50%. The fecundity of WHM was lower than the KN/KTK. Changes in sex skin of WHM were more prone to Mm's pattern than the KN/KTK.

CONCLUSION

The introgressive gene flow from Mm to Mf does not affect Mf's sexual behaviors; however, it can impact fecundity and physiological (sex skin) changes.

Population Structure of Macaca fascicularis aurea, and their Genetic Relationships with M. f. fascicularis and M. mulatta Determined by 868 RADseq-Derived Autosomal SNPs-A consideration for biomedical research

BACKGROUND

This study examined the population structure of Macaca fascicularis aurea and their genetic relationships with M. f. fascicularis and M. mulatta.

METHODS

The study analyzed 868 RADseq-derived SNPs from samples representing the entire distribution range of M. f. aurea, including their inter- and intraspecific hybrid zones.

RESULTS

The study supports a M. mulatta/Indochinese M. f. fascicularis, Sundaic M. f. fascicularis, and M. f. aurea trichotomy; M. f. aurea was genetically distinct from both forms of M. f. fascicularis and M. mulatta. Hybridization between M. f. aurea and M. f. fascicularis occurred in two directions: south-north (8°25' to 15°56') and west-east (98°28' to 99°02'). Low levels of M. mulatta introgression were also detected in M. f. aurea.

CONCLUSION

This study showcases a complicated scenario of genetic relationships between the M. fascicularis subspecies and between M. fascicularis and M. mulatta and underscores the importance of these taxa's population structure and genetic relationships for biomedical research.

Analysis of Simian Endogenous Retrovirus (SERV) Full-Length Proviruses in Old World Monkey Genomes

Simian endogenous retrovirus, SERV, is a successful germ line invader restricted to Old World monkey (OWM) species. (1) Background: The availability of high-quality primate genomes warrants a study of the characteristics, evolution, and distribution of SERV proviruses. (2) Methods: Cercopithecinae OWM genomes from public databases were queried for the presence of full-length SERV proviruses. A dataset of 81 Cer-SERV genomes was generated and analyzed. (3) Results: Full-length Cer-SERV proviruses were mainly found in terrestrial OWM, and less so in arboreal, forest- dwelling monkeys. Phylogenetic analysis confirmed the existence of two genotypes, Cer-SERV-1 and Cer-SERV-2, with Cer-SERV-1 showing evidence of recent germ-line expansions. Long Terminal Repeat (LTR) variation indicated that most proviruses were of a similar age and were estimated to be between <0.3 and 10 million years old. Integrations shared between species were relatively rare. Sequence analysis further showed extensive CpG methylation-associated mutations, variable Primer Binding Site (PBS) use with Cer-SERV-1 using PBSlys3 and Cer-SERV-2 using PBSlys1,2, and the recent gain of LTR motifs for transcription factors active during embryogenesis in Cer-SERV-1. (4) Conclusions: sequence analysis of 81 SERV proviruses from Cercopithecinae OWM genomes provides evidence for the adaptation of this retrovirus to germ line reproduction.

Dramatic transcriptomic differences in Macaca mulatta and Macaca fascicularis with Plasmodium knowlesi infections

Plasmodium knowlesi, a model malaria parasite, is responsible for a significant portion of zoonotic malaria cases in Southeast Asia and must be controlled to avoid disease severity and fatalities. However, little is known about the host-parasite interactions and molecular mechanisms in play during the course of P. knowlesi malaria infections, which also may be relevant across Plasmodium species. Here we contrast P. knowlesi sporozoite-initiated infections in Macaca mulatta and Macaca fascicularis using whole blood RNA-sequencing and transcriptomic analysis. These macaque hosts are evolutionarily close, yet malaria-naïve M. mulatta will succumb to blood-stage infection without treatment, whereas malaria-naïve M. fascicularis controls parasitemia without treatment. This comparative analysis reveals transcriptomic differences as early as the liver phase of infection, in the form of signaling pathways that are activated in M. fascicularis, but not M. mulatta. Additionally, while most immune responses are initially similar during the acute stage of the blood infection, significant differences arise subsequently. The observed differences point to prolonged inflammation and anti-inflammatory effects of IL10 in M. mulatta, while M. fascicularis undergoes a transcriptional makeover towards cell proliferation, consistent with its recovery. Together, these findings suggest that timely detection of P. knowlesi in M. fascicularis, coupled with control of inflammation while initiating the replenishment of key cell populations, helps contain the infection. Overall, this study points to specific genes and pathways that could be investigated as a basis for new drug targets that support recovery from acute malaria.

Evolutionary pattern of Macaca fascicularis in Southeast Asia inferred using Y-chromosomal gene

BACKGROUND

We analyzed a combined segment (2032-bp) of the sex-determining region and the testis-specific protein of the Y-chromosome (Y-DNA) gene to clarify the gene flow and phylogenetic relationships of the long-tailed macaques (Macaca fascicularis) in Southeast Asia. Phylogenetic relationships were constructed using the maximum likelihood, Bayesian inference, and the median-joining network from a total of 164 adult male M. fascicularis from 62 localities in Malaysia, including sequences from the other regions from previous studies.

RESULTS

Based on Y-DNA, we confirm the presence of two lineages of M. fascicularis: the Indochinese and Sundaic lineages. The Indochinese lineage is represented by M. fascicularis located northwards of the Surat Thani-Krabi depression region and is introgressed by the Macaca mulatta Y-DNA. The Sundaic lineage is free from such hybridization event, thus defined as the original carrier of the M. fascicularis Y-DNA. We further revealed that the Sundaic lineage differentiated into two forms: the insular and the continental forms. The insular form, which represents the ancestral form of M. fascicularis, consists of two haplotypes: a single homogenous haplotype occupying the island of Borneo, Philippines, and southern Sumatra; and the Javan haplotype. The more diverse continental form consists of 17 haplotypes in which a dominant haplotype was shared by individuals from southern Thai Peninsular (south of Surat Thani-Krabi depression), Peninsular Malaysia, and Sumatra. Uniquely, Sumatra contains both the continental and insular Y-DNA which can be explained by a secondary contact hypothesis.

CONCLUSIONS

Overall, the findings in this study are important: (1) to help authority particularly in Malaysia on the population management activities including translocation and culling of conflict M. fascicularis, (2) to identify the unknown origin of captive M. fascicularis used in biomedical research, and; (3) the separation between the continental and insular forms warrants for the treatment as separate management units.

Secondary contact and genomic admixture between rhesus and long-tailed macaques in the Indochina Peninsula

Understanding the process and consequences of hybridization is one of the major challenges in evolutionary biology. A growing body of literature has reported evidence of ancient hybridization events or natural hybrid zones in primates, including humans; however, we still have relatively limited knowledge about the pattern and history of admixture because there have been little studies that simultaneously achieved genome-scale analysis and a geographically wide sampling of wild populations. Our study applied double-digest restriction site-associated DNA sequencing to samples from the six localities in and around the provisional hybrid zone of rhesus and long-tailed macaques and evaluated population structure, phylogenetic relationships, demographic history, and geographic clines of morphology and allele frequencies. A latitudinal gradient of genetic components was observed, highlighting the transition from rhesus (north) to long-tailed macaque distribution (south) as well as the presence of one northern population of long-tailed macaques exhibiting unique genetic structure. Interspecific gene flow was estimated to have recently occurred after an isolation period, and the migration rate from rhesus to long-tailed macaques was slightly greater than in the opposite direction. Although some rhesus macaque-biased alleles have widely introgressed into long-tailed macaque populations, the inflection points of allele frequencies have been observed as concentrated around the traditionally recognized interspecific boundary where morphology discontinuously changed; this pattern was more pronounced in the X chromosome than in autosomes. Thus, due to geographic separation before secondary contact, reproductive isolation could have evolved, contributing to the maintenance of an interspecific boundary and species-specific morphological characteristics.

MHC class I allele diversity in cynomolgus macaques of Vietnamese origin

Cynomolgus macaques (Macaca fascicularis, Mafa) have been used as important experimental animal models for carrying out biomedical researches. The results of biomedical experiments strongly depend on the immunogenetic background of animals, especially on the diversity of major histocompatibility complex (MHC) alleles. However, there is much less information available on the polymorphism of MHC class I genes in cynomolgus macaques, than is currently available for humans. In this study, we have identified 40 Mafa-A and 60 Mafa-B exons 2 and 3 sequences from 30 unrelated cynomolgus macaques of Vietnamese origin. Among these alleles, 28 are novel. As for the remaining 72 known alleles, 15 alleles are shared with other cynomolgus macaque populations and 32 are identical to alleles previously reported in other macaque species. A potential recombination event was observed between Mafa-A1*091:02 and Mafa-A1*057:01. In addition, the Mafa-A1 genes were found to be more diverse than human HLA-A and the functional residues for peptide binding sites (PBS) or TCR binding sites (TBS) in Mafa-A1 have greater variability than that for non-PBS or non-TBS regions. Overall, this study provides important information on the diversity of Mafa-A and Mafa-B alleles from Vietnamese origin, which may help researchers to choose the most appropriate animals for their studies.

DNA-based Determination of Ancestry in Cynomolgus Macaques (Macaca fascicularis)

Interest in the genetic composition of cynomolgus macaques (Macaca fascicularis) has increased due to the rising demand for NHP models in human biomedical research. Significant genetic differences among regional populations of cynomolgus macaques can confound interpretations of research results because they do not solely reflect differences in experimental treatment effects. Therefore, the common origin of cynomolgus macaques used as research subjects should be verified by using region-specific genetic markers to minimize the influence of underlying genetic variation among animals selected as research subjects on phenotypes under study. We compared the effectiveness of 18 short tandem repeat (STR) markers with that of 83 single-nucleotide polymorphism (SNP) markers to differentiate the ancestry of cynomolgus macaques from 6 different populations (Cambodia, Sumatra, Mauritius, Singapore, and the islands of Luzon and Zamboanga in the Philippines). Genetic diversity indices such as allele numbers and expected heterozygosity based on SNP were lower and exhibited lower standard errors than those provided by STR, probably because, unlike STR, most SNP are biallelic and consequently exhibit maximal expected heterozygosity values of 0.50. However, the standard error of estimates of observed heterozygosity based on SNP was higher than that for STR, perhaps reflecting sampling errors. Only 27 SNP were required to match the resolving power of 17 STR to detect population structure, that is, 1.6 SNP:1 STR. Whereas STR only differentiated the Mauritian population from all other populations, SNP detected 4 genetically distinct groups (Cambodia, Singapore-Sumatra, Mauritius, and Zamboanga). SNP are poised to become as valuable as STR for understanding and detecting genetic structure among cynomolgus macaques. Although STR will remain an important tool for cynomolgus macaque population studies, SNP have the potential to become the mainstream marker type.

Malaria parasites of long-tailed macaques in Sarawak, Malaysian Borneo: a novel species and demographic and evolutionary histories

Background

Non-human primates have long been identified to harbour different species of Plasmodium. Long-tailed macaques (Macaca fascicularis), in particular, are reservoirs for P. knowlesi, P. inui, P. cynomolgi, P. coatneyi and P. fieldi. A previous study conducted in Sarawak, Malaysian Borneo, however revealed that long-tailed macaques could potentially harbour novel species of Plasmodium based on sequences of small subunit ribosomal RNA and circumsporozoite genes. To further validate this finding, the mitochondrial genome and the apicoplast caseinolytic protease M genes of Plasmodium spp. were sequenced from 43 long-tailed macaque blood samples.

Results

Apart from several named species of malaria parasites, long-tailed macaques were found to be potentially infected with novel species of Plasmodium, namely one we refer to as “P. inui-like.” This group of parasites bifurcated into two monophyletic clades indicating the presence of two distinct sub-populations. Further analyses, which relied on the assumption of strict co-phylogeny between hosts and parasites, estimated a population expansion event of between 150,000 to 250,000 years before present of one of these sub-populations that preceded that of the expansion of P. knowlesi. Furthermore, both sub-populations were found to have diverged from a common ancestor of P. inui approximately 1.5 million years ago. In addition, the phylogenetic analyses also demonstrated that long-tailed macaques are new hosts for P. simiovale.

Conclusions

Malaria infections of long-tailed macaques of Sarawak, Malaysian Borneo are complex and include a novel species of Plasmodium that is phylogenetically distinct from P. inui. These macaques are new natural hosts of P. simiovale, a species previously described only in toque monkeys (Macaca sinica) in Sri Lanka. The results suggest that ecological factors could affect the evolution of malaria parasites.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1170-9) contains supplementary material, which is available to authorized users.

The Cynomolgus Macaque Natural History Model of Pneumonic Tularemia for Predicting Clinical Efficacy Under the Animal Rule

Francisella tularensis is a highly infectious Gram-negative bacterium that is the etiologic agent of tularemia in animals and humans and a Tier 1 select agent. The natural incidence of pneumonic tularemia worldwide is very low; therefore, it is not feasible to conduct clinical efficacy testing of tularemia medical countermeasures (MCM) in human populations. Development and licensure of tularemia therapeutics and vaccines need to occur under the Food and Drug Administration's (FDA's) Animal Rule under which efficacy studies are conducted in well-characterized animal models that reflect the pathophysiology of human disease. The Tularemia Animal Model Qualification (AMQ) Working Group is seeking qualification of the cynomolgus macaque (Macaca fascicularis) model of pneumonic tularemia under Drug Development Tools Qualification Programs with the FDA based upon the results of studies described in this manuscript. Analysis of data on survival, average time to death, average time to fever onset, average interval between fever and death, and bacteremia; together with summaries of clinical signs, necropsy findings, and histopathology from the animals exposed to aerosolized F. tularensis Schu S4 in five natural history studies and one antibiotic efficacy study form the basis for the proposed cynomolgus macaque model. Results support the conclusion that signs of pneumonic tularemia in cynomolgus macaques exposed to 300-3,000 colony forming units (cfu) aerosolized F. tularensis Schu S4, under the conditions described herein, and human pneumonic tularemia cases are highly similar. Animal age, weight, and sex of animals challenged with 300-3,000 cfu Schu S4 did not impact fever onset in studies described herein. This study summarizes critical parameters and endpoints of a well-characterized cynomolgus macaque model of pneumonic tularemia and demonstrates this model is appropriate for qualification, and for testing efficacy of tularemia therapeutics under Animal Rule.

Spontaneous Findings in the Eyes of Cynomolgus Monkeys (Macaca fascicularis) of Mauritian Origin

Spontaneous findings noted in the eyes of Mauritian cynomolgus monkeys are described and descriptions are supplemented with illustrations. Findings observed after extensive histopathologic examinations (20 to 44 sections per eye) from 20 control, 17 treatment-naive stock monkeys, and 2 findings noted in drug-treated monkeys that were considered to be spontaneous are included. Also included are findings from 361 control monkeys of routine toxicity studies performed at our laboratories, for most of which a standard histopathological examination of 1 section per eye was conducted. Common observations in monkeys examined extensively and in historical controls were limited to lymphocytic or mononuclear cell infiltrations of the uvea and/or conjunctiva/sclera and, less commonly observed, melanocytoma of the ciliary body or iris. Findings noted only in monkeys examined extensively consisted of inflammation of the conjunctiva, ora serrata cysts, glial nodules, focal degeneration of the retina, cystoid degeneration of the central retina, ballooning degeneration of the ciliary epithelium, cyst of the ciliary body, and decreased pigmentation of the retinal pigment epithelium. Changes recorded only in historical controls included retinal atrophy and nuclear displacement in the retina. Lesions are discussed and compared with pertinent literature.

Genetic analysis of samples from wild populations opens new perspectives on hybridization between long-tailed (Macaca fascicularis) and rhesus macaques (Macaca mulatta)

In the past decade, many researchers have published papers about hybridization between long-tailed and rhesus macaques. These previous works have proposed unidirectional gene flow with the Isthmus of Kra as the zoogeographical barrier of hybridization. However, these reports analyzed specimens of unknown origin and/or did not include specimens from Thailand, the center of the proposed area of hybridization. Collected specimens of long-tailed and rhesus macaques representing all suspected hybridization areas were examined. Blood samples from four populations each of long-tailed and rhesus macaques inhabiting Thailand, Myanmar, and Laos were collected and analyzed with conspecific references from China (for rhesus macaques) and multiple countries from Sundaic regions (for long-tailed macaques). Ninety-six single nucleotide polymorphism (SNP) markers specifically designed to interrogate admixture and ancestry were used in genotyping. We found genetic admixture maximized at the hybrid zone (15-20°N), as well as admixture signals of varying strength in both directions outside of the hybrid zone. These findings show that the Isthmus of Kra is not a barrier to gene flow from rhesus to long-tailed populations. However, to precisely identify a southernmost barrier, if in fact a boundary rather than simple isolation by distance exists, the samples from peninsular Malaysia must be included in the analysis. Additionally, a long-tailed to rhesus gene flow boundary was found between northern Thailand and Myanmar. Our results suggest that selection of long-tailed and rhesus macaques, the two most commonly used non-human primates for biomedical research, should take into account not only the species identification but also the origin of and genetic admixture within and between the species.

Evidence for a primate origin of zoonotic Helicobacter suis colonizing domesticated pigs

Helicobacter suis is the second most prevalent Helicobacter species in the stomach of humans suffering from gastric disease. This bacterium mainly inhabits the stomach of domesticated pigs, in which it causes gastric disease, but it appears to be absent in wild boars. Interestingly, it also colonizes the stomach of asymptomatic rhesus and cynomolgus monkeys. The origin of modern human-, pig- or non-human primate-associated H. suis strains in these respective host populations was hitherto unknown. Here we show that H. suis in pigs possibly originates from non-human primates. Our data suggest that a host jump from macaques to pigs happened between 100 000 and 15 000 years ago and that pig domestication has had a significant impact on the spread of H. suis in the pig population, from where this pathogen occasionally infects humans. Thus, in contrast to our expectations, H. suis appears to have evolved in its main host in a completely different way than its close relative Helicobacter pylori in humans.

ABO blood group phenotype frequency estimation using molecular phenotyping in rhesus and cynomolgus macaques

A much larger sample (N = 2369) was used to evaluate a previously reported distribution of the A, AB and B blood group phenotypes in rhesus and cynomolgus macaques from six different regional populations. These samples, acquired from 15 different breeding and research facilities in the United States, were analyzed using a real-time quantitative polymerase chain reaction (qPCR) assay that targets single nucleotide polymorphisms (SNPs) responsible for the macaque A, B and AB phenotypes. The frequency distributions of blood group phenotypes of the two species differ significantly from each other and significant regional differentiation within the geographic ranges of each species was also observed. The B blood group phenotype was prevalent in rhesus macaques, especially those from India, while the frequencies of the A, B and AB phenotypes varied significantly among cynomolgus macaques from different geographic regions. The Mauritian cynomolgus macaques, despite having originated in Indonesia, showed significant (P ≪ .01) divergence from the Indonesian animals at the ABO blood group locus. Most Mauritian animals belonged to the B blood group while the Indonesian animals were mostly A. The close similarity in blood group frequency distributions between the Chinese rhesus and Indochinese cynomolgus macaques demonstrates that the introgression between these two species extends beyond the zone of intergradation in Indochina. This study underscores the importance of ABO blood group phenotyping of the domestic supply of macaques and their biospecimens.

Tracing the phylogeographic history of Southeast Asian long-tailed macaques through mitogenomes of museum specimens

The biogeographical history of Southeast Asia is complicated due to the continuous emergences and disappearances of land bridges throughout the Pleistocene. Here, we use long-tailed macaques (Macaca fascicularis), which are widely distributed throughout the mainland and islands of Southeast Asia, asa model for better understanding the biogeographical patterns of diversification in this geographically complex region. A reliable intraspecific phylogeny including individuals from localities on oceanic islands, continental islands, and the mainland is needed to trace relatedness along with the pattern and timing of colonization in this region. We used high-throughput sequencing techniques to sequence mitochondrial genomes (mitogenomes) from 95 Southeast Asian M. fascicularis specimens housed at natural history museums around the world. To achieve a comprehensive picture, we more than tripled the mitogenome sample size for M. fascicularis from previous studies, and for the first time included documented samples from the Philippines and several small Indonesian islands. Confirming the result from a previous, recent intraspecific phylogeny for M. fascicularis, the newly reconstructed phylogeny of 135 specimens divides the samples into two major clades: Clade A includes haplotypes from the mainland and some from northern Sumatra, while Clade B includes all insular haplotypes along with lineages from southern Sumatra. This study resolves a previous disparity by revealing a disjunction in the origin of Sumatran macaques, with separate lineages originating within the two major clades, suggesting that at least two major migrations to Sumatra occurred. However, our dated phylogeny reveals that the two major clades split ∼1.88Ma, which is earlier than in previously published phylogenies. Our new data reveal that most Philippine macaque lineages diverged from the Borneo stock within the last ∼0.06-0.43Ma. Finally, our study provides insight into successful sequencing of DNA across museums and shotgun sequencing of DNA specimens asa method to sequence the mitogenome.

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.

Ancestry, Plasmodium cynomolgi prevalence and rhesus macaque admixture in cynomolgus macaques (Macaca fascicularis) bred for export in Chinese breeding farms

BACKGROUND

Most cynomolgus macaques (Macaca fascicularis) used in the United States as animal models are imported from Chinese breeding farms without documented ancestry. Cynomolgus macaques with varying rhesus macaque ancestry proportions may exhibit differences, such as susceptibility to malaria, that affect their suitability as a research model.

METHODS

DNA of 400 cynomolgus macaques from 10 Chinese breeding farms was genotyped to characterize their regional origin and rhesus ancestry proportion. A nested PCR assay was used to detect Plasmodium cynomolgi infection in sampled individuals.

RESULTS

All populations exhibited high levels of genetic heterogeneity and low levels of inbreeding and genetic subdivision. Almost all individuals exhibited an Indochinese origin and a rhesus ancestry proportion of 5%-48%. The incidence of P. cynomolgi infection in cynomolgus macaques is strongly associated with proportion of rhesus ancestry.

CONCLUSIONS

The varying amount of rhesus ancestry in cynomolgus macaques underscores the importance of monitoring their genetic similarity in malaria research.

Genome-wide mining and comparative analysis of microsatellites in three macaque species

Microsatellites are found in taxonomically different organisms, and such repeats are related with genomic structure, function and certain diseases. To characterize microsatellites for macaques, we searched and compared SSRs with 1-6 bp nucleotide motifs in rhesus, cynomolgus and pigtailed macaque. A total of 1395671, 1284929 and 1266348 perfect SSRs were mined, respectively. The most frequent perfect SSRs were mononucleotide SSRs. The most GC-content was in dinucleotide SSRs and the least was in the mononucleotide SSRs. Chromosome size was positively correlated with SSR number and negatively correlated with the relative frequency and density of SSRs. The GC content of chromosome SSRs were negatively correlated with relative frequency of SSRs and GC content of chromosome sequences. The features of microsatellite distribution in assembled genomes of the three species were greatly similar, which revealed that the distributional pattern of microsatellites is probably conservative in genus Macaca. The degenerated number of repeat motifs was found to be different in pentanucleotide and hexanucleotide repeats. Species-specific motifs for each macaque were significantly underrepresented. Overall, SSR frequencies of each chromosome in rhesus macaque were higher than in cynomolgus macaque. The maximum repeat times of mono- to pentanucleotide repeats in cynomolgus macaque was more than other two macaques. These results emphasize the genetic diversity and phylogenetic relationship of genus Macaca species. Our data will be beneficial for comparative genome mapping, understanding the distribution of SSRs and genome structure between these animal models, and provide a foundation for further development and identification of more macaque-specific SSRs.

Cross-Species Rhesus Cytomegalovirus Infection of Cynomolgus Macaques

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.

Impact of the Genetics and Source of Preclinical Safety Animal Models on Study Design, Results, and Interpretation

It has been long established that not only the species but also the strain and supplier of rodents used in preclinical safety studies can have a significant impact on the outcome of studies due to variability in their genetic background and thus spontaneous pathologic findings. In addition, local husbandry, housing, and other environmental conditions may have effects on the development and expression of comorbidities, particularly in longer-term or chronic studies. More recently, similar effects related to the source, including genetic and environmental variability, have been recognized in cynomolgus macaques ( Macaca fascicularis). The increased use of cynomolgus macaques from various sources of captive-bred animals (including nonnative, U.S./European Union-based breeding facilities or colonies) can affect study design and study results and outcome. It is important to acknowledge and understand the impact of this variability on the results and interpretation of research studies. This review includes recent examples where variability of preclinical animal models (rats and monkeys) affected the postmortem observations highlighting its relevance to study design or interpretation in safety studies.

Mitochondrial DNA and two Y-chromosome genes of common long-tailed macaques (Macaca fascicularis fascicularis) throughout Thailand and vicinity

Macaca fascicularis fascicularis is distributed over a wide area of Southeast Asia. Thailand is located at the center of their distribution range and is the bridge connecting the two biogeographic regions of Indochina and Sunda. However, only a few genetic studies have explored the macaques in this region. To shed some light on the evolutionary history of M. f. fascicularis, including hybridization with M. mulatta, M. f. fascicularis and M. mulatta samples of known origins throughout Thailand and the vicinity were analyzed by molecular phylogenetics using mitochondrial DNA (mtDNA), including the hypervariable region 1, and Y-chromosomal DNA, including SRY and TSPY genes. The mtDNA phylogenetic analysis divided M. f. fascicularis into five subclades (Insular Indonesia, Sundaic Thai Gulf, Vietnam, Sundaic Andaman sea coast, and Indochina) and revealed genetic differentiation between the two sides of the Thai peninsula, which had previously been reported as a single group of Malay peninsular macaques. From the estimated divergence time of the Sundaic Andaman sea coast subclade, it is proposed that after M. f. fascicularis dispersed throughout Southeast Asia, some populations on the south-easternmost Indochina (eastern Thailand, southern Cambodia and southern Vietnam at the present time) migrated south-westwards across the land bridge, which was exposed during the glacial period of the late Pleistocene epoch, to the southernmost Thailand/northern peninsular Malaysia. Then, some of them migrated north and south to colonize the Thai Andaman sea coast and northern Sumatra, respectively. The SRY-TSPY phylogenetic analysis suggested that male-mediated gene flow from M. mulatta southward to M. f. fascicularis was restricted south of, but close to, the Isthmus of Kra. There was a strong impact of the geographical factors in Thailand, such as the Isthmus of Kra, Nakhon Si Thammarat, and Phuket ranges and Sundaland, on M. f. fascicularis biogeography and their hybridization with M. mulatta.

Emergence and evolution of inter-specific segregating retrocopies in cynomolgus monkey (Macaca fascicularis) and rhesus macaque (Macaca mulatta)

Retroposition is an RNA-mediated mechanism to generate gene duplication, and is believed to play an important role in genome evolution and phenotypic adaptation in various species including primates. Previous studies suggested an elevated rate of recent retroposition in the rhesus macaque genome. To better understand the impact of retroposition on macaque species which have undergone an adaptive radiation approximately 3–6 million years ago, we developed a bioinformatics pipeline to identify recently derived retrocopies in cynomolgus monkeys. As a result, we identified seven experimentally validated young retrocopies, all of which are polymorphic in cynomolgus monkeys. Unexpectedly, five of them are also present in rhesus monkeys and are still segregating. Molecular evolutionary analysis indicates that the observed inter-specific polymorphism is attribute to ancestral polymorphism. Further population genetics analysis provided strong evidence of balancing selection on at least one case (Crab-eating monkey retrocopy 6, or CER6) in both species. CER6 is in adjacent with an immunoglobulin related gene and may be involved in host-pathogen interaction, a well-known target of balancing selection. Altogether, our data support that retroposition is an important force to shape genome evolution and species adaptation.

An SNP marker at the STAT6 locus can identify the hybrids between rhesus (Macaca mulatta) and long-tailed macaques (M. fascicularis) in Thailand: a rapid and simple screening method and its application

A polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay was developed to genetically discriminate rhesus (Macaca mulatta) macaques from long-tailed (M. fascicularis) macaques. The 745 bp PCR amplicon of the STAT6 locus that spans a potentially species-diagnostic single nucleotide polymorphism (SNP) marker was digested with ApaI and gel electrophoresed to give (1) two (234 and 511 bp), (2) one (745 bp) and (3) three (234, 511 and 745 bp) band patterns that correspond to the genotypes G/G (long-tailed macaque specific homozygote), A/A (rhesus macaque specific homozygote) and A/G (hybrid specific heterozygote), respectively. The diagnostic robustness and efficiency of this PCR-RFLP assay was tested on wild rhesus and long-tailed macaques inhabiting Thailand and a known hybrid population. The Indochinese and Sundaic long-tailed macaque samples (n = 18) all showed a homozygous G/G pattern, while the Indochinese rhesus macaques (n = 10) all showed a homozygous A/A pattern. The rhesus/long-tailed hybrid population at Khao Khieow Open Zoo, which resulted from an introduced group of rhesus macaques that hybridized with the indigenous long-tailed macaques about 20 years ago, revealed 47% (56/118 samples analyzed) with the heterogenous A/G genotype. In addition, the frequency of the rhesus-specific allele A significantly decreased in the hybrid population during 2006-2014, where a strong association between the STAT6 genotype and the morphology of the individuals was detected. In conclusion, a robust PCR-RFLP assay allows a simple, effective and inexpensive approach, in particular for field studies, to assess hybrid individuals between rhesus and long-tailed macaques. Although this assay cannot conclusively identify all the hybrids over two or more generations, it at least can allow the evaluation of the process of hybridization, and so it is applicable to the assessment of the status of natural or anthropogenic hybridization between the two species across their geographic range.

Assessment of the hybridization between rhesus (Macaca mulatta) and long-tailed macaques (M. fascicularis) based on morphological characters

OBJECTIVES

Hybridization between rhesus (Macaca mulatta) and long-tailed (M. fascicularis) macaques has become a focal point of interest. The majority of such studies have evaluated their genetics, but not their morphological characters.

MATERIALS AND METHODS

We analyzed morphological characters of eight free-ranging populations of Indochinese rhesus and long-tailed macaques distributed at the proposed hybrid zone (15.75-21.58° N) in comparison with one population each of Chinese and Indian-derived rhesus macaques and three populations of Sundaic long-tailed macaques.

RESULTS

Chinese and Indian rhesus macaques had a heavier body mass, longer crown-rump length, shorter relative facial length and relative tail length, and a greater contrast of reddish and yellowish dorsal pelage color than the Sundaic long-tailed macaques for which the latter three parameters could be used to visually discriminate between the two species. Although the morphological characters of Indochinese rhesus and long-tailed macaques were intermediate between the Chinese/Indian rhesus and Sundaic long-tailed macaques, they were more similar to their respective conspecifics. The species-specific characters of a shorter tail (<70%) and a bipartite pelage color pattern were retained in the Indochinese rhesus macaques while the longer tail (>90%) and no bipartite pattern was found in the Indochinese long-tailed macaques. No morphological cline was observed across the species and the variations were abrupt to some extent.

DISCUSSION

The hybridization between rhesus and long-tailed macaques may be results of multiple contacts and isolations over a long period of time, thus their evolutionary history should not be drawn solely by genetic or morphological analysis.

Examining the species-specificity of rhesus macaque cytomegalovirus (RhCMV) in cynomolgus macaques

Cytomegalovirus (CMV) is a highly species-specific virus that has co-evolved with its host over millions of years and thus restricting cross-species infection. To examine the extent to which host restriction may prevent cross-species research between closely related non-human primates, we evaluated experimental infection of cynomolgus macaques with a recombinant rhesus macaque-derived CMV (RhCMV-eGFP). Twelve cynomolgus macaques were randomly allocated to three groups: one experimental group (RhCMV-eGFP) and two control groups (UV-inactivated RhCMV-eGFP or media alone). The animals were given two subcutaneous inoculations at week 0 and week 8, and a subset of animals received an intravenous inoculation at week 23. No overt clinical or haematological changes were observed and PBMCs isolated from RhCMV-eGFP inoculated animals had comparable eGFP- and IE-1-specific cellular responses to the control animals. Following inoculation with RhCMV-eGFP, we were unable to detect evidence of infection in any blood or tissue samples up to 4 years post-inoculation, using sensitive viral co-culture, qPCR, and Western blot assays. Co-culture of urine and saliva samples demonstrated the presence of endogenous cynomolgus CMV (CyCMV) cytopathic effect, however no concomitant eGFP expression was observed. The absence of detectable RhCMV-eGFP suggests that the CyCMV-seropositive cynomolgus macaques were not productively infected with RhCMV-eGFP under these inoculation conditions. In a continued effort to develop CMV as a viral vector for an HIV/SIV vaccine, these studies demonstrate that CMV is highly restricted to its host species and can be highly affected by laboratory cell culture. Consideration of the differences between lab-adapted and primary viruses with respect to species range and cell tropism should be a priority in evaluating CMV as vaccine vector for HIV or other pathogens at the preclinical development stage.

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.

Mitogenomic phylogeny of the common long-tailed macaque (Macaca fascicularis fascicularis)

Background

Long-tailed macaques (Macaca fascicularis) are an important model species in biomedical research and reliable knowledge about their evolutionary history is essential for biomedical inferences. Ten subspecies have been recognized, of which most are restricted to small islands of Southeast Asia. In contrast, the common long-tailed macaque (M. f. fascicularis) is distributed over large parts of the Southeast Asian mainland and the Sundaland region. To shed more light on the phylogeny of M. f. fascicularis, we sequenced complete mitochondrial (mtDNA) genomes of 40 individuals from all over the taxon’s range, either by classical PCR-amplification and Sanger sequencing or by DNA-capture and high-throughput sequencing.

Results

Both laboratory approaches yielded complete mtDNA genomes from M. f. fascicularis with high accuracy and/or coverage. According to our phylogenetic reconstructions, M. f. fascicularis initially diverged into two clades 1.70 million years ago (Ma), with one including haplotypes from mainland Southeast Asia, the Malay Peninsula and North Sumatra (Clade A) and the other, haplotypes from the islands of Bangka, Java, Borneo, Timor, and the Philippines (Clade B). The three geographical populations of Clade A appear as paraphyletic groups, while local populations of Clade B form monophyletic clades with the exception of a Philippine individual which is nested within the Borneo clade. Further, in Clade B the branching pattern among main clades/lineages remains largely unresolved, most likely due to their relatively rapid diversification 0.93-0.84 Ma.

Conclusions

Both laboratory methods have proven to be powerful to generate complete mtDNA genome data with similarly high accuracy, with the DNA-capture and high-throughput sequencing approach as the most promising and only practical option to obtain such data from highly degraded DNA, in time and with relatively low costs. The application of complete mtDNA genomes yields new insights into the evolutionary history of M. f. fascicularis by providing a more robust phylogeny and more reliable divergence age estimations than earlier studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1437-0) contains supplementary material, which is available to authorized users.

Genotyping of non-human primate models: perspectives and challenges for the implementation of the &quot;three R's&quot;

Abstract. Although non-human primates (NHPs) represent only a minor fraction of all animals used in biomedical research, there is a continuous effort to further reduce, refine and replace research with NHPs in accordance with the principles of the three R's. Most of the NHP model species are genetically highly diverse, and significant variation occurs among populations of different geographic origins, particularly in macaques. Since such differences can considerably affect the outcome of biomedical experiments, genotyping represents a promising tool to refine research approaches and to reduce the number of NHPs in biomedical research. Accordingly, the European Primate Network (EUPRIM-Net) developed an anonymous online survey to evaluate possibilities and potential hindrances for the application of genotyping in NHP research. On the one hand, our results point to the importance of genetic variation in NHPs and the need to consider the genetic background for future research approaches. On the other hand, our survey identified several hindrances and limiting factors for the application of genotyping and its incorporation in research, primate husbandry and breeding. We provide some fundamental recommendations on how to meet these challenges and how genotyping can be efficiently used to refine NHP research and to reduce the number of NHPs in biomedical research in the long term.

A large-scale SNP-based genomic admixture analysis of the captive rhesus macaque colony at the California National Primate Research Center

Some breeding facilities in the United States have crossbred Chinese and Indian rhesus macaque (Macaca mulatta) founders either purposefully or inadvertently. Genetic variation that reflects geographic origins among research subjects has the potential to influence experimental outcomes. The use of animals from different geographic regions, their hybrids, and animals of varying degrees of kinship in an experiment can obscure treatment effects under study because high interanimal genetic variance can increase phenotypic variance among the research subjects. The intent of this study, based on a broad genomic analysis of 2,808 single nucleotide polymorphisms (SNPs), is to ensure that only animals estimated to be of pure Indian or Chinese ancestry, based on both demographic and genetic information, are used as sources of infants for derivation and expansion of the California National Primate Research Center's (CNPRC) super-Specific Pathogen Free (SSPF) rhesus macaque colony. Studies of short tandem repeats (STRs) in Indian and Chinese rhesus macaques have reported that heterozygosity of STRs is higher in Chinese rhesus macaques than in Indian rhesus macaques. The present study shows that heterozygosity of SNPs is actually higher in Indian than in Chinese rhesus macaques and that the Chinese SSPF rhesus macaque colony is far less differentiated from their founders compared to the Indian-origin animals. The results also reveal no evidence of recent gene flow from long-tailed and pig-tailed macaques into the source populations of the SSPF rhesus macaques. This study indicates that many of the long-tailed macaques held in the CNPRC are closely related individuals. Most polymorphisms shared among the captive rhesus, long-tailed, and pig-tailed macaques likely predate the divergence among these groups.

Utility, limitations, and future of non-human primates for dengue research and vaccine development

Dengue is considered the most important emerging, human arboviruses, with worldwide distribution in the tropics. Unfortunately, there are no licensed dengue vaccines available or specific anti-viral drugs. The development of a dengue vaccine faces unique challenges. The four serotypes co-circulate in endemic areas, and pre-existing immunity to one serotype does not protect against infection with other serotypes, and actually may enhance severity of disease. One foremost constraint to test the efficacy of a dengue vaccine is the lack of an animal model that adequately recapitulates the clinical manifestations of a dengue infection in humans. In spite of this limitation, non-human primates (NHP) are considered the best available animal model to evaluate dengue vaccine candidates due to their genetic relatedness to humans and their ability to develop a viremia upon infection and a robust immune response similar to that in humans. Therefore, most dengue vaccines candidates are tested in primates before going into clinical trials. In this article, we present a comprehensive review of published studies on dengue vaccine evaluations using the NHP model, and discuss critical parameters affecting the usefulness of the model. In the light of recent clinical data, we assess the ability of the NHP model to predict immunological parameters of vaccine performances in humans and discuss parameters that should be further examined as potential correlates of protection. Finally, we propose some guidelines toward a more standardized use of the model to maximize its usefulness and to better compare the performance of vaccine candidates from different research groups.

A genetic comparison of two alleged subspecies of Philippine cynomolgus macaques

Two subspecies of cynomolgus macaques (Macaca fascicularis) are alleged to co-exist in the Philippines, M. f. philippensis in the north and M. f. fascicularis in the south. However, genetic differences between the cynomolgus macaques in the two regions have never been studied to document the propriety of their subspecies status. We genotyped samples of cynomolgus macaques from Batangas in southwestern Luzon and Zamboanga in southwestern Mindanao for 15 short tandem repeat (STR) loci and sequenced an 835 bp fragment of the mtDNA of these animals. The STR genotypes were compared with those of cynomolgus macaques from southern Sumatra, Singapore, Mauritius and Cambodia, and the mtDNA sequences of both Philippine populations were compared with those of cynomolgus macaques from southern Sumatra, Indonesia and Sarawak, Malaysia. We conducted STRUCTURE and PCA analyses based on the STRs and constructed a median joining network based on the mtDNA sequences. The Philippine population from Batangas exhibited much less genetic diversity and greater genetic divergence from all other populations, including the Philippine population from Zamboanga. Sequences from both Batangas and Zamboanga were most closely related to two different mtDNA haplotypes from Sarawak from which they are apparently derived. Those from Zamboanga were more recently derived than those from Batangas, consistent with their later arrival in the Philippines. However, clustering analyses do not support a sufficient genetic distinction of cynomolgus macaques from Batangas from other regional populations assigned to subspecies M. f. fascicularis to warrant the subspecies distinction M. f. philippensis.

Development and validation of a SNP-based assay for inferring the genetic ancestry of rhesus macaques (Macaca mulatta)

Rhesus macaques (Macaca mulatta) are an important primate model species in several areas of biomedical research. The wide geographic distribution of this species has led to significant genetic differentiation among local and regional populations. These regional differences can be important factors in the selection of the most appropriate subjects for particular research studies, as animals from different populations can respond differently to the same experimental treatment. Consequently, it is valuable to confirm the ancestry of individual rhesus monkeys from geographically distinct populations. Using DNA samples obtained from rhesus macaques from six National Primate Research Centers, we tested a set of 384 potential ancestry informative single nucleotide polymorphisms (SNPs) and identified a final panel of 91 SNPs that can reliably distinguish Indian-origin from Chinese-origin rhesus monkeys. This genetic test can be used to determine the ancestral origin of animals and to detect individuals that are hybrids between these two regional populations. To demonstrate use of the SNP panel, we investigated the ancestry of 480 animals from the Yerkes NPRC (YNPRC) for which the colony records were insufficient to clearly establish ancestry. Three of the YNPRC animals tested were determined to be hybrids. This SNP ancestry tool will be useful to researchers, colony managers, and others who wish to evaluate the ancestral origin of individual rhesus macaques, and therefore will facilitate more effective and efficient use of these animals in biomedical research.

Genetic diversity and population structure of long-tailed macaque (Macaca fascicularis) populations in Peninsular Malaysia

BACKGROUND

The genetic diversity and structure of long-tailed macaques (Macaca fascicularis) in Peninsular Malaysia, a widely used non-human primate species in biomedical research, have not been thoroughly characterized.

METHODS

Thirteen sites of wild populations of long-tailed macaques representing six states were sampled and analyzed with 18 STR markers.

RESULTS

The Sunggala and Penang Island populations showed the highest genetic diversity estimates, while the Jerejak Island population was the most genetically discrete due to isolation from the mainland shelf. Concordant with pairwise F(st) estimates, STRUCTURE analyses of the seven PCA-correlated clusters revealed low to moderate differentiation among the sampling sites. No association between geographic and genetic distances exists, suggesting that the study sites, including island study sites, are genetically if not geographically contiguous.

CONCLUSIONS

The status of the genetic structure and composition of long-tailed macaque populations require further scrutiny to develop this species as an important animal model in biomedical research.

Whole-genome sequencing of tibetan macaque (Macaca Thibetana) provides new insight into the macaque evolutionary history

Macaques are the most widely distributed nonhuman primates and used as animal models in biomedical research. The availability of full-genome sequences from them would be essential to both biomedical and primate evolutionary studies. Previous studies have reported whole-genome sequences from rhesus macaque (Macaca mulatta) and cynomolgus macaque (M. fascicularis, CE), both of which belong to the fascicularis group. Here, we present a 37-fold coverage genome sequence of the Tibetan macaque (M. thibetana; TM). TM is an endemic species to China belonging to the sinica group. On the basis of mapping to the rhesus macaque genome, we identified approximately 11.9 million single-nucleotide variants), of which 3.9 million were TM specific, as assessed by comparison two Chinese rhesus macaques (CR) and two CE genomes. Some genes carried TM-specific homozygous nonsynonymous variants (TSHNVs), which were scored as deleterious in human by both PolyPhen-2 and SIFT (Sorting Tolerant From Intolerant) and were enriched in the eye disease genes. In total, 273 immune response and disease-related genes carried at least one TSHNV. The heterozygosity rates of two CRs (0.002617 and 0.002612) and two CEs (0.003004 and 0.003179) were approximately three times higher than that of TM (0.000898). Polymerase chain reaction resequencing of 18 TM individuals showed that 29 TSHNVs exhibited high allele frequencies, thus confirming their low heterozygosity. Genome-wide genetic divergence analysis demonstrated that TM was more closely related to CR than to CE. We further detected unusual low divergence regions between TM and CR. In addition, after applying statistical criteria to detect putative introgression regions (PIRs) in the TM genome, up to 239,620 kb PIRs (8.84% of the genome) were identified. Given that TM and CR have overlapping geographical distributions, had the same refuge during the Middle Pleistocene, and show similar mating behaviors, it is highly likely that there was an ancient introgression event between them. Moreover, demographic inferences revealed that TM exhibited a similar demographic history as other macaques until 0.5 Ma, but then it maintained a lower effective population size until present time. Our study has provided new insight into the macaque evolutionary history, confirming hybridization events between macaque species groups based on genome-wide data.

Differential response of the cynomolgus macaque gut microbiota to Shigella infection

Little is known about the role of gut microbiota in response to live oral vaccines against enteric pathogens. We examined the effect of immunization with an oral live-attenuated Shigella dysenteriae 1 vaccine and challenge with wild-type S. dysenteriae 1 on the fecal microbiota of cynomolgus macaques using 16 S rRNA analysis of fecal samples. Multi-dimensional cluster analysis identified different bacterial community types within macaques from geographically distinct locations. The fecal microbiota of Mauritian macaques, observed to be genetically distinct, harbored a high-diversity community and responded differently to Shigella immunization, as well as challenge compared to the microbiota in non-Mauritian macaques. While both macaque populations exhibited anti-Shigella antibody responses, clinical shigellosis was observed only among non-Mauritian macaques. These studies highlight the importance of further investigation into the possible protective role of the microbiota against enteric pathogens and consideration of host genetic backgrounds in conducting vaccine studies.

The genetic composition of populations of cynomolgus macaques (Macaca fascicularis) used in biomedical research

BACKGROUND

The genetic composition of cynomolgus macaques used in biomedical research is not as well-characterized as that of rhesus macaques.

METHODS

Populations of cynomolgus macaques from Sumatra, Corregidor, Mauritius, Singapore, Cambodia, and Zamboanga were analyzed using 24 STRs.

RESULTS

The Sumatran and Cambodian populations exhibited the highest allelic diversity, while the Mauritian population exhibited the lowest. Sumatran cynomolgus macaques were the most genetically similar to all others, consistent with an Indonesian origin of the species. The high diversity among Cambodian animals may result from interbreeding with rhesus macaques. The Philippine and Mauritian samples were the most divergent from other populations, the former due to separation from the Sunda Shelf by deepwater and the latter due to anthropogenic translocation and extreme founder effects.

CONCLUSIONS

Investigators should verify their research subjects' origin, ancestry, and pedigree to minimize risks to biomedical experimentation from genetic variance stemming from close kinship and mixed ancestry as these can obscure treatment effects.

Characterization of the major histocompatibility complex class I A alleles in cynomolgus macaques of Vietnamese origin

Cynomolgus macaques (Macaca fascicularis, Mafa) have emerged as an important animal model for infectious disease and transplantation research. Extensive characterization of their major histocompatibility complex (MHC) polymorphism regions therefore becomes urgently required. In this study, we identified 41 MHC class I A nucleotide sequences in 34 unrelated cynomolgus macaques of Vietnamese origin farmed in Southern China, including eight novel Mafa-A sequences. We found two sequences with perfect identity and six sequences with close similarity to previously defined MHC class I alleles from other populations, especially from Indonesian-origin macaques. We also found three Vietnamese-origin cynomolgus macaque MHC class I sequences for which the predicted protein sequences identical throughout their B and F binding pockets to Mamu-A1*001:01 and Mamu-A3*13:03, respectively. This is important because Mamu-A1*001:01 and Mamu-A3*13:03 are associated with longer survival and lower set-point viral load in simian immunodeficiency virus (SIV)-infected rhesus monkeys. These findings have implications for the evolutionary history of Vietnamese-origin cynomolgus macaque as well as for the use of this model in SIV/SHIV (a virus combining parts of the HIV and SIV genomes) research.

Population and landscape genetics of an introduced species (M. fascicularis) on the island of Mauritius

The cynomolgus macaque, Macaca fascicularis, was introduced onto the island of Mauritius in the early 17(th) century. The species experienced explosive population growth, and currently exists at high population densities. Anecdotes collected from nonhuman primate trappers on the island of Mauritius allege that animals from the northern portion of the island are larger in body size than and superior in condition to their conspecifics in the south. Although previous genetic studies have reported Mauritian cynomolgus macaques to be panmictic, the individuals included in these studies were either from the southern/central or an unknown portion of the island. In this study, we sampled individuals broadly throughout the entire island of Mauritius and used spatial principle component analysis to measure the fine-scale correlation between geographic and genetic distance in this population. A stronger correlation between geographic and genetic distance was found among animals in the north than in those in the southern and central portions of the island. We found no difference in body weight between the two groups, despite anecdotal evidence to the contrary. We hypothesize that the increased genetic structure among populations in the north is related to a reduction in dispersal distance brought about by human habitation and tourist infrastructure, but too recent to have produced true genetic differentiation.

Preliminary observations of MHC class I A region polymorphism in three populations of Chinese-origin rhesus macaques

Rhesus macaques are an animal model for the study of a variety of human diseases. The Chinese rhesus macaques have been widely used in biomedical research in recent years. However, the polymorphism of major histocompatibility complex (MHC) class I A region among different local populations of Chinese rhesus macaques has never been investigated. In this study, we identified 46 Mamu-A alleles by cDNA cloning and sequencing on a cohort of 53 Chinese rhesus monkeys including Zhiming, Chuanxi, and Fujian populations, of which 5 were first reported in rhesus monkeys. The frequencies of alleles were identified for each population. The result suggests that the repertoire of allelic variants of MHC class I A region found in different populations of Chinese macaques is largely non-overlapping. The frequencies of alleles and the popular allele are also different for different populations. PCR-SSP experiment further confirms the different frequencies of two alleles, Mamu-A*026:01 and Mamu-A*022:01, in additional 99 Zhiming monkeys and 191 Chuanxi monkeys. Our findings have important practical implications in that the origin of the individuals and the genetic polymorphism of the monkeys need to be considered at the level of local populations for Chinese rhesus monkeys in biomedical research. Further immunogenetic work is needed to investigate the MHC polymorphism among different populations of Chinese rhesus macaques and to reveal the functional implication of such polymorphism and disease outcome correlations.

Molecular ABO phenotyping in cynomolgus macaques using real-time quantitative PCR

Macaques are commonly used in biomedical research as animal models of human disease. The ABO phenotype of donors and recipients plays an important role in the success of transplantation and stem cell research of both human and macaque tissue. Traditional serological methods for ABO phenotyping can be time consuming, provide ambiguous results and/or require tissue that is unavailable or unsuitable. We developed a novel method to detect the A, B, and AB phenotypes of macaques using real-time quantitative polymerase chain reaction. This method enables the simple and rapid screening of these phenotypes in macaques without the need for fresh blood or saliva. This study reports the distribution of the A, B, and AB phenotypes of captive cynomolgus macaques that, while regionally variable, closely resembles that of rhesus macaques. Blood group B, as in rhesus macaques, predominates in cynomolgus macaques and its frequency distribution leads to a probability of major incompatibility of 41%. No silencing mutations have been identified in exon 6 or 7 in macaques that could be responsible for the O phenotype, that, although rare, have been reported. The excess homozygosity of rhesus and cynomolgus macaque genotypes in this study, that assumes the absence of the O allele, suggests the possibility of some mechanism preventing the expression of the A and B transferases.

Population genetics of the Washington National Primate Research Center's (WaNPRC) captive pigtailed macaque (Macaca nemestrina) population

Pigtailed macaques (Macaca nemestrina) provide an important model for biomedical research on human disease and for studying the evolution of primate behavior. The genetic structure of captive populations of pigtailed macaques is not as well described as that of captive rhesus (M. mulatta) or cynomolgus (M. fascicularis) macaques. The Washington National Primate Research Center houses the largest captive colony of pigtailed macaques located in several different housing facilities. Based on genotypes of 18 microsatellite (short tandem repeat [STR]) loci, these pigtailed macaques are more genetically diverse than captive rhesus macaques and exhibit relatively low levels of inbreeding. Colony genetic management facilitates the maintenance of genetic variability without compromising production goals of a breeding facility. The periodic introduction of new founders from specific sources to separate housing facilities at different times influenced the colony's genetic structure over time and space markedly but did not alter its genetic diversity significantly. Changes in genetic structure over time were predominantly due to the inclusion of animals from the Yerkes National Primate Research Center in the original colony and after 2005. Strategies to equalize founder representation in the colony have maximized the representation of the founders' genomes in the extant population. Were exchange of animals among the facilities increased, further differentiation could be avoided. The use of highly differentiated animals may confound interpretations of phenotypic differences due to the inflation of the genetic contribution to phenotypic variance of heritable traits.

Whole-genome sequencing and analysis of the Malaysian cynomolgus macaque (Macaca fascicularis) genome

BACKGROUND

The genetic background of the cynomolgus macaque (Macaca fascicularis) is made complex by the high genetic diversity, population structure, and gene introgression from the closely related rhesus macaque (Macaca mulatta). Herein we report the whole-genome sequence of a Malaysian cynomolgus macaque male with more than 40-fold coverage, which was determined using a resequencing method based on the Indian rhesus macaque genome.

RESULTS

We identified approximately 9.7 million single nucleotide variants (SNVs) between the Malaysian cynomolgus and the Indian rhesus macaque genomes. Compared with humans, a smaller nonsynonymous/synonymous SNV ratio in the cynomolgus macaque suggests more effective removal of slightly deleterious mutations. Comparison of two cynomolgus (Malaysian and Vietnamese) and two rhesus (Indian and Chinese) macaque genomes, including previously published macaque genomes, suggests that Indochinese cynomolgus macaques have been more affected by gene introgression from rhesus macaques. We further identified 60 nonsynonymous SNVs that completely differentiated the cynomolgus and rhesus macaque genomes, and that could be important candidate variants for determining species-specific responses to drugs and pathogens. The demographic inference using the genome sequence data revealed that Malaysian cynomolgus macaques have experienced at least three population bottlenecks.

CONCLUSIONS

This list of whole-genome SNVs will be useful for many future applications, such as an array-based genotyping system for macaque individuals. High-quality whole-genome sequencing of the cynomolgus macaque genome may aid studies on finding genetic differences that are responsible for phenotypic diversity in macaques and may help control genetic backgrounds among individuals.

Microsatellite variation in two subspecies of cynomolgus monkeys (Macaca fascicularis)

To estimate the genetic variability of two subspecies of cynomolgus monkeys (Macaca fascicularis fascicularis and M. f. aurea) using microsatellite markers, 26 microsatellite markers were selected from previous reports. Seventeen markers showed high polymorphism in a subset of monkeys and were used for the assessment of genetic diversity in the larger sample. The effective number of alleles, the polymorphism information content (PIC) and the expected heterozygosity of M. f. aurea monkeys were all statistically significantly higher than those of M. f. fascicularis monkeys (P < 0.05), suggesting the M. f. aurea monkeys had a higher degree of genetic variation than the M. f. fascicularis monkeys. Substantial differences in allele distribution were also detected between the two subspecies of cynomolgus monkeys. Private alleles restricted to the M. f. fascicularis or the M. f. aurea monkeys were found throughout the selected 17 loci. These private alleles may allow the discrimination of the two subspecies of cynomolgus monkeys. The selected markers could also be used to estimate the genetic variation for other subspecies of cynomolgus monkeys. Further work using additional animals obtained from native or independent sources will be important for a more complete understanding of the genetic differences between these two subgroups.

Genome-based analysis of the nonhuman primate Macaca fascicularis as a model for drug safety assessment

The long-tailed macaque, also referred to as cynomolgus monkey (Macaca fascicularis), is one of the most important nonhuman primate animal models in basic and applied biomedical research. To improve the predictive power of primate experiments for humans, we determined the genome sequence of a Macaca fascicularis female of Mauritian origin using a whole-genome shotgun sequencing approach. We applied a template switch strategy that uses either the rhesus or the human genome to assemble sequence reads. The sixfold sequence coverage of the draft genome sequence enabled discovery of about 2.1 million potential single-nucleotide polymorphisms based on occurrence of a dimorphic nucleotide at a given position in the genome sequence. Homology-based annotation allowed us to identify 17,387 orthologs of human protein-coding genes in the M. fascicularis draft genome, and the predicted transcripts enabled the design of a M. fascicularis-specific gene expression microarray. Using liver samples from 36 individuals of different geographic origin we identified 718 genes with highly variable expression in liver, whereas the majority of the transcriptome shows relatively stable and comparable expression. Knowledge of the M. fascicularis draft genome is an important contribution to both the use of this animal in disease models and the safety assessment of drugs and their metabolites. In particular, this information allows high-resolution genotyping and microarray-based gene-expression profiling for animal stratification, thereby allowing the use of well-characterized animals for safety testing. Finally, the genome sequence presented here is a significant contribution to the global "3R" animal welfare initiative, which has the goal to reduce, refine, and replace animal experiments.

Variable prevalence and functional diversity of the antiretroviral restriction factor TRIMCyp in Macaca fascicularis

The retroviral restriction factor TRIMCyp, derived from the TRIM5 gene, blocks replication at a postentry step. TRIMCyp has so far been found in four species of Asian macaques, Macaca fascicularis, M. mulatta, M. nemestrina, and M. leonina. M. fascicularis is commonly used as a model for AIDS research, but TRIMCyp has not been analyzed in detail in this species. We analyzed the prevalence of TRIMCyp in samples from Indonesia, Indochina, the Philippines, and Mauritius. We found that TRIMCyp is present at a higher frequency in Indonesian than in Indochinese M. fascicularis macaques and is also present in samples from the Philippines. TRIMCyp is absent in Mauritian M. fascicularis macaques. We then analyzed the restriction specificity of TRIMCyp derived from three animals of Indonesian origin. One allele, like the prototypic TRIMCyp alleles described for M. mulatta and M. nemestrina, restricts human immunodeficiency virus type 2 (HIV-2) and feline immunodeficiency virus (FIV) but not HIV-1. The others restrict HIV-1 and FIV but not HIV-2. Mutagenesis studies confirmed that polymorphisms at amino acid residues 369 and 446 in TRIMCyp (or residues 66 and 143 in the cyclophilin A [CypA] domain) confer restriction specificity. Additionally, we identified a polymorphism in the coiled-coil domain that appears to affect TRIMCyp expression or stability. Taken together, these data show that M. fascicularis has the most diverse array of TRIM5 restriction factors described for any primate species to date. These findings are relevant to our understanding of the evolution of retroviral restriction factors and the use of M. fascicularis models in AIDS research.