Genetic characterization of rhesus macaques (Macaca mulatta) in Nepal

The natural history of model organisms: the rhesus macaque as a success story of the Anthropocene

Of all the non-human primate species studied by researchers, the rhesus macaque (Macaca mulatta) is likely the most widely used across biological disciplines. Rhesus macaques have thrived during the Anthropocene and now have the largest natural range of any non-human primate. They are highly social, exhibit marked genetic diversity, and display remarkable niche flexibility (which allows them to live in a range of habitats and survive on a variety of diets). These characteristics mean that rhesus macaques are well-suited for understanding the links between sociality, health and fitness, and also for investigating intra-specific variation, adaptation and other topics in evolutionary ecology.

Staphylococcus aureus and Methicillin Resistant S. aureus in Nepalese Primates: Resistance to Antimicrobials, Virulence, and Genetic Lineages

Staphylococcus aureus is a ubiquitous pathogen and colonizer in humans and animals. There are few studies on the molecular epidemiology of S. aureus in wild monkeys and apes. S. aureus carriage in rhesus macaques (Macaca mulatta) and Assam macaques (Macaca assamensis) is a species that has not previously been sampled and lives in remote environments with limited human contact. Forty Staphylococcus aureus isolates including 33 methicillin-susceptible S. aureus (MSSA) and seven methicillin-resistant S. aureus (MRSA) were characterized. Thirty-four isolates were from rhesus macaques and six isolates (five MSSA, one MRSA) were from Assam macaques. Isolates were characterized using StaphyType DNA microarrays. Five of the MRSA including one from Assam macaque were CC22 MRSA-IV (PVL+/tst+), which is a strain previously identified in Nepalese rhesus. One MRSA each were CC6 MRSA-IV and CC772 MRSA-V (PVL+). One MSSA each belonged to CC15, CC96, and CC2990. Six MRSA isolates carried the blaZ, while ten known CC isolates (seven MRSA, three MSSA) carried a variety of genes including aacA-aphD, aphA3, erm(C), mph(C), dfrA, msrA, and/or sat genes. The other 30 MSSA isolates belonged to 17 novel clonal complexes, carried no antibiotic resistance genes, lacked Panton–Valentine Leukocidin (PVL), and most examined exotoxin genes. Four clonal complexes carried egc enterotoxin genes, and four harbored edinB, which is an exfoliative toxin homologue.

Mitochondrial DNA analyses and ecological niche modeling reveal post-LGM expansion of the Assam macaque (Macaca assamensis) in the foothills of Nepal Himalaya

Genetic diversity of a species is influenced by multiple factors, including the Quaternary glacial-interglacial cycles and geophysical barriers. Such factors are not yet well documented for fauna from the southern border of the Himalayan region. This study used mitochondrial DNA (mtDNA) sequences and ecological niche modeling (ENM) to explore how the late Pleistocene climatic fluctuations and complex geography of the Himalayan region have shaped genetic diversity, population genetic structure, and demographic history of the Nepalese population of Assam macaques (Macaca assamensis) in the Himalayan foothills. A total of 277 fecal samples were collected from 39 wild troops over almost the entire distribution of the species in Nepal. The mtDNA fragment encompassing the complete control region (1121 bp) was recovered from 208 samples, thus defining 54 haplotypes. Results showed low nucleotide diversity (0.0075 ± SD 0.0001) but high haplotype diversity (0.965 ± SD 0.004). The mtDNA sequences revealed a shallow population genetic structure with a moderate but statistically significant effect of isolation by distance. Demographic history analyses using mtDNA sequences suggested a post-pleistocene population expansion. Paleodistribution reconstruction projected that the potential habitat of the Assam macaque was confined to the lower elevations of central Nepal during the Last Glacial Maximum. With the onset of the Holocene climatic optimum, the glacial refugia population experienced eastward range expansion to higher elevations. We conclude that the low genetic diversity and shallow population genetic structure of the Assam macaque population in the Nepal Himalaya region are the consequence of recent demographic and spatial expansion.

Performing monkeys of Bangladesh: characterizing their source and genetic variation

The acquisition and training of monkeys to perform is a centuries-old tradition in South Asia, resulting in a large number of rhesus macaques kept in captivity for this purpose. The performing monkeys are reportedly collected from free-ranging populations, and may escape from their owners or may be released into other populations. In order to determine whether this tradition involving the acquisition and movement of animals has influenced the population structure of free-ranging rhesus macaques in Bangladesh, we first characterized the source of these monkeys. Biological samples from 65 performing macaques collected between January 2010 and August 2013 were analyzed for genetic variation using 716 base pairs of mitochondrial DNA. Performing monkey sequences were compared with those of free-ranging rhesus macaque populations in Bangladesh, India and Myanmar. Forty-five haplotypes with 116 (16 %) polymorphic nucleotide sites were detected among the performing monkeys. As for the free-ranging rhesus population, most of the substitutions (89 %) were transitions, and no indels (insertion/deletion) were observed. The estimate of the mean number of pair-wise differences for the performing monkey population was 10.1264 ± 4.686, compared to 14.076 ± 6.363 for the free-ranging population. Fifteen free-ranging rhesus macaque populations were identified as the source of performing monkeys in Bangladesh; several of these populations were from areas where active provisioning has resulted in a large number of macaques. The collection of performing monkeys from India was also evident.

Macaque species susceptibility to simian immunodeficiency virus: increased incidence of SIV central nervous system disease in pigtailed macaques versus rhesus macaques

Immune pressure exerted by MHC class I-restricted cytotoxic T cells drives the development of viral escape mutations, thereby regulating HIV disease progression. Nonetheless, the relationship between host immunity and HIV central nervous system (CNS) disease remains poorly understood. The simian immunodeficiency virus (SIV) macaque model recapitulates key features of HIV infection including development of AIDS and CNS disease. To investigate cell-mediated immunity regulating SIV CNS disease progression, we compared the incidence of SIV encephalitis and the influence of MHC class I allele expression on the development of CNS disease in rhesus macaques (Macaca mulatta) versus pigtailed macaques (Macaca nemestrina). After inoculation with the immunosuppressive swarm SIV/DeltaB670 and the neurovirulent molecular clone SIV/17E-Fr, pigtailed macaques progressed more rapidly to AIDS, had higher plasma and cerebrospinal fluid (CSF) viral loads, and were more likely to progress to SIV-associated encephalitis (SIVE) compared to rhesus macaques. In addition, MHC class I alleles were neuroprotective in both species (Mamu-A*001 in rhesus macaques and Mane-A1*084:01:01 in pigtailed macaques); animals expressing these alleles were less likely to develop SIV encephalitis and correspondingly had lower viral replication in the brain. Species-specific differences in susceptibility to SIV disease demonstrated that cell mediated immune responses are critical to SIV CNS disease progression.

Diversity and molecular phylogeny of mitochondrial DNA of rhesus macaques (Macaca mulatta) in Bangladesh

While studies of rhesus macaques (Macaca mulatta) in the eastern (e.g., China) and western (e.g., India) parts of their geographic range have revealed major genetic differences that warrant the recognition of two different subspecies, little is known about genetic characteristics of rhesus macaques in the transitional zone extending from eastern India and Bangladesh through the northern part of Indo-China, the probable original homeland of the species. We analyzed genetic variation of 762 base pairs of mitochondrial DNA from 86 fecal swab samples and 19 blood samples from 25 local populations of rhesus macaque in Bangladesh collected from January 2010 to August 2012. These sequences were compared with those of rhesus macaques from India, China, and Myanmar. Forty-six haplotypes defined by 200 (26%) polymorphic nucleotide sites were detected. Estimates of gene diversity, expected heterozygosity, and nucleotide diversity for the total population were 0.9599 ± 0.0097, 0.0193 ± 0.0582, and 0.0196 ± 0.0098, respectively. A mismatch distribution of paired nucleotide differences yielded a statistically significantly negative value of Tajima's D, reflecting a population that rapidly expanded after the terminal Pleistocene. Most haplotypes throughout regions of Bangladesh, including an isolated region in the southwestern area (Sundarbans), clustered with haplotypes assigned to the minor haplogroup Ind-2 from India reflecting an east to west dispersal of rhesus macaques to India. Haplotypes from the southeast region of Bangladesh formed a cluster with those from Myanmar, and represent the oldest rhesus macaque haplotypes of Bangladesh. These results are consistent with the hypothesis that rhesus macaques first entered Bangladesh from the southeast, probably from Indo-China, then dispersed westward throughout eastern and central India.

Ecological genetics of Chinese rhesus macaque in response to mountain building: all things are not equal

BACKGROUND

Pliocene uplifting of the Qinghai-Tibetan Plateau (QTP) and Quaternary glaciation may have impacted the Asian biota more than any other events. Little is documented with respect to how the geological and climatological events influenced speciation as well as spatial and genetic structuring, especially in vertebrate endotherms. Macaca mulatta is the most widely distributed non-human primate. It may be the most suitable model to test hypotheses regarding the genetic consequences of orogenesis on an endotherm.

METHODOLOGY AND PRINCIPAL FINDINGS

Using a large dataset of maternally inherited mitochondrial DNA gene sequences and nuclear microsatellite DNA data, we discovered two maternal super-haplogroups exist, one in western China and the other in eastern China. M. mulatta formed around 2.31 Ma (1.51-3.15, 95%), and divergence of the two major matrilines was estimated at 1.15 Ma (0.78-1.55, 95%). The western super-haplogroup exhibits significant geographic structure. In contrast, the eastern super-haplogroup has far greater haplotypic variability with little structure based on analyses of six variable microsatellite loci using Structure and Geneland. Analysis using Migrate detected greater gene flow from WEST to EAST than vice versa. We did not detect signals of bottlenecking in most populations.

CONCLUSIONS

Analyses of the nuclear and mitochondrial datasets obtained large differences in genetic patterns for M. mulatta. The difference likely reflects inheritance mechanisms of the maternally inherited mtDNA genome versus nuclear biparentally inherited STRs and male-mediated gene flow. Dramatic environmental changes may be responsible for shaping the matrilineal history of macaques. The timing of events, the formation of M. mulatta, and the divergence of the super-haplogroups, corresponds to both the uplifting of the QTP and Quaternary climatic oscillations. Orogenesis likely drove divergence of western populations in China, and Pleistocene glaciations are likely responsible for genetic structuring in the eastern super-haplogroup via geographic isolation and secondary contact.

Extensive DRB region diversity in cynomolgus macaques: recombination as a driving force

The DR region of primate species is generally complex and displays diversity concerning the number and combination of distinct types of DRB genes present per region configuration. A highly variable short tandem repeat (STR) present in intron 2 of nearly all primate DRB genes can be utilized as a quick and accurate high through-put typing procedure. This approach resulted previously in the description of unique and haplotype-specific DRB-STR length patterns in humans, chimpanzees, and rhesus macaques. For the present study, a cohort of 230 cynomolgus monkeys, including self-sustaining breeding groups, has been examined. MtDNA analysis showed that most animals originated from the Indonesian islands, but some are derived from the mainland, south and north of the Isthmus of Kra. Haplotyping and subsequent sequencing resulted in the detection of 118 alleles, including 28 unreported ones. A total of 49 Mafa-DRB region configurations were detected, of which 28 have not yet been described. Humans and chimpanzees possess a low number of different DRB region configurations in concert with a high degree of allelic variation. In contrast, however, allelic heterogeneity within a given Mafa-DRB configuration is even less frequently observed than in rhesus macaques. Several of these region configurations appear to have been generated by recombination-like events, most probably propagated by a retroviral element mapping within DRB6 pseudogenes, which are present on the majority of haplotypes. This undocumented high level of DRB region configuration-associated diversity most likely represents a species-specific strategy to cope with various pathogens.

Comparative genetics of a highly divergent DRB microsatellite in different macaque species

The DRB region of the major histocompatibility complex (MHC) of cynomolgus and rhesus macaques is highly plastic, and extensive copy number variation together with allelic polymorphism makes it a challenging enterprise to design a typing protocol. All intact DRB genes in cynomolgus monkeys (Mafa) appear to possess a compound microsatellite, DRB-STR, in intron 2, which displays extensive length polymorphism. Therefore, this STR was studied in a large panel of animals, comprising pedigreed families as well. Sequencing analysis resulted in the detection of 60 Mafa-DRB exon 2 sequences that were unambiguously linked to the corresponding microsatellite. Its length is often allele specific and follows Mendelian segregation. In cynomolgus and rhesus macaques, the nucleotide composition of the DRB-STR is in concordance with the phylogeny of exon 2 sequences. As in humans and rhesus monkeys, this protocol detects specific combinations of different DRB-STR lengths that are unique for each haplotype. In the present panel, 22 Mafa-DRB region configurations could be defined, which exceeds the number detected in a comparable cohort of Indian rhesus macaques. The results suggest that, in cynomolgus monkeys, even more frequently than in rhesus macaques, new haplotypes are generated by recombination-like events. Although both macaque species are known to share several identical DRB exon 2 sequences, the lengths of the corresponding microsatellites often differ. Thus, this method allows not only fast and accurate DRB haplotyping but may also permit discrimination between highly related macaque species.

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

Genotypes for 13 short tandem repeats (STRs) were used to assess the genetic diversity within and differentiation among populations of rhesus macaques (Macaca mulatta) from mainland Asia and long-tailed macaques (M. fascicularis) from mainland and insular Southeast Asia. These animals were either recently captured in the wild or derived from wild-caught founders maintained in captivity for biomedical research.A large number of alleles is shared between the two macaque species but a significant genetic division between them persists. This distinction is more clear-cut among populations that are not, or are unlikely to have recently been, geographically contiguous. Our results suggest there has been significant interspecies nuclear gene flow between rhesus macaques and long-tailed macaques on the mainland. Comparisons of mainland and island populations of long-tailed macaques reflect marked genetic subdivisions due to barriers to migration. Geographic isolation has restricted gene flow, allowing island populations to become subdivided and genetically differentiated. Indonesian long-tailed macaques show evidence of long-term separation and genetic isolation from the mainland populations, while long-tailed macaques from the Philippines and Mauritius both display evidence of founder effects and subsequent isolation, with the impact from genetic drift being more profound in the latter.

Genetic characterization of wild and captive rhesus macaques in China

The genetic structures of wild and captive rhesus macaque populations within China were compared by analyzing the mtDNA sequences of 203 captive-bred Chinese rhesus macaques with 77 GenBank sequences from wild-caught animals trapped throughout China. The genotypes of 22 microsatellites of captive Chinese rhesus macaques were also compared with those of captive Indian animals. The Chinese population is significantly differentiated from the Indian population and is more heterogeneous. Thus, compared with Indian rhesus macaques the phenotypic variance of traits with high heritability will be inflated in Chinese animals. Our data suggest that the western Chinese provinces have more subdivided populations than the eastern and southern Chinese provinces. The southern Chinese populations are the least structured and might have been more recently established. Human-mediated interbreeding among captive Chinese populations has occurred, implying that Chinese breeding strategies can influence the interpretation of biomedical research in the USA.

Characterization of the major histocompatibility complex class II DQB (MhcMamu-DQB1) alleles in a cohort of Chinese rhesus macaques (Macaca mulatta)

Rhesus macaques have long been used in animal models for various human diseases, the susceptibility and/or resistance to some of which have been associated with the major histocompatibilty complex (MHC). To gain insight into the MHC background and to facilitate the experimental use of Chinese rhesus macaques, the second exon of MhcMamu-DQB1 genes in 105 rhesus macaques were characterized by cloning and sequencing. A total of 37 MhcMamu-DQB1 alleles were identified, illustrating a marked allelic polymorphism at DQB1 in these monkeys. In addition to 10 alleles were novel sequences that had not been documented in earlier reports, at least 14 alleles reported in earlier studies were not detected in this study. Most of the sequences (73%) observed in this study belong to DQB1 06 (13 alleles) and DQB1 18 (14 alleles) lineages, and the rest (27%) belong to DQB1 15, DQB1 16 and DQB1 17 lineages. The most frequent allele detected among these monkeys was MhcMamu-DQB1 06111 (22%), followed by DQB1 1503 (19%); and most of the novel alleles were present at a frequency of less than 2.5%. As for individual animals, 24 of 105 (23%) were homozygous whereas 81 of 105 (77%) were heterozygous at the MhcMamu-DQB1 locus. These data indicated significant differences in MhcMamu-DQB1 allele distribution between the Chinese rhesus macaques and the previously reported rhesus macaques, which were mostly of Indian origin. This information will not only promote the understanding of rhesus macaque MHC diversity and polymorphism but will also facilitate the use of Chinese rhesus macaques in human disease studies, especially those that may be associated with HLA-DQB genes.

Pyrosequencing as a method for SNP identification in the rhesus macaque (Macaca mulatta)

Background

Rhesus macaques (Macaca mulatta) are the primate most used for biomedical research, but phenotypic differences between Indian-origin and Chinese rhesus macaques have encouraged genetic methods for identifying genetic differences between these two populations. The completion of the rhesus genome has led to the identification of many single nucleotide polymorphisms (SNPs) in this species. These single nucleotide polymorphisms have many advantages over the short tandem repeat (STR) loci currently used to assay genetic variation. However, the number of currently identified polymorphisms is too small for whole genome analysis or studies of quantitative trait loci. To that end, we tested a combination of methods to identify large numbers of high-confidence SNPs, and screen those with high minor allele frequencies (MAF).

Results

By testing our previously reported single nucleotide polymorphisms, we identified a subset of high-confidence, high-MAF polymorphisms. Resequencing revealed a large number of regionally specific SNPs not identified through a single pyrosequencing run. By resequencing a pooled sample of four individuals, we reliably identified loci with a MAF of at least 12.5%. Finally, we found that when applied to a larger, geographically variable sample of rhesus, a large proportion of our loci were variable in both populations, and very few loci were ancestry informative. Despite this fact, the SNP loci were more effective at discriminating Indian and Chinese rhesus than STR loci.

Conclusion

Pyrosequencing and pooled resequencing are viable methods for the identification of high-MAF SNP loci in rhesus macaques. These SNP loci are appropriate for screening both the inter- and intra-population genetic variation.

Flow cytometric characterization of T lymphocyte subsets in the peripheral blood of Chinese rhesus macaques: normal range, age- and sex-related differences

Available data on the normal levels of white blood cell populations in healthy rhesus macaques (Macaca mulatta) originated and living in China is scanty. To obtain such data, blood samples from 150 Chinese rhesus macaques were collected and the normal range of white blood cells and their subsets were analyzed according to age and sex by flow cytometry. CBC data showed that the count of total white blood cells and lymphocytes decreased with age. Phenotypic analysis of CD4 and CD8 expression on CD3+ T lymphocytes showed that the percentage of CD4+ T cells (51.4+/-9.6%), CD4-CD8- T cells (8.5+/-4.1%) and the ratio of CD4+ T to CD8+ T cells (1.26+/-0.55) decreased with age; and the percentage of CD8+ T cells (42.0+/-9.7%), CD4+CD8+ T cells (1.3+/-0.9%) and CD3+ lymphocytes (55.3+/-13.3%) increased with age. However, no statistically significant difference was observed between the male and female groups in most parameters in these monkeys except for the percentage of CD4+CD8+ T cells. This study provided basic information about blood cell count and T lymphocyte subsets in Chinese rhesus macaques. It may be useful for comparative studies using Indian and Chinese rhesus macaques.

Mitochondrial DNA variation within and among regional populations of longtail macaques (Macaca fascicularis) in relation to other species of the fascicularis group of macaques

An 835 base pair (bp) fragment of mitochondrial DNA (mtDNA) was sequenced to characterize genetic variation within and among 1,053 samples comprising five regional populations each of longtail macaques (Macaca fascicularis) and rhesus macaques (Macaca mulatta), and one sample each of Japanese (M. fuscata) and Taiwanese (M. cyclopis) macaques. The mtDNA haplotypes of longtail macaques clustered in two large highly structured clades (Fas1 and Fas2) of a neighbor-joining tree that were reciprocally monophyletic with respect to those representing rhesus macaques, Japanese macaques, and Taiwanese macaques. Both clades exhibited haplotypes of Indonesian and Malaysian longtail macaques widely dispersed throughout them; however, longtail macaques from Indochina, Philippines, and Mauritius each clustered in a separate well-defined clade together with one or a few Malaysian and/or Indonesian longtail macaques, suggesting origins on the Sunda shelf. Longtail macaques from Malaysia and Indonesia were far more genetically diverse, and those from Mauritius were far less diverse than any other population studied. Nucleotide diversity between mtDNA sequences of longtail macaques from different geographic regions is, in some cases, greater than that between Indian and Chinese rhesus macaques. Approximately equal amounts of genetic diversity are due to differences among animals in the same regional population, different regional populations, and different species. A greater proportion of genetic variance was explained by interspecies differences when Japanese and Taiwanese macaques were regarded as regional populations of rhesus macaques than when they were treated as separate species. Rhesus macaques from China were more closely related to both Taiwanese and Japanese macaques than to their own conspecifics from India.

Single nucleotide polymorphisms (SNPs) distinguish Indian-origin and Chinese-origin rhesus macaques (Macaca mulatta)

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.

Considerations in the selection and conditioning of Old World monkeys for laboratory research: animals from domestic sources

Nonhuman primates from domestic sources constitute an important resource for the research community. The life history of the Old World monkey species that comprise the bulk of this resource is described, and issues that colony managers and researchers alike should consider regarding animal selection (e.g., species, age, sex, rearing history, temperament, genotype, viral status, geographic origin) are discussed. Preparation of domestically bred animals for research usually involves some combination of social separation, relocation, resocialization, alterations in physical space, photoperiod, and diet, as well as exposure to novel environments. The research literature that has focused on these issues is reviewed, and authors suggest that once animals have been assigned to their project housing situation, a period ranging up to 3 mo (depending on the magnitude of the change in housing) might be warranted before an experimental protocol should begin. Attention to issues of animal selection and conditioning by both researchers and colony managers can lead to the shared goal of high-quality research that utilizes the minimal number of animals.