Neural signatures of natural behaviour in socializing macaques.

Our understanding of the neurobiology of primate behaviour largely derives from artificial tasks in highly controlled laboratory settings, overlooking most natural behaviours that primate brains evolved to produce1-3. How primates navigate the multidimensional social relationships that structure daily life4 and shape survival and reproductive success5 remains largely unclear at the single-neuron level. Here we combine ethological analysis, computer vision and wireless recording technologies to identify neural signatures of natural behaviour in unrestrained, socially interacting pairs of rhesus macaques. Single-neuron and population activity in the prefrontal and temporal cortex robustly encoded 24 species-typical behaviours, as well as social context. Male-female partners demonstrated near-perfect reciprocity in grooming, a key behavioural mechanism supporting friendships and alliances6, and neural activity maintained a running account of these social investments. Confronted with an aggressive intruder, behavioural and neural population responses reflected empathy and were buffered by the presence of a partner. Our findings reveal a highly distributed neurophysiological ledger of social dynamics, a potential computational foundation supporting communal life in primate societies, including our own.

Copy number variants and fixed duplications among 198 rhesus macaques (Macaca mulatta).

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

Long-read assembly of the Chinese rhesus macaque genome and identification of ape-specific structural variants.

We present a high-quality de novo genome assembly (rheMacS) of the Chinese rhesus macaque (Macaca mulatta) using long-read sequencing and multiplatform scaffolding approaches. Compared to the current Indian rhesus macaque reference genome (rheMac8), rheMacS increases sequence contiguity 75-fold, closing 21,940 of the remaining assembly gaps (60.8 Mbp). We improve gene annotation by generating more than two million full-length transcripts from ten different tissues by long-read RNA sequencing. We sequence resolve 53,916 structural variants (96% novel) and identify 17,000 ape-specific structural variants (ASSVs) based on comparison to ape genomes. Many ASSVs map within ChIP-seq predicted enhancer regions where apes and macaque show diverged enhancer activity and gene expression. We further characterize a subset that may contribute to ape- or great-ape-specific phenotypic traits, including taillessness, brain volume expansion, improved manual dexterity, and large body size. The rheMacS genome assembly serves as an ideal reference for future biomedical and evolutionary studies.

Genetic contributions to dental dimensions in brown-mantled tamarins (Saguinus fuscicollis) and rhesus macaques (Macaca mulatta).

OBJECTIVES: The use of dental metrics in phylogenetic reconstructions of fossil primates assumes variation in tooth size is highly heritable. Quantitative genetic studies in humans and baboons have estimated high heritabilities for dental traits, providing a preliminary view of the variability of dental trait heritability in nonhuman primate species. To expand upon this view, the heritabilities and evolvabilities of linear dental dimensions are estimated in brown-mantled tamarins (Saguinus fuscicollis) and rhesus macaques (Macaca mulatta). MATERIALS AND METHODS: Quantitative genetic analyses were performed on linear dental dimensions collected from 302 brown-mantled tamarins and 364 rhesus macaques. Heritabilities were estimated in SOLAR using pedigrees from each population, and evolvabilities were calculated manually. RESULTS: Tamarin heritability estimates range from 0.19 to 0.99, and 25 of 26 tamarin estimates are significantly different from zero. Macaque heritability estimates range from 0.08 to 1.00, and 25 out of 28 estimates are significantly different from zero. DISCUSSION: Dental dimensions are highly heritable in captive brown-mantled tamarins and free-ranging rhesus macaques. The range of heritability estimates in these populations is broadly similar to those of baboons and humans. Evolvability tends to increase with heritability, although evolvability is high relative to heritability in some dimensions. Estimating evolvability helps to contextualize differences in heritability, and the observed relationship between evolvability and heritability in dental dimensions requires further investigation.

Comparison of cellular localisation of the Ca2+ -binding proteins calbindin, calretinin and parvalbumin in the retina of four different Macaca species.

Ca2+ -binding proteins are differentially expressed in the nervous system; their functional role often remains unclear. This immunohistochemical study aimed at characterising and comparing the expression pattern of the Ca2+ -binding proteins calbindin (Calb), calretinin (Calr) and parvalbumin (Parv) in the retina of four species of macaque monkeys: Macaca fascicularis (cynomolgus macaque), M. mulatta (rhesus macaque), M. thibetana (Tibetan macaque) and M. fuscata (Japanese macaque). Calb was found in cone photoreceptors and in a subset of bipolar cells. Calr was expressed in a subpopulation of amacrine cells. Parv was present in horizontal and ganglion cells. In addition, Müller cells were stained using antibodies against the specific marker cellular retinaldehyde-binding protein (CRALBP). Immunostainings were used for calculation of the density of different cell populations. The expression pattern was similar between the examined species and between retinal regions.

A cerebellar substrate for cognition evolved multiple times independently in mammals.

Given that complex behavior evolved multiple times independently in different lineages, a crucial question is whether these independent evolutionary events coincided with modifications to common neural systems. To test this question in mammals, we investigate the lateral cerebellum, a neurobiological system that is novel to mammals, and is associated with higher cognitive functions. We map the evolutionary diversification of the mammalian cerebellum and find that relative volumetric changes of the lateral cerebellar hemispheres (independent of cerebellar size) are correlated with measures of domain-general cognition in primates, and are characterized by a combination of parallel and convergent shifts towards similar levels of expansion in distantly related mammalian lineages. Results suggest that multiple independent evolutionary occurrences of increased behavioral complexity in mammals may at least partly be explained by selection on a common neural system, the cerebellum, which may have been subject to multiple independent neurodevelopmental remodeling events during mammalian evolution.

Phylogenetic relationship of a fossil macaque (Macaca cf. robusta) from the Korean Peninsula to extant species of macaques based on zygomaxillary morphology.

Little is known about the biogeographical and evolutionary histories of macaques (Macaca spp.) in East Asia because the phylogenetic positions of fossil species remain unclear. Here we examined the zygomaxillary remains of a fossil macaque (M. cf. robusta) from the Durubong Cave Complex, South Korea, that dates back to the late Middle to Late Pleistocene, to infer its phylogenetic relationship to extant species. We took 195 fixed- and semi-landmarks from the zygomaxillary regions of the fossil specimen and from 147 specimens belonging to 14 extant species. We then conducted a generalized Procrustes analysis followed by a multivariate statistical analysis to evaluate the phenetic affinities of the fossil to the extant species and reconstructed the most parsimonious phylogenetic tree using a phylogenetic morphometric approach. We found that the fossil was most similar to Macaca fuscata (Japanese macaque) in the zygomaxillary morphospace although it was at the limit of the range of variation for this species. The second closest in the morphospace was the continental Macaca mulatta (rhesus macaque). Parsimonious reconstruction confirmed that the fossil was most closely related to M. fuscata, even after controlling for the effects of allometry. These findings suggest that in the late Middle to Late Pleistocene, close relatives of M. fuscata that looked like the extant species were distributed on the Korean Peninsula, where no species of macaques are found today. Thus, some morphological characteristics of M. fuscata may have developed before its ancestor dispersed into the Japanese archipelago.

Predominant patterns of splicing evolution on human, chimpanzee and macaque evolutionary lineages.

Although splicing is widespread and evolves rapidly among species, the mechanisms driving this evolution, as well as its functional implications, are not yet fully understood. We analyzed the evolution of splicing patterns based on transcriptome data from five tissues of humans, chimpanzees, rhesus macaques and mice. In total, 1526 exons and exon sets from 1236 genes showed significant splicing differences among primates. More than 60% of these differences represent constitutive-to-alternative exon transitions while an additional 25% represent changes in exon inclusion frequency. These two dominant evolutionary patterns have contrasting conservation, regulation and functional features. The sum of these features indicates that, despite their prevalence, constitutive-to-alternative exon transitions do not substantially contribute to long-term functional transcriptome changes. Conversely, changes in exon inclusion frequency appear to be functionally relevant, especially for changes taking place in the brain on the human evolutionary lineage.

The orthologs of HLA-DQ and -DP genes display abundant levels of variability in macaque species.

The human major histocompatibility complex (MHC) region encodes three types of class II molecules designated HLA-DR, -DQ, and -DP. Both the HLA-DQ and -DP gene region comprise a duplicated tandem of A and B genes, whereas in macaques, only one set of genes is present per region. A substantial sequencing project on the DQ and DP genes in various macaque populations resulted in the detection of previously 304 unreported full-length alleles. Phylogenetic studies showed that humans and macaques share trans-species lineages for the DQA1 and DQB1 genes, whereas the DPA1 and DPB1 lineages in macaques appear to be species-specific. Amino acid variability plot analyses revealed that each of the four genes displays more allelic variation in macaques than is encountered in humans. Moreover, the numbers of different amino acids at certain positions in the encoded proteins are higher than in humans. This phenomenon is remarkably prominent at the contact positions of the peptide-binding sites of the deduced macaque DPß-chains. These differences in the MHC class II DP regions of macaques and humans suggest separate evolutionary mechanisms in the generation of diversity.

The population genomics of rhesus macaques (Macaca mulatta) based on whole-genome sequences.

Rhesus macaques (Macaca mulatta) are the most widely used nonhuman primate in biomedical research, have the largest natural geographic distribution of any nonhuman primate, and have been the focus of much evolutionary and behavioral investigation. Consequently, rhesus macaques are one of the most thoroughly studied nonhuman primate species. However, little is known about genome-wide genetic variation in this species. A detailed understanding of extant genomic variation among rhesus macaques has implications for the use of this species as a model for studies of human health and disease, as well as for evolutionary population genomics. Whole-genome sequencing analysis of 133 rhesus macaques revealed more than 43.7 million single-nucleotide variants, including thousands predicted to alter protein sequences, transcript splicing, and transcription factor binding sites. Rhesus macaques exhibit 2.5-fold higher overall nucleotide diversity and slightly elevated putative functional variation compared with humans. This functional variation in macaques provides opportunities for analyses of coding and noncoding variation, and its cellular consequences. Despite modestly higher levels of nonsynonymous variation in the macaques, the estimated distribution of fitness effects and the ratio of nonsynonymous to synonymous variants suggest that purifying selection has had stronger effects in rhesus macaques than in humans. Demographic reconstructions indicate this species has experienced a consistently large but fluctuating population size. Overall, the results presented here provide new insights into the population genomics of nonhuman primates and expand genomic information directly relevant to primate models of human disease.

The Population Genetic Composition of Conventional and SPF Colonies of Rhesus Macaques (Macaca mulatta) at the Caribbean Primate Research Center.

The SPF breeding program at the Caribbean Primate Research Center supplies Indian-origin rhesus macaques of known genetic and virologic background for biomedical research. In this study, population genetic analyses using 14 short tandem-repeat sequences showed that the SPF colony has remained genetically homogenous over time, with sufficient amounts of heterozygosity and minimal stratification from its founders. Intergenerational studies indicated that an average of 7 alleles have been retained, inbreeding levels have remained low, and the degree of Indian ancestry is one of the highest among several national primate research centers. The relative low genetic diversity in the free-ranging population as well as in the captive SPF and conventional colonies when compared with that of other primate centers indicates that the free-ranging population, from which the captive-colony animals were derived, has experienced significant founder effects and genetic drift during the years after its establishment. This study supports the historical origin of the free-ranging population and confirms the high value of this resource for biomedical research. Current genetic diversity levels within the SPF colony can be ensured with the practice of colony management approaches such as equalizing male:female ratios in each SPF breeding group and increasing breeding group sizes. Introducing new Indian-origin macaques from other captive colonies might help to maximize the genetic diversity of the breeding stock. Furthermore, genetic estimates must be used to rank breeders according to their genetic value or their genome uniqueness to increase founder-genome representation and curb future genetic bottlenecks and allele loss.

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.

Divergence and diversity of ULBP2 genes in rhesus and cynomolgus macaques.

Non-human primates such as rhesus macaque and cynomolgus macaque are important animals for medical research fields and they are classified as Old World monkey, in which genome structure is characterized by gene duplications. In the present study, we investigated polymorphisms in two genes for ULBP2 molecules that are ligands for NKG2D. A total of 15 and 11 ULBP2.1 alleles and 11 and 10 ULBP2.2 alleles were identified in rhesus macaques and cynomolgus macaques, respectively. Nucleotide sequences of exons for extra cellular domain were highly polymorphic and more than 70 % were non-synonymous variations in both ULBP2.1 and ULBP2.2. In addition, phylogenetic analyses revealed that the ULBP2.2 was diverged from a branch of ULBP2.1 along with ULBP2s of higher primates. Moreover, when 3D structural models were constructed for the rhesus ULBP2 molecules, residues at presumed contact sites with NKG2D were polymorphic in ULBP2.1 and ULBP2.2 in the rhesus macaque and cynomolgus macaque, respectively. These observations suggest that amino acid replacements at the interaction sites with NKG2D might shape a specific nature of ULBP2 molecules in the Old World monkeys.

Mamu-B genes and their allelic repertoires in different populations of Chinese-origin rhesus macaques.

Since rhesus monkeys of Chinese origin have gained greater utilization in recent years, it is urgent to investigate the major histocompatibility complex (MHC) immunogenetics of Chinese rhesus macaques. In this study, we identified 81 Mamu-B sequences using complementary DNA cloning and sequencing on a cohort of 58 rhesus monkeys derived from three local populations of China. Twenty of these Mamu-B alleles are novel and four of them represent new lineages. Although more alleles are shared among different populations than Mamu-A locus, the Mamu-B allelic repertoires found in these three populations of Chinese macaques are largely independent, which underscores the MHC polymorphism among different populations of Chinese rhesus macaques. Our results are an important addition to the limited MHC immunogenetic information available for rhesus macaques of Chinese origin.

Monkey models for brain-machine interfaces: the need for maintaining diversity.

Brain-machine interfaces (BMIs) aim to help disabled patients by translating neural signals from the brain into control signals for guiding prosthetic arms, computer cursors, and other assistive devices. Animal models are central to the development of these systems and have helped enable the successful translation of the first generation of BMIs. As we move toward next-generation systems, we face the question of which animal models will aid broader patient populations and achieve even higher performance, robustness, and functionality. We review here four general types of rhesus monkey models employed in BMI research, and describe two additional, complementary models. Given the physiological diversity of neurological injury and disease, we suggest a need to maintain the current diversity of animal models and to explore additional alternatives, as each mimic different aspects of injury or disease.

Bioinformatic approaches to identifying orthologs and assessing evolutionary relationships.

Non-human primate genetic research defines itself through comparisons to humans; few other species require the implicit comparative genomics approaches. Because of this, errors in the identification of non-human primate orthologs can have profound effects. Gene prediction algorithms can and have produced false transcripts that have become incorporated into commonly used databases and genomics portals. These false transcripts can arise from deficiencies in the algorithms themselves as well as through gaps and other problems in the genome assembly. Putative genes generated can not only miss microexons, but improperly incorporate non-coding sequence resulting in pseudogenes or other transcripts without biological relevance. False transcripts then become identified as orthologs to established human genes and are too often taken as gospel by unwary researchers. Here, the processes through which these errors propagate are isolated and methods are described for identifying false orthologs in databases with several representative errors illustrated. Through these steps any researcher seeking to make use of non-human primate genetic information will have the tools at their disposal to ascertain where errors exist and to remedy them once encountered.

The efficacy of diazepam treatment for the management of acute wounding episodes in captive rhesus macaques.

The spontaneous development of self-injurious behavior (SIB) in singly housed monkeys poses a challenge for their management and well-being in captivity. Relatively little information is available on effective treatments for SIB. This study examined the effects of diazepam (Valium) on self-wounding and other abnormal behaviors in eight individually housed male rhesus monkeys (Macaca mulatta). Each monkey's response to an anxiolytic dose of diazepam (1 mg/kg or greater orally) was compared with the animal's behavior during drug-free periods. When examined across all animals, treatment with diazepam did not significantly alter wounding frequency or rates of self-directed biting without wounding. However, closer examination of the data revealed that four of the animals showed significant decreases in self-biting and wounding frequency (positive responders, PR group), whereas the remaining monkeys showed a trend towards increased wounding frequency (negative responders, NR group). Subsequent examination of colony and veterinary records demonstrated that compared with NR monkeys, PR monkeys had spent significantly more years in individual cage housing and had experienced a greater number of minor veterinary procedures. PR animals also were significantly less likely to have a documented history of self-biting behavior. Our findings suggest that SIB is not a homogeneous disorder in rhesus monkeys; rather, distinct subtypes exist that require different treatment approaches.

Genetic characterization of Indian-origin and Chinese-origin rhesus macaques (Macaca mulatta).

Genetic differences between Indian-origin and Chinese-origin rhesus macaques are as great as those between some primate species and can influence the results of experiments in which both are used as animal models for the study of the same human diseases. Unfortunately, many breeding facilities do not know with certainty the origin of the founders of their rhesus breeding colonies. Here I summarize the most definitive of the genetic traits among the microsatellite (STR) loci and mitochondrial DNA sequences that my laboratory previously reported to characterize Indian-origin and Chinese-origin rhesus macaques and then estimate the frequencies of these traits and their reliability as indicators of country of origin. The expression of diagnostic traits at two or more of four different unlinked loci provides a nearly 100% reliability in distinguishing rhesus macaques of Indian and Chinese origin.

Ontogeny of limb mass distribution in infant baboons (Papio cynocephalus).

Primates have more distally distributed limb muscle mass compared to most nonprimate mammals. The heavy distal limbs of primates are likely related to their strong manual and pedal grasping abilities, and interspecific differences in limb mass distributions among primates are correlated with the amount of time spent on arboreal supports. Within primate species, individuals at different developmental stages appear to differ in limb mass distribution patterns. For example infant macaques have more distally distributed limb mass at young ages. A shift from distal to proximal limb mass concentrations coincides with a shift from dependent travel (grasping their mother's hair) to independent locomotion. Because the functional demands placed on limbs may differ between taxa, understanding the ontogeny of limb mass distribution patterns is likely an essential element in interpreting the diversity of limb mass distribution patterns present in adult primates. This study examines changes in limb inertial properties during ontogeny in a longitudinal sample of infant baboons (Papio cynocephalus). The results of this study show that infant baboons undergo a transition from distal to proximal limb mass distribution patterns. This transition in limb mass distribution coincides with the transition from dependent to independent locomotion during infant development. Compared to more arboreal macaques, infant baboons undergo a faster transition to more proximal limb mass distribution patterns. These results suggest that functional demands placed on the limbs during ontogeny have a strong impact on the development of limb mass distribution patterns.