Whole-genome sequencing of six Mauritian Cynomolgus macaques (Macaca fascicularis) reveals a genome-wide pattern of polymorphisms under extreme population bottleneck

Genomic Evidence for the Complex Evolutionary History of Macaques (Genus Macaca)

The genus Macaca is widely distributed, occupies a variety of habitats, shows diverse phenotypic characteristics, and is one of the best-studied genera of nonhuman primates. Here, we reported five re-sequencing Macaca genomes, including one M. cyclopis, one M. fuscata, one M. thibetana, one M. silenus, and one M. sylvanus. Together with published genomes of other macaque species, we combined 20 genome sequences of 10 macaque species to investigate the gene introgression and genetic differences among the species. The network analysis of the SNV-fragment trees indicates a reticular phylogeny of macaque species. Combining the results from various analytical methods, we identified extensive ancient introgression events among macaque species. The multiple introgression signals between different species groups were also observed, such as between fascicularis group species and silenus group species. However, gene flow signals between fascicularis and sinica group were not as strong as those between fascicularis group and silenus group. On the other hand, the unidirect gene flow in M. arctoides probably occurred between the progenitor of M. arctoides and the common ancestor of fascicularis group. Our study also shows that the genetic backgrounds and genetic diversity of different macaques vary dramatically among species, even among populations of the same species. In conclusion, using whole genome sequences and multiple methods, we have studied the evolutionary history of the genus Macaca and provided evidence for extensive introgression among the species.

The phylogenetic relationship and demographic history of rhesus macaques (Macaca mulatta) in subtropical and temperate regions, China

Pleistocene climatic oscillations exerted significant influences on the genetic structure and demography of rhesus macaque (Macaca mulatta) in eastern China. However, the evolutionary history of rhesus macaques in subtropical and temperate China remained unclear and/or controversial. Herein, we analyzed the autosomes, mitochondrial genomes, and Y-chromosomes from 84 individuals of Chinese rhesus macaque. The results revealed that (1) all individuals were clustered into pan-west and pan-east genetic groups, which exhibited Shaanxi Province as the northernmost region of western dispersal route of rhesus macaques in China; (2) in subtropical and temperate China, rhesus macaques were divided into four lineages (TH, DB, HS, and QL), and their divergence times corresponded to the Penultimate Glaciation (300-130 kya) and Last Glaciation (70-10 kya), respectively; (3) the individuals from Mt. Taihangshan (TH) are closely related to individuals from Mt. Dabashan (DB) in the autosomal tree, rather than individuals from Mt. Huangshan (HS) as indicated by the mitogenome tree, which supports the hypothesis that the ancestral rhesus macaques radiated into Mt. Taihangshan from Mt. Huangshan via Mt. Dabashan; and (4) the demographic scenario of the four lineages showed the ancestral rhesus macaques bottleneck and expansion corresponding to the suitable habitat reduction and expansion, which confirmed they had experienced northward recolonization and southward retreat events from Mt. Huangshan area via Northern China Plain to Northernmost China along with Pleistocene glacial cycles. This study provides a new insight into understanding how Pleistocene glaciation has influenced faunal diversity in subtropical and temperate China, especially for those exhibiting differential patterns of sex dispersal.

Phylogeography and demographic history of macaques, fascicularis species group, in East Asia: Inferred from multiple genomic markers

Climate changes at larger scales have influenced dispersal and range shifts of many taxa in East Asia. The fascicularis species group of macaques is composed of four species and is widely distributed in Southeast and East Asia. However, its phylogeography and demographic histories are currently poorly understood. Herein, we assembled autosomal, mitogenome, and Y-chromosome data for 106 individuals, and combined them with 174 mtDNA dloop haplotypes of this species group, with particular focus on the demographic histories and dispersal routes of Macaca fuscata, M. cyclopis, and M. mulatta. The results showed: (1) three monophyletic clades for M. fuscata, M. cyclopis, and M. mulatta based on the multiple genomics analyses; (2) the disparate demographic trajectories of the three species after their split ∼1.0 Ma revealed that M. cyclopis and M. fuscata were derived from an ancestral M. mulatta population; (3) the speciation time of M. cyclopis was later than that of M. fuscata, and their divergence time occurred at the beginning of "Ryukyu Coral Sea Stage" (1.0-0.2 Ma) when the East China Sea land bridge was completely submerged by the sea level rose; and (4) the three parallel rivers (Nujiang, Lancangjiang, and Jinshajiang) of Southwestern China divided M. mulatta into Indian and Chinese genetic populations ∼200 kya. These results shed light on understanding not only the evolutionary history of the fascicularis species group but also the formation mechanism of faunal diversity in East Asia during the Pleistocene.

DNA methylation and transcriptome analysis reveal epigenomic differences among three macaque species

Macaques (genus Macaca) are the most widely distributed non-human primates, and their evolutionary history, gene expression profiles, and genetic differences have been extensively studied. However, the DNA methylomes of macaque species are not available in public databases, which hampers understanding of epigenetic differences among macaque species. Epigenetic modifications can potentially affect development, physiology, behavior, and evolution. Here, we investigated the methylation patterns of the Tibetan macaque (M. thibetana; TM), Chinese rhesus macaque (M. mulatta lasiota; CR), and crab-eating macaque (M. fascicularis; CE) through whole-genome bisulfite sequencing from peripheral blood. We compared genome-wide methylation site information for the three species. We identified 12,128 (CR vs. CE), 59,165 (CR vs. TM), and 39,751 (CE vs. TM) differentially methylated regions (DMRs) in the three macaques. Furthermore, we obtained the differentially expressed genes (DEGs) among the three macaque species. The differences between CR and CE were smaller at both the methylome and transcriptome levels than compared with TM (CR vs. TM and CE vs. TM). We also found a change in the density of single nucleotide mutations in DMRs relative to their flanking regions, indicating a potential mechanism through which genomic alterations may modulate methylation landscapes, thereby influencing the transcriptome. Functional enrichment analyses showed the DMR-related genes were enriched in developmental processes and neurological functions, such as the growth hormone-related pathway, insulin secretion pathway, thyroid hormone synthesis pathway, morphine addiction, and GABAergic synapses. These differences may be associated with variations in physiology and habitat among the macaques. Our study provides one of the first genome-wide comparisons of genetic, gene expression, and epigenetic variations across different macaques. Our results should facilitate further research on comparative genomic and genetic differences in macaque species.

Special Issue on the Pathobiology of Laboratory Nonhuman Primates: A Review of Species, Substrain, Geographical Origin, Age, and Modality-Related Factors

Nonhuman primates (NHPs) are utilized in nonclinical safety testing due to their phylogenetic proximity to humans and similarity in physiology and anatomy. However, ethical considerations and the increased demand for NHPs, coupled with the current shortage in their supply, have increased the calls to minimize their use. In addition, the increased demand and supply shortage of NHPs have increased the use of animals sourced from different geographical origins, and animals of different ages, which can complicate the interpretation of study results. Coupled with the relative uniqueness of findings induced by novel therapeutic modalities, there is an increasing need for a deeper understanding of the systemic pathobiology of NHPs. Here we provide a brief preview of the two main themes discussed in this special issue, which include the influence of geographical origin, age, and sex on background pathology, clinical pathology reference values, other relevant toxicology endpoints, and organ system pathology.

The Influence of Geographical Origin, Age, Sex, and Animal Husbandry on the Spontaneous Histopathology of Laboratory Cynomolgus Macaques (Macaca Fascicularis): A Contemporary Global and Multisite Review of Historical Control Data

To investigate the influence of geographical origin, age, and sex on toxicologically relevant spontaneous histopathology findings in cynomolgus macaques (Macaca fascicularis), we performed a comparative analysis of historical control data (HCD) from 13 test sites that included 3351 animals (1645 females and 1706 males) sourced from Mauritius, China, Vietnam, and Cambodia, aged from 2 to 9.5 years, and from 446 toxicology studies evaluated between 2016 and 2021. The most common findings were mononuclear infiltrates in the kidney, liver, brain, and lung, which showed highest incidences in Mauritian macaques, and heart, salivary glands, and gastrointestinal tract (GIT), which showed highest incidences of mononuclear infiltrates in mainland Asian macaques. Developmental and degenerative findings were more common in Mauritian macaques, while lymphoid hyperplasia and lung pigment showed higher incidences in Asian macaques. Various sex and age-related differences were also present. Despite origin-related differences, the similarities in the nature and distribution of background lesions indicate that macaques from all geographical regions are suitable for toxicity testing and show comparable lesion spectrum. However, in a toxicity study, it is strongly recommended to use animals from a single geographical origin and to follow published guidelines when using HCD to evaluate and interpretate commonly diagnosed spontaneous lesions.

VDJ Gene Usage in IgM Repertoires of Rhesus and Cynomolgus Macaques

Macaques are frequently used to evaluate candidate vaccines and to study infection-induced antibody responses, requiring an improved understanding of their naïve immunoglobulin (IG) repertoires. Baseline gene usage frequencies contextualize studies of antigen-specific immune responses, providing information about how easily one may stimulate a response with a particular VDJ recombination. Studies of human IgM repertoires have shown that IG VDJ gene frequencies vary several orders of magnitude between the most and least utilized genes in a manner that is consistent across many individuals but to date similar analyses are lacking for macaque IgM repertoires. Here, we quantified VDJ gene usage levels in unmutated IgM repertoires of 45 macaques, belonging to two species and four commonly used subgroups: Indian and Chinese origin rhesus macaques and Indonesian and Mauritian origin cynomolgus macaques. We show that VDJ gene frequencies differed greatly between the most and least used genes, with similar overall patterns observed in macaque subgroups and individuals. However, there were also clear differences affecting the use of specific V, D and J genes. Furthermore, in contrast to humans, macaques of both species utilized IGHV4 family genes to a much higher extent and showed evidence of evolutionary expansion of genes of this family. Finally, we used the results to inform the analysis of a broadly neutralizing HIV-1 antibody elicited in SHIV-infected rhesus macaques, RHA1.V2.01, which binds the apex of the Env trimer in a manner that mimics the binding mode of PGT145. We discuss the likelihood that similar antibodies could be elicited in different macaque subgroups.

Demographic History of the Brown Bear (Ursus arctos) on Hokkaido Island, Japan, Based on Whole-Genomic Sequence Analysis

Previous studies of the brown bear (Ursus arctos) on Hokkaido Island, Japan, have detected three geographically distinct subpopulations representing different mitochondrial lineages and shown that gene flow between subpopulations has occurred due to male-biased dispersal. In this study, we determined whole-genomic sequences for six Hokkaido brown bears and analyzed these data along with previously published genomic sequences of 17 brown bears from other parts of the world. We found that the Hokkaido population is genetically distinct from the other populations, keeping genetic diversity higher than the endangered populations in western Europe but lower than most populations on the continents. A reconstruction of historical demography showed no increase in population size for the Hokkaido population during the Eemian interglacial period (130,000–114,000 years ago). In a phylogenetic analysis of the autosomal data, the Hokkaido population formed a clade distinct from North American and European populations, showing that it has maintained genetic diversity independently from continental populations following geographical isolation on the island. This autosomal genetic similarity contrasts with the geographically separate mitochondrial lineages on Hokkaido and indicates the occurrence of male-driven gene flow between subpopulations.

Population genomics reveal rapid genetic differentiation in a recently invasive population of Rattus norvegicus

Background

Invasive species bring a serious effect on local biodiversity, ecosystems, and even human health and safety. Although the genetic signatures of historical range expansions have been explored in an array of species, the genetic consequences of contemporary range expansions have received little attention, especially in mammal species. In this study, we used whole-genome sequencing to explore the rapid genetic change and introduction history of a newly invasive brown rat (Rattus norvegicus) population which invaded Xinjiang Province, China in the late 1970s.

Results

Bayesian clustering analysis, principal components analysis, and phylogenetic analysis all showed clear genetic differentiation between newly introduced and native rat populations. Reduced genetic diversity and high linkage disequilibrium suggested a severe population bottleneck in this colonization event. Results of TreeMix analyses revealed that the introduced rats were derived from an adjacent population in geographic region (Northwest China). Demographic analysis indicated that a severe bottleneck occurred in XJ population after the split off from the source population, and the divergence of XJ population might have started before the invasion of XJ. Moreover, we detected 42 protein-coding genes with allele frequency shifts throughout the genome for XJ rats and they were mainly associated with lipid metabolism and immunity, which could be seen as a prelude to future selection analyses in the novel environment of XJ.

Conclusions

This study presents the first genomic evidence on genetic differentiation which developed rapidly, and deepens the understanding of invasion history and evolutionary processes of this newly introduced rat population. This would add to our understanding of how invasive species become established and aid strategies aimed at the management of this notorious pest that have spread around the world with humans.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12983-021-00387-z.

Rhesus and cynomolgus macaque immunoglobulin heavy-chain genotyping yields comprehensive databases of germline VDJ alleles

Definition of the specific germline immunoglobulin (Ig) alleles present in an individual is a critical first step to delineate the ontogeny and evolution of antigen-specific antibody responses. Rhesus and cynomolgus macaques are important animal models for pre-clinical studies, with four main sub-groups being used: Indian- and Chinese-origin rhesus macaques and Mauritian and Indonesian cynomolgus macaques. We applied the (Ig) gene inference tool IgDiscover and performed extensive Sanger sequencing-based genomic validation to define germline VDJ alleles in these 4 sub-groups, comprising 45 macaques in total. There was allelic overlap between Chinese- and Indian-origin rhesus macaques and also between the two macaque species, which is consistent with substantial admixture. The island-restricted Mauritian cynomolgus population displayed the lowest number of alleles of the sub-groups, yet maintained high individual allelic diversity. These comprehensive databases of germline IGH alleles for rhesus and cynomolgus macaques provide a resource toward the study of B cell responses in these important pre-clinical models.

Identifying a Minor Histocompatibility Antigen in Mauritian Cynomolgus Macaques Encoded by APOBEC3C

Allogeneic hematopoietic stem cell transplants can lead to dramatic reductions in human immunodeficiency virus (HIV) reservoirs. This effect is partially mediated by donor T cells recognizing lymphocyte-expressed minor histocompatibility antigens (mHAgs). The potential to mark malignant and latently infected cells for destruction makes mHAgs attractive targets for cellular immunotherapies. However, testing such HIV reservoir reduction strategies will likely require preclinical studies in non-human primates (NHPs). In this study, we used a combination of alloimmunization, whole exome sequencing, and bioinformatics to identify an mHAg in Mauritian cynomolgus macaques (MCMs). We mapped the minimal optimal epitope to a 10-mer peptide (SW10) in apolipoprotein B mRNA editing enzyme catalytic polypeptide-like 3C (APOBEC3C) and determined the major histocompatibility complex class I restriction element as Mafa-A1∗063, which is expressed in almost 90% of MCMs. APOBEC3C SW10-specific CD8+ T cells recognized immortalized B cells but not fibroblasts from an mHAg-positive MCM. These results provide a framework for identifying mHAgs in a non-transplant setting and suggest that APOBEC3C SW10 could be used as a model antigen to test mHAg-targeted therapies in NHPs.

MACSNVdb: a high-quality SNV database for interspecies genetic divergence investigation among macaques

Macaques are the most widely used non-human primates in biomedical research. The genetic divergence between these animal models is responsible for their phenotypic differences in response to certain diseases. However, the macaque single nucleotide polymorphism resources mainly focused on rhesus macaque (Macaca mulatta), which hinders the broad research and biomedical application of other macaques. In order to overcome these limitations, we constructed a database named MACSNVdb that focuses on the interspecies genetic diversity among macaque genomes. MACSNVdb is a web-enabled database comprising ~74.51 million high-quality non-redundant single nucleotide variants (SNVs) identified among 20 macaque individuals from six species groups (muttla, fascicularis, sinica, arctoides, silenus, sylvanus). In addition to individual SNVs, MACSNVdb also allows users to browse and retrieve groups of user-defined SNVs. In particular, users can retrieve non-synonymous SNVs that may have deleterious effects on protein structure or function within macaque orthologs of human disease and drug-target genes. Besides position, alleles and flanking sequences, MACSNVdb integrated additional genomic information including SNV annotations and gene functional annotations. MACSNVdb will facilitate biomedical researchers to discover molecular mechanisms of diverse responses to diseases as well as primatologist to perform population genetic studies. We will continue updating MACSNVdb with newly available sequencing data and annotation to keep the resource up to date.

Database URL: http://big.cdu.edu.cn/macsnvdb/

Resources for functional genomic studies of health and development in nonhuman primates

Primates display a wide range of phenotypic variation underlaid by complex genetically regulated mechanisms. The links among DNA sequence, gene function, and phenotype have been of interest from an evolutionary perspective, to understand functional genome evolution and its phenotypic consequences, and from a biomedical perspective to understand the shared and human-specific roots of health and disease. Progress in methods for characterizing genetic, transcriptomic, and DNA methylation (DNAm) variation is driving the rapid development of extensive omics resources, which are now increasingly available from humans as well as a growing number of nonhuman primates (NHPs). The fast growth of large-scale genomic data is driving the emergence of integrated tools and databases, thus facilitating studies of gene functionality across primates. This review describes NHP genomic resources that can aid in exploration of how genes shape primate phenotypes. It focuses on the gene expression trajectories across development in different tissues, the identification of functional genetic variation (including variants deleterious for protein function and regulatory variants modulating gene expression), and DNAm profiles as an emerging tool to understand the process of aging. These resources enable comparative functional genomics approaches to identify species-specific and primate-shared gene functionalities associated with health and development.

The Cynomolgus Macaque MHC Polymorphism in Experimental Medicine

Among the non-human primates used in experimental medicine, cynomolgus macaques (Macaca fascicularis hereafter referred to as Mafa) are increasingly selected for the ease with which they are maintained and bred in captivity. Macaques belong to Old World monkeys and are phylogenetically much closer to humans than rodents, which are still the most frequently used animal model. Our understanding of the Mafa genome has progressed rapidly in recent years and has greatly benefited from the latest technical advances in molecular genetics. Cynomolgus macaques are widespread in Southeast Asia and numerous studies have shown a distinct genetic differentiation of continental and island populations. The major histocompatibility complex of cynomolgus macaque (Mafa MHC) is organized in the same way as that of human, but it differs from the latter by its high degree of classical class I gene duplication. Human polymorphic MHC regions play a pivotal role in allograft transplantation and have been associated with more than 100 diseases and/or phenotypes. The Mafa MHC polymorphism similarly plays a crucial role in experimental allografts of organs and stem cells. Experimental results show that the Mafa MHC class I and II regions influence the ability to mount an immune response against infectious pathogens and vaccines. MHC also affects cynomolgus macaque reproduction and impacts on numerous biological parameters. This review describes the Mafa MHC polymorphism and the methods currently used to characterize it. We discuss some of the major areas of experimental medicine where an effect induced by MHC polymorphism has been demonstrated.

Cynomolgus macaque IL37 polymorphism and control of SIV infection

The association between gene polymorphisms and plasma virus load at the set point (SP-PVL) was investigated in Mauritian macaques inoculated with SIV. Among 44 macaques inoculated with 50 AID50, six individuals were selected: three with SP-PVL among the highest and three with SP-PVL among the lowest. The exons of 390 candidate genes of these six animals were sequenced. Twelve non-synonymous single nucleotide polymorphisms (NS-SNPs) lying in nine genes potentially associated with PVL were genotyped in 23 animals. Three NS-SNPs with probabilities of association with PVL less than 0.05 were genotyped in a total of 44 animals. One NS-SNP lying in exon 1 of the IL37 gene displayed a significant association (p = 3.33 × 10⁻⁴) and a strong odds ratio (19.52). Multiple linear regression modeling revealed three significant predictors of SP-PVL, including the IL37 exon 1 NS-SNP (p = 0.0004) and the MHC Class IB haplotypes M2 (p = 0.0007) and M6 (p = 0.0013). These three factors in conjunction explained 48% of the PVL variance (p = 4.8 × 10⁻⁶). The potential role of IL37 in the control of SIV infection is discussed.

RNA sequencing of kidney and liver transcriptome obtained from wild cynomolgus macaque (Macaca fascicularis) originating from Peninsular Malaysia

Objective

Using high-throughput RNA sequencing technology, this study aimed to sequence the transcriptome of kidney and liver tissues harvested from Peninsular Malaysia cynomolgus macaque (Macaca fascicularis). M. fascicularis are significant nonhuman primate models in the biomedical field, owing to the macaque’s biological similarities with humans. The additional transcriptomic dataset will supplement the previously described Peninsular Malaysia M. fascicularis transcriptomes obtained in a past endeavour.

Results

A total of 75,350,240 sequence reads were obtained via Hi-seq 2500 sequencing technology. A total of 5473 significant differentially expressed genes were called. Gene ontology functional categorisation showed that cellular process, catalytic activity, and cell part categories had the highest number of expressed genes, while the metabolic pathways category possessed the highest number of expressed genes in the KEGG pathway analysis. The additional sequence dataset will further enrich existing M. fascicularis transcriptome assemblies, and provide a dataset for further downstream studies.

Electronic supplementary material

The online version of this article (10.1186/s13104-018-4014-1) contains supplementary material, which is available to authorized users.

Identification and characterization of short tandem repeats in the Tibetan macaque genome based on resequencing data

The Tibetan macaque, which is endemic to China, is currently listed as a Near Endangered primate species by the International Union for Conservation of Nature (IUCN). Short tandem repeats (STRs) refer to repetitive elements of genome sequence that range in length from 1-6 bp. They are found in many organisms and are widely applied in population genetic studies. To clarify the distribution characteristics of genome-wide STRs and understand their variation among Tibetan macaques, we conducted a genome-wide survey of STRs with next-generation sequencing of five macaque samples. A total of 1 077 790 perfect STRs were mined from our assembly, with an N50 of 4 966 bp. Mono-nucleotide repeats were the most abundant, followed by tetra- and di-nucleotide repeats. Analysis of GC content and repeats showed consistent results with other macaques. Furthermore, using STR analysis software (lobSTR), we found that the proportion of base pair deletions in the STRs was greater than that of insertions in the five Tibetan macaque individuals (P<0.05, t-test). We also found a greater number of homozygous STRs than heterozygous STRs (P<0.05, t-test), with the Emei and Jianyang Tibetan macaques showing more heterozygous loci than Huangshan Tibetan macaques. The proportion of insertions and mean variation of alleles in the Emei and Jianyang individuals were slightly higher than those in the Huangshan individuals, thus revealing differences in STR allele size between the two populations. The polymorphic STR loci identified based on the reference genome showed good amplification efficiency and could be used to study population genetics in Tibetan macaques. The neighbor-joining tree classified the five macaques into two different branches according to their geographical origin, indicating high genetic differentiation between the Huangshan and Sichuan populations. We elucidated the distribution characteristics of STRs in the Tibetan macaque genome and provided an effective method for screening polymorphic STRs. Our results also lay a foundation for future genetic variation studies of macaques.

Whole genome sequencing in the search for genes associated with the control of SIV infection in the Mauritian macaque model

In the Mauritian macaque experimentally inoculated with SIV, gene polymorphisms potentially associated with the plasma virus load at a set point, approximately 100 days post inoculation, were investigated. Among the 42 animals inoculated with 50 AID50 of the same strain of SIV, none of which received any preventive or curative treatment, nine individuals were selected: three with a plasma virus load (PVL) among the lowest, three with intermediate PVL values and three among the highest PVL values. The complete genomes of these nine animals were then analyzed. Initially, attention was focused on variants with a potential functional impact on protein encoding genes (non-synonymous SNPs (NS-SNPs) and splicing variants). Thus, 424 NS-SNPs possibly associated with PVL were detected. The 424 candidates SNPs were genotyped in these 42 SIV experimentally infected animals (including the nine animals subjected to whole genome sequencing). The genes containing variants most probably associated with PVL at a set time point are analyzed herein.

Ancient hybridization and admixture in macaques (genus Macaca) inferred from whole genome sequences

The evolutionary history of the stump-tailed macaque (Macaca arctoides) and its genetic relationship to other macaques is a subject of continuing controversy. Here, we have reported the first genome sequences of two stump-tailed macaques and one Assamese macaque (M. assamensis). Additionally, we have investigated the genetic diversity between macaque species and analyzed ancient hybridization events. Genome-wide analyses demonstrated that the stump-tailed macaque is more closely related to sinica species than to fascicularis/mulatta species. This topology contradicts the mitochondrial sequence-based phylogeny that places the stump-tailed macaque into the fascicularis/mulatta group. However, our results further show that stump-tailed macaques have genetic backgrounds distinct from sinica species, and present evidence of gene flows with rhesus macaques. We suggest that an ancient introgression occurred after stump-tailed macaques diverged from sinica species. The distinct gene flow between proto-arctoides and proto-mulatta resulted in the transfer of rhesus macaque-type mitochondria into proto-arctoides. The rhesus macaque-type mitochondria remained in populations because of genetic drift during the bottleneck. The PSMC results and morphological and geographic evidence are consistent with the mitochondria capture pattern in the stump-tailed macaque. The molecular clock estimates suggest that the mitochondrial transference into stump-tailed macaques occurred 0.4-1.4 million years ago. Furthermore, we detected extensive admixtures between different macaque species, indicating that gene flow has played an important role in the evolutionary history of the genus Macaca.

RNA sequencing (RNA-Seq) of lymph node, spleen, and thymus transcriptome from wild Peninsular Malaysian cynomolgus macaque (Macaca fascicularis)

The cynomolgus macaque (Macaca fascicularis) is an extensively utilised nonhuman primate model for biomedical research due to its biological, behavioural, and genetic similarities to humans. Genomic information of cynomolgus macaque is vital for research in various fields; however, there is presently a shortage of genomic information on the Malaysian cynomolgus macaque. This study aimed to sequence, assemble, annotate, and profile the Peninsular Malaysian cynomolgus macaque transcriptome derived from three tissues (lymph node, spleen, and thymus) using RNA sequencing (RNA-Seq) technology. A total of 174,208,078 paired end 70 base pair sequencing reads were obtained from the Illumina Hi-Seq 2500 sequencer. The overall mapping percentage of the sequencing reads to the M. fascicularis reference genome ranged from 53–63%. Categorisation of expressed genes to Gene Ontology (GO) and KEGG pathway categories revealed that GO terms with the highest number of associated expressed genes include Cellular process, Catalytic activity, and Cell part, while for pathway categorisation, the majority of expressed genes in lymph node, spleen, and thymus fall under the Global overview and maps pathway category, while 266, 221, and 138 genes from lymph node, spleen, and thymus were respectively enriched in the Immune system category. Enriched Immune system pathways include Platelet activation pathway, Antigen processing and presentation, B cell receptor signalling pathway, and Intestinal immune network for IgA production. Differential gene expression analysis among the three tissues revealed 574 differentially expressed genes (DEG) between lymph and spleen, 5402 DEGs between lymph and thymus, and 7008 DEGs between spleen and thymus. Venn diagram analysis of expressed genes revealed a total of 2,630, 253, and 279 tissue-specific genes respectively for lymph node, spleen, and thymus tissues. This is the first time the lymph node, spleen, and thymus transcriptome of the Peninsular Malaysian cynomolgus macaque have been sequenced via RNA-Seq. Novel transcriptomic data will further enrich the present M. fascicularis genomic database and provide future research potentials, including novel transcript discovery, comparative studies, and molecular markers development.

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.

Diversity and regulatory impact of copy number variation in the primate Macaca fascicularis

BACKGROUND

Copy number variations (CNVs) are a significant source of genetic diversity and commonly found in mammalian genomes. We have generated a genome-wide CNV map for Cynomolgus monkeys (Macaca fascicularis). This crab-eating macaque is the closest animal model to humans that is used in biomedical research.

RESULTS

We show that Cynomolgus monkey CNVs are in general much smaller in size than gene loci and are specific to the population of origin. Genome-wide expression data from five vitally important organs demonstrates that CNVs in close proximity to transcription start sites associate strongly with expression changes. Among these eQTL genes we find an overrepresentation of genes involved in metabolism, receptor activity, and transcription.

CONCLUSION

These results provide evidence that CNVs shape tissue transcriptomes in monkey populations, potentially offering an adaptive advantage. We suggest that this genetic diversity should be taken into account when using Cynomolgus macaques as models.

Single nucleotide polymorphisms in the FcγR3A and TAP1 genes impact ADCC in cynomolgus monkey PBMCs

Phenotypic variability is often observed in cynomolgus monkeys on preclinical studies and may, in part, be driven by genetic variability. However, the role of monkey genetic variation remains largely unexplored in the context of drug response. This study evaluated genetic variation in cynomolgus monkey FcγR3A and TAP1 genes and the potential impact of identified polymorphisms on antibody-dependent cell-mediated cytotoxicity (ADCC) in vitro. Studies in humans have demonstrated that a single nucleotide polymorphism (SNP), F158V, in FcγR3A can influence response to rituximab through altered ADCC and that SNPs in TAP1/2 decrease natural killer (NK) cell activity against major histocompatibility complex (MHC) class I deficient cells, potentially through altered ADCC. Monkeys were genotyped for FcγR3A and TAP1 SNPs, and ADCC was assessed in vitro using peripheral blood mononuclear cells (PBMCs) treated with trastuzumab in the presence of NCI-N87 cells. FcγR3A g.1134A>C (exonic S42R), FcγR3A g.5027A>G (intronic), and TAP1 g.1A>G (start codon loss) SNPs were all significantly associated with decreased ADCC for at least one trastuzumab concentration ≥0.0001 μM when compared with wild type (WT). Regression analysis demonstrated significant association of the SNP-SNP pairs FcγR3A g.1134A>C/TAP1 g.1A>G and FcγR3A g.5027A>G/TAP1 g.1A>G with a combinatorial decrease on ADCC. Mechanisms underlying the decreased ADCC were investigated by measuring FcγR3A/IgG binding affinity and expression of FcγR3A and TAP1 in PBMCs; however, no functional associations were observed. These data demonstrate that genetic variation in cynomolgus monkeys is reflective of known human genetic variation and may potentially contribute to variable drug response in preclinical studies.

Genetic diversity in humans and non-human primates and its evolutionary consequences

Genetic diversity is a key parameter in population genetics and is important for understanding the process of evolution and for the development of appropriate conservation strategies. Recent advances in sequencing technology have enabled the measurement of genetic diversity of various organisms at the nucleotide level and on a genome-wide scale, yielding more precise estimates than were previously achievable. In this review, I have compiled and summarized the estimates of genetic diversity in humans and non-human primates based on recent genome-wide studies. Although studies on population genetics demonstrated fluctuations in population sizes over time, general patterns have emerged. As shown previously, genetic diversity in humans is one of the lowest among primates; however, certain other primate species exhibit genetic diversity that is comparable to or even lower than that in humans. There exists greater than 10-fold variation in genetic diversity among primate species, and I found weak correlation with species fecundity but not with body or propagule size. I further discuss the potential evolutionary consequences of population size decline on the evolution of primate species. The level of genetic diversity negatively correlates with the ratio of non-synonymous to synonymous polymorphisms in a population, suggesting that proportionally greater numbers of slightly deleterious mutations segregate in small rather than large populations. Although population size decline is likely to promote the fixation of slightly deleterious mutations, there are molecular mechanisms, such as compensatory mutations at various molecular levels, which may prevent fitness decline at the population level. The effects of slightly deleterious mutations from theoretical and empirical studies and their relevance to conservation biology are also discussed in this review.

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.

Translational research on influenza virus infection using a nonhuman primate model

Influenza virus infection is a seasonal infectious disease for humans, whereas it is also a zoonosis that is originally transmitted from animals to humans. Therefore, several animal models are used in research on influenza virus infection. We have used a nonhuman primate (NHP) model to extrapolate pathogenicity of various influenza viruses and efficacy of vaccines and antiviral drugs against the influenza viruses in humans. NHPs have genes, anatomical structure, and immune responses similar to those of humans as compared to other animal models. Using an NHP model, we revealed that the pandemic 2009 influenza A virus caused viral pneumonia as reported in human patients. Thus, it is thought that NHP models can be used to predict replication of emerging viruses in humans. We also examined the pathogenicity of highly pathogenic avian influenza viruses and evaluated a new therapeutic antibody in macaques under an immunocompromised condition. NHP models have provided promising results in research on other infectious diseases including Ebola virus and human/simian immunodeficiency virus infections. Thus, NHPs are important in biomedical research for determining the pathogenesis and for development of treatments, especially when clinical trials are difficult. We summarize the characteristics and advantages of research using NHP models in this review.

The conservation genetics juggling act: integrating genetics and ecology, science and policy

The field of conservation genetics, when properly implemented, is a constant juggling act integrating molecular genetics, ecology, and demography with applied aspects concerning managing declining species or implementing conservation laws and policies. This young field has grown substantially since the 1980s following the development of polymerase chain reaction and now into the genomics era. Our laboratory has ‘grown up’ with the field, having worked on these issues for over three decades. Our multidisciplinary approach entails understanding the behavior and ecology of species as well as the underlying processes that contribute to genetic viability. Taking this holistic approach provides a comprehensive understanding of factors that influence species persistence and evolutionary potential while considering annual challenges that occur throughout their life cycle. As a federal laboratory, we are often addressing the needs of the U.S. Fish and Wildlife Service in their efforts to list, de‐list, or recover species. Nevertheless, there remains an overall communication gap between research geneticists and biologists who are charged with implementing their results. Therefore, we outline the need for a National Center for Small Population Biology to ameliorate this problem and provide organizations charged with making status decisions firmer ground from which to make their critical decisions.

Identification of microRNAs in Macaca fascicularis (Cynomolgus Monkey) by Homology Search and Experimental Validation by Small RNA-Seq and RT-qPCR Using Kidney Cortex Tissues

MicroRNAs (miRNAs) present in tissues and biofluids are emerging as sensitive and specific safety biomarkers. MiRNAs have not been thoroughly described in M. fascicularis, an animal model used in pharmaceutical industry especially in drug safety evaluation. Here we investigated the miRNAs in M. fascicularis. For Macaca mulatta, a closely related species of M. fascicularis, 619 stem-loop precursor miRNAs (pre-miRNAs) and 914 mature miRNAs are available in miRBase version 21. Using M. mulatta miRNAs as a reference list and homology search tools, we identified 604 pre-miRNAs and 913 mature miRNAs in the genome of M. fascicularis. In order to validate the miRNAs identified by homology search we attempted to sequence miRNAs expressed in kidney cortex from M. fascicularis. MiRNAs expressed in kidney cortex may indeed be released in urine upon kidney cortex damage and be potentially used to monitor drug induced kidney injury. Hence small RNA sequencing libraries were prepared using kidney cortex tissues obtained from three naive M. fascicularis and sequenced. Analysis of sequencing data indicated that 432 out of 913 mature miRNAs were expressed in kidney cortex tissues. Assigning these 432 miRNAs to pre-miRNAs revealed that 273 were expressed from both the -5p and -3p arms of 150 pre-miRNAs and 159 miRNAs expressed from either the -5p or -3p arm of 176 pre-miRNAs. Mapping sequencing reads to pre-miRNAs also facilitated the detection of twenty-two new miRNAs. To substantiate miRNAs identified by small RNA sequencing, 313 miRNAs were examined by RT-qPCR. Expression of 262 miRNAs in kidney cortex tissues ware confirmed by TaqMan microRNA RT-qPCR assays. Analysis of kidney cortex miRNA targeted genes suggested that they play important role in kidney development and function. Data presented in this study may serve as a valuable resource to assess the renal safety biomarker potential of miRNAs in Cynomolgus monkeys.

Mx1 and Mx2 key antiviral proteins are surprisingly lost in toothed whales

Viral outbreaks in dolphins and other Delphinoidea family members warrant investigation into the integrity of the cetacean immune system. The dynamin-like GTPase genes Myxovirus 1 (Mx1) and Mx2 defend mammals against a broad range of viral infections. Loss of Mx1 function in human and mice enhances infectivity by multiple RNA and DNA viruses, including orthomyxoviruses (influenza A), paramyxoviruses (measles), and hepadnaviruses (hepatitis B), whereas loss of Mx2 function leads to decreased resistance to HIV-1 and other viruses. Here we show that both Mx1 and Mx2 have been rendered nonfunctional in Odontoceti cetaceans (toothed whales, including dolphins and orcas). We discovered multiple exon deletions, frameshift mutations, premature stop codons, and transcriptional evidence of decay in the coding sequence of both Mx1 and Mx2 in four species of Odontocetes. We trace the likely loss event for both proteins to soon after the divergence of Odontocetes and Mystocetes (baleen whales) ∼33-37 Mya. Our data raise intriguing questions as to what drove the loss of both Mx1 and Mx2 genes in the Odontoceti lineage, a double loss seen in none of 56 other mammalian genomes, and suggests a hitherto unappreciated fundamental genetic difference in the way these magnificent mammals respond to viral infections.