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

ABSTRACT: BACKGROUND: The genetic background of the cynomolgus macaque (Macaca fascicularis) is made complex by the high genetic diversity, population structure, and gene introgression from the closely related rhesus macaque (Macaca mulatta). Herein we report the whole-genome sequence of a Malaysian cynomolgus macaque male with more than 40-fold coverage, which was determined using a resequencing method based on the Indian rhesus macaque genome. RESULTS: We identified approximately 9.7 million single nucleotide variants (SNVs) between the Malaysian cynomolgus and the Indian rhesus macaque genomes. Compared with humans, a smaller nonsynonymous/synonymous SNV ratio in the cynomolgus macaque suggests more effective removal of slightly deleterious mutations. Comparison of two cynomolgus (Malaysian and Vietnamese) and two rhesus (Indian and Chinese) macaque genomes, including previously published macaque genomes, suggests that Indochinese cynomolgus macaques have been more affected by gene introgression from rhesus macaques. We further identified 60 nonsynonymous SNVs that completely differentiated the cynomolgus and rhesus macaque genomes, and that could be important candidate variants for determining species-specific responses to drugs and pathogens. The demographic inference using the genome sequence data revealed that Malaysian cynomolgus macaques have experienced at least three population bottlenecks. CONCLUSIONS: This list of whole-genome SNVs will be useful for many future applications, such as an array-based genotyping system for macaque individuals. High-quality whole-genome sequencing of the cynomolgus macaque genome may aid studies on finding genetic differences that are responsible for phenotypic diversity in macaques and may help control genetic backgrounds among individuals.     

Molecular Evidence for Distinct Genotypes of Monkey B Virus (Herpesvirus Simiae) Which Are Related to the Macaque Host Species

Although monkey B virus (herpesvirus simiae; BV) is common in all macaque species, fatal human infections appear to be associated with exposure to rhesus macaques (Macaca mulatta), suggesting that BV isolates from rhesus monkeys may be more lethal to nonmacaques than are BV strains indigenous to other macaque species. To determine if significant differences that would support this supposition exist among BV isolates, we compared multiple BV strains isolated from rhesus, cynomolgus, pigtail, and Japanese macaques. Antigenic analyses indicated that while the isolates were very closely related to one another, there are some antigenic determinants that are specific to BV isolates from different macaque species. Restriction enzyme digest patterns of viral DNA revealed marked similarities between rhesus and Japanese macaque isolates, while pigtail and cynomolgus macaque isolates had distinctive cleavage patterns. To further compare genetic diversity among BV isolates, DNA sequences from two regions of the viral genome containing genes that are conserved (UL27 and US6) and variable (US4 and US5) among primate alphaherpesviruses, as well as from two noncoding intergenic regions, were determined. From these sequence data and a phylogenetic analysis of them it was evident that while all isolates were closely related strains of BV, there were three distinct genotypes. The three BV genotypes were directly related to the macaque species of origin and were composed of (i) isolates from rhesus and Japanese macaques, (ii) cynomolgus monkey isolates, and (iii) isolates from pigtail macaques. This study demonstrates the existence of different BV genotypes which are related to the macaque host species and thus provides a molecular basis for the possible existence of BV isolates which vary in their levels of pathogenicity for nonmacaque species.     

A BAC-based contig map of the cynomolgus macaque (Macaca fascicularis) major histocompatibility complex genomic region

The construction of a cynomolgus macaque (Macaca fascicularis, Mafa) BAC library for genomic comparison between rhesus and cynomolgus macaques is necessary to promote the cynomolgus macaque as one of the important experimental animals for future medical and biological research. In this paper, we constructed a cynomolgus macaque BAC library and a map of the MHC (Mafa) genomic region for comparison of the genomic organization and nucleotide similarities between the human, the chimpanzee, and the rhesus macaque. The BAC library consists of 221,184 clones with an average insert size of 83 kb, providing a sixfold coverage of the haploid genome. A total of 114 BAC clones and 54 PCR primer sets were used to construct a 4.3-Mb contig of the MHC region. Diversity analysis of genomic sequence from selected subregions of the MHC revealed that the cynomolgus sequence varied compared to rhesus macaque, human, and chimpanzee sequences by 0.48, 4.15, and 4.10%, respectively. From these findings, we conclude that the BAC library and Mafa genomic map are useful tools for genome analysis and will have important applications for comparative genomics and identifying regions of consequence in medical research.     

Pathogenicity of simian-human immunodeficiency virus SHIV-89.6P and SIVmac is attenuated in cynomolgus macaques and associated with early T-lymphocyte responses

Because most studies of AIDS pathogenesis in nonhuman primates have been performed in Indian-origin rhesus macaques (Macaca mulatta), little is known about lentiviral pathogenicity and control of virus replication following infection of alternative macaque species. Here, we report the consequences of simian-human immunodeficiency virus SHIV-89.6P and SIVmac251 infection in cynomolgus (Macaca fascicularis) and rhesus macaques of Chinese origin. Compared to the pathogenicity of the same viruses in Indian rhesus macaques, both cynomolgus and Chinese rhesus macaques showed lower levels of plasma virus. By 9 to 10 months after infection, both viruses became undetectable in plasma more frequently in cynomolgus than in either Chinese or Indian rhesus macaques. Furthermore, after SHIV-89.6P infection, CD4+ T-cell numbers declined less and survival was longer in cynomolgus and Chinese rhesus macaques than in Indian rhesus macaques. This attenuated pathogenicity was associated with gamma interferon ELISPOT responses to Gag and Env that were generated earlier and of higher frequency in cynomolgus than in Indian rhesus macaques. Cynomolgus macaques also developed higher titer neutralizing antibodies against SHIV-89.6 at 10 and 20 weeks postinoculation than Indian rhesus macaques. These studies demonstrate that the pathogenicity of nonhuman primate lentiviruses varies markedly based on the species or geographic origin of the macaques infected and suggest that the cellular immune responses may contribute to the control of pathogenicity in cynomolgus macaques. While cynomolgus and Chinese rhesus macaques provide alternative animal models of lentiviral infection, the lower levels of viremia in cynomolgus macaques limit the usefulness of infection of this species for vaccine trials that utilize viral load as an experimental endpoint.     

猕猴MHC-DPB1基因外显子2多态性研究

猕猴(Macaca mulatta)是最理想的医学实验灵长类动物, 且为国家二级保护动物。为了解中国猕猴主要组织相容复合体(Major histocompatibility complex, MHC)基因的遗传多态性背景, 为它们在生物医学研究中的应用及其遗传资源的保护提供一定的科学依据, 文章采用变性梯度凝胶电泳(Denaturing gradient gel electrophoresis, DGGE)和克隆测序技术分析了106个四川野生猕猴MHC-DPB1基因的exon 2, 共检测到21个Mamu-DPB1等位基因, 其中有15个为本研究中首次发现的新等位基因; 从整个大的猕猴群体(106个个体)来看, 等位基因频率最高的是Mamu-DPB1*30(0.1120); 单独从不同地理群体来看, 最高等位基因频率分别为: 小金-DPB1*30 (0.1120), 黑水-DPB1*04 (0.1702), 巴中-DPB1*32 (0.1613), 汉源-DPB1*30(0.1120), 九龙-DPB1*04(0.1139); 氨基酸序列比对发现, 猕猴Mamu-DPB1等位基因编码的氨基酸序列中, 有12个氨基酸残基变异位点表现出物种特异性, 其中有9个位于新发现的15个Mamu-DPB1等位基因氨基酸序列中; 不同物种来源的DPB1等位基因系统发生树表明, 猕猴与其近缘物种食蟹猴(Macaca fascicularis)的DPB1等位基因间存在着跨种多态(Trans-species polymorphism)现象。研究还表明, MHC-DPB1等位基因在中国猕猴群体和先前为主要研究对象的印度猕猴群体间具有较大的差异。     

Development of an integrative database with 499 novel microsatellite markers for Macaca fascicularis

Background

Cynomolgus macaques (Macaca fascicularis) are a valuable resource for linkage studies of genetic disorders, but their microsatellite markers are not sufficient. In genetic studies, a prerequisite for mapping genes is development of a genome-wide set of microsatellite markers in target organisms. A whole genome sequence and its annotation also facilitate identification of markers for causative mutations. The aim of this study is to establish hundreds of microsatellite markers and to develop an integrative cynomolgus macaque genome database with a variety of datasets including marker and gene information that will be useful for further genetic analyses in this species.

Results

We investigated the level of polymorphisms in cynomolgus monkeys for 671 microsatellite markers that are covered by our established Bacterial Artificial Chromosome (BAC) clones. Four hundred and ninety-nine (74.4%) of the markers were found to be polymorphic using standard PCR analysis. The average number of alleles and average expected heterozygosity at these polymorphic loci in ten cynomolgus macaques were 8.20 and 0.75, respectively.

Conclusion

BAC clones and novel microsatellite markers were assigned to the rhesus genome sequence and linked with our cynomolgus macaque cDNA database (QFbase). Our novel microsatellite marker set and genomic database will be valuable integrative resources in analyzing genetic disorders in cynomolgus macaques.     

Single‐Nucleotide Polymorphisms Reveal Patterns of Allele Sharing Across the Species Boundary Between Rhesus (Macaca mulatta) and Cynomolgus (M. fascicularis) Macaques

Both phenotypic and genetic evidence for asymmetric hybridization between rhesus (Macaca mulatta) and cynomolgus (Macaca fascicularis) macaques has been observed in the region of Indochina where both species are sympatric. The large‐scale sharing of major histocompatibility complex (MHC) class II alleles between the two species in this region supports the hypothesis that genes, and especially genes involved in immune response, are being transferred across the species boundary. This differential introgression has important implications for the incorporation of cynomolgus macaques of unknown geographic origin in biomedical research protocols. Our study found that for 2,808 single‐nucleotide polymorphism (SNP) markers, the minor allele frequencies (MAF) and observed heterozygosity calculated from a sample of Vietnamese cynomolgus macaques was significantly different from those calculated from samples of both Chinese rhesus and Indonesian cynomolgus macaques. SNP alleles from Chinese rhesus macaques were overrepresented in a sample of Vietnamese cynomolgus macaques relative to their Indonesian conspecifics and located in genes functionally related to the primary immune system. These results suggest that Indochinese cynomolgus macaques represent a genetically and immunologically distinct entity from Indonesian cynomolgus macaques. Am. J. Primatol. 75:135‐144, 2013. © 2012 Wiley Periodicals, Inc.     

MHC class I A region diversity and polymorphism in macaque species

The HLA-A locus represents a single copy gene that displays abundant allelic polymorphism in the human population, whereas, in contrast, a nonhuman primate species such as the rhesus macaque (Macaca mulatta) possesses multiple HLA-A-like (Mamu-A) genes, which parade varying degrees of polymorphism. The number and combination of transcribed Mamu-A genes present per chromosome display diversity in a population of Indian animals. At present, it is not clearly understood whether these different A region configurations are evolutionarily stable entities. To shed light on this issue, rhesus macaques from a Chinese population and a panel of cynomolgus monkeys (Macaca fascicularis) were screened for various A region-linked variations. Comparisons demonstrated that most A region configurations are old entities predating macaque speciation, whereas most allelic variation (>95%) is of more recent origin. The latter situation contrasts the observations of the major histocompatibility complex class II genes in rhesus and cynomolgus macaques, which share a high number of identical alleles (>30%) as defined by exon 2 sequencing.     

Assessing natural introgression in 2 biomedical model species, the rhesus macaque (Macaca mulatta) and the long-tailed macaque (Macaca fascicularis)

Rhesus macaque (Macaca mulatta) and long-tailed macaque (Macaca fascicularis) are the 2 most commonly used primate model species in biomedical sciences. Although morphological studies have revealed a weak hybridization at the interspecific contact zone, in the north of Indochina, a molecular study has suggested an ancient introgression from rhesus to long-tailed macaque into the Indo-Chinese peninsula. However, the gene flow between these 2 taxa has never been quantified using genetic data and theoretical models. In this study, we have examined genetic variation within and between the parapatric Chinese rhesus macaque and Indo-Chinese long-tailed macaque populations, using 13 autosomal, 5 sex-linked microsatellite loci and mitochondrial DNA sequence data. From these data, we assessed genetic structure and estimated gene flow using a Bayesian clustering approach and the "Isolation with Migration" model. Our results reveal a weak interspecific genetic differentiation at both autosomal and sex-linked loci, suggesting large population sizes and/or gene flow between populations. According to the Bayesian clustering, Chinese rhesus macaque is a highly homogeneous gene pool that contributes strongly to the current Indo-Chinese long-tailed macaque genetic makeup, whether or not current admixture is assumed. Coalescent simulations, which integrated the characteristics of the loci, pointed out 1) a higher effective population size in rhesus macaque, 2) no mitochondrial gene flow, and 3) unilateral and male-mediated nuclear gene flow of approximately 10 migrants per generation from rhesus to long-tailed macaque. These patterns of genetic structure and gene flow suggest extensive ancient introgression from Chinese rhesus macaque into the Indo-Chinese long-tailed macaque population.     

Detecting lineage-specific adaptive evolution of brain-expressed genes in human using rhesus macaque as outgroup

Comparative genetic analysis between human and chimpanzee may detect genetic divergences responsible for human-specific characteristics. Previous studies have identified a series of genes that potentially underwent Darwinian positive selection during human evolution. However, without a closely related species as outgroup, it is difficult to identify human-lineage-specific changes, which is critical in delineating the biological uniqueness of humans. In this study, we conducted phylogeny-based analyses of 2633 human brain-expressed genes using rhesus macaque as the outgroup. We identified 47 candidate genes showing strong evidence of positive selection in the human lineage. Genes with maximal expression in the brain showed a higher evolutionary rate in human than in chimpanzee. We observed that many immune-defense-related genes were under strong positive selection, and this trend was more prominent in chimpanzee than in human. We also demonstrated that rhesus macaque performed much better than mouse as an outgroup in identifying lineage-specific selection in humans.     

Use of Genome-wide Heterospecific Single-Nucleotide Polymorphisms to Estimate Linkage Disequilibrium in Rhesus and Cynomolgus Macaques

Rhesus and cynomolgus macaques are frequently used in biomedical research, and the availability of their reference genomes now provides for their use in genome-wide association studies. However, little is known about linkage disequilibrium (LD) in their genomes, which can affect the design and success of such studies. Here we studied LD by using 1781 conserved single-nucleotide polymorphisms (SNPs) in 183 rhesus macaques (Macaca mulatta), including 97 purebred Chinese and 86 purebred Indian animals, and 96 cynomolgus macaques (M. fascicularis fascicularis). Correlation between loci pairs decayed to 0.02 at 1146.83, 2197.92, and 3955.83 kb for Chinese rhesus, Indian rhesus, and cynomolgus macaques, respectively. Differences between the observed heterozygosity and minor allele frequency (MAF) of pairs of these 3 taxa were highly statistically significant. These 3 nonhuman primate taxa have significantly different genetic diversities (heterozygosity and MAF) and rates of LD decay. Our study confirms a much lower rate of LD decay in Indian than in Chinese rhesus macaques relative to that previously reported. In contrast, the especially low rate of LD decay in cynomolgus macaques suggests the particular usefulness of this species in genome-wide association studies. Although conserved markers, such as those used here, are required for valid LD comparisons among taxa, LD can be assessed with less bias by using species-specific markers, because conserved SNPs may be ancestral and therefore not informative for LD.Abbreviations: GWAS, genome-wide association study; LD, linkage disequilibrium; MAF, minor allele frequencyContributing to the widespread use of nonhuman primates in biomedical research, captive-breeding programs such as those of the National Primate Research Center system in the United States were established initially by using animals imported from Asia. The 2 most commonly used primates are rhesus macaques (Macaca mulatta) and long-tailed or cynomolgus macaques (M. fascicularis fascicularis).After humans, rhesus macaques are the most widely distributed primate species.^37,38 This species is found throughout mainland Asia, ranging from Afghanistan to India and eastward through Thailand and southern China to the Yellow Sea.^31,34 In addition to their significant morphological differences,⁹ rhesus macaques of Indian and Chinese origins have been demonstrated to exhibit significant phenotypic differences that are directly relevant to their use as biomedical models in experimental studies.^2,23,42 Cynomolgus macaques are found south of the subtropical and temperate geographic distributions of rhesus macaques, in the south and southeast Indo-Malayan regions.^8,10The 2 species share a common ancestor that lived 1 to 2 million years ago.^3,13,25 This ancestral population of rhesus macaques diverged from a fascicularis-like ancestor shared in common with both rhesus and cynomolgus macaques after cynomolgus macaques expanded from their homeland in Indonesia.³⁶ For this reason, genetic markers present in Indian rhesus macaques are either highly derived or are conserved as ancestral markers shared with Chinese rhesus macaques. The interspecific boundaries of rhesus and cynomolgus macaques are delineated by a narrow zone of parapatry in northern Indochina,^7,8,10 within which male-biased gene flow^37,39 and relatively high, but highly variable, levels of introgression of genes³² have occurred from rhesus to cynomolgus macaque groups.^37,39 Because cynomolgus macaques originated in Indonesia³⁶ and because rhesus macaques probably diverged from cynomolgus macaques in southwestern China,¹¹ genetic markers shared between Indonesian cynomolgus macaques and Chinese rhesus macaques comprise a unique set of markers that are conserved in both macaque species.The wide assortment of morphometric differences^8,9 and the broad geographic distribution of these 2 macaque species foster an expectation of high genetic diversity within and between them that could be exploited for mapping genes responsible for phenotypic differences between taxa. A better understanding of linkage disequilibrium (LD) in these nonhuman primate species can lead to a more informed selection of study subjects for, and more efficient conduct of, genome-wide association studies (GWAS) of particular diseases that macaques share in common with humans. LD is the nonrandom association of alleles at 2 or more adjacent loci that descend from single, ancestral chromosomes.²⁹ LD plays a critical role in gene mapping, both as a tool for fine mapping of complex disease genes and in GWAS-based approaches. GWAS facilitate the identification of genes associated with complex and common traits or diseases by examining LD estimates among large numbers of common genetic variants, typically single-nucleotide polymorphisms (SNPs), between pairs of different groups of subjects to determine whether any variant is associated with a trait or disease of interest. LD data make tightly linked variants strongly correlated to produce successful association studies. For instance, LD reduces the number of markers and sample size of study subjects required to map genes influencing phenotypes to the genome because markers in LD are linked and inherited together.¹³ In addition, differences in LD can be used to identify orthologs for detecting the signatures of selective sweeps,²¹ as defined by d_N/d_S ratios obtained through the McDonald–Kreitman neutrality test.²⁴ Furthermore, LD assessments can provide a more complete understanding of genome structure by defining the boundaries of haplotype blocks, within which recombination is rare or absent and which are separated by recombination ‘hotspots,’ in genomes.⁴³Evidence from a study based on 1476 SNPs identified in ENCODE regions of the Indian rhesus macaque genome¹³ indicated that the rate of LD decay is higher in Chinese than in Indian rhesus macaques due to an hypothesized genetic bottleneck experienced by Indian rhesus macaques after diverging from the eastern subspecies, and, therefore, that Indian rhesus macaques, having higher LD, may be more useful for GWAS than Chinese rhesus macaques. In that study,¹³ only 33% of the SNPs were shared in common between the 2 subspecies, with Chinese rhesus macaques contributing to more than 60% of the remaining rhesus SNPs. Conversely, another study⁴¹ reported a slower rate of decay of LD in 25 Chinese than in 25 Indian rhesus macaques on the basis of 4040 SNPs, only 2% of which fell in coding regions, but 68% of those SNPs were shared between the 2 subspecies, with Indian rhesus macaques contributing almost 60% of the remaining SNPs. The marked disparity between the 2 studies in the proportions of shared SNPs used, the subspecies with the most genetic diversity, the sample size of Chinese rhesus macaques, the proportions of SNPs located in or near coding regions that are subject to functional constraints, and the greater disparity in LD decay between the 2 subspecies of rhesus macaques might reflect biases in either or both studies. For example, the use of markers whose frequencies are uncharacteristically low in one subspecies relative to the other can underestimate the rate of LD decay because lower frequency alleles, on average, are younger and have experienced less time for recombination.²⁶ To avoid the influence of such ascertainment biases, comparisons of LD between 2 taxa should involve only SNPs conserved in both taxa. Moreover, because 2 points do not provide a phylogenetic or cladistic analysis to assign specific SNPs to origin on one phylogenetic line or another, comparing just the Indian and Chinese rhesus macaques without an additional primate taxon makes it is difficult to establish polarity and distinguish between derived and conserved SNPs. This limitation likely led to the contradictory conclusions of the 2 previously cited studies^13,41 regarding the rate of LD decay in Chinese and Indian rhesus macaques.Because rhesus and cynomolgus macaques share a common fascicularis-like ancestor, a comparison of heterospecific SNPs among cynomolgus, Indian rhesus, and Chinese rhesus macaques would likely be fundamental to inferences regarding genome-wide LD estimates. The objective of the present study was to evaluate the conclusions of previous studies^13,41 by using our panel of 1781 autosomal SNPs that are conserved in both rhesus and cynomolgus macaques to estimate the rates at which genome-wide LD decays in Indian and Chinese rhesus macaques and cynomolgus macaques, the species ancestral to rhesus macaques, and to evaluate the suitability of these populations for GWAS.     

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.     

Metagenomic comparison of the rectal microbiota between rhesus macaques(Macaca mulatta) and cynomolgus macaques(Macaca fascicularis)

Rhesus macaques(Macaca mulatta) and cynomolgus macaques(Macaca fascicularis) are frequently used in establishing animal models for human diseases. To determine the differences in gut microbiota between these species, rectal swabs from 20 rhesus macaques and 21 cynomolgus macaques were collected, and the microbial composition was examined by deep sequencing of the 16 S rR NA gene. We found that the rectal microbiota of cynomolgus macaques exhibited significantly higher alpha diversity than that of rhesus macaques, although the observed number of operational taxonomic units(OTUs) was almost the same. The dominant taxa at both the phylum and genus levels were similar between the two species, although the relative abundances of these dominant taxa were significantly different between them. Phylogenetic Investigation of Communities by Reconstruction of Unobserved States(PICRUSt) showed significant differences in the functional components between the microbiota of the two species, in particular the lipopolysaccharide(LPS) synthesis proteins. The above data indicated significant differences in microbial composition and function between these two closely related macaque species, which should be taken into consideration in the future selection of these animals for disease models.     

A radiation hybrid mapping panel for the rhesus macaque

The genomes of nonhuman primates have recently become highly visible candidates for full genome analysis, as they provide powerful models of human disease and a better understanding of the evolution of the human genome. We describe the creation of a 5000 rad radiation hybrid (RH) mapping panel for the rhesus macaque. Duplicate genotypes of 84 microsatellite and coding gene sequence tagged sites from six macaque chromosomes produced an estimated whole genome retention frequency of 0.33. To test the mapping ability of the panel, we constructed RH maps for macaque chromosomes 7 and 9 and compared them to orthologous locus orders in existing human and baboon maps derived from different methodologies. Concordant marker order between all three species maps suggests that the current panel represents a powerful mapping resource for generating high-density comparative maps of the rhesus macaque and other species genomes.     

High-throughput single-nucleotide polymorphism discovery and the search for candidate genes for long-term SIVmac nonprogression in Chinese rhesus macaques (Macaca mulatta)

Background Genetic differences between Indian and Chinese rhesus macaques contribute to the phenotypic variance of clinical trials, including infection with SIVmac. The completion of the rhesus genome has facilitated the discovery of several thousand markers. Methods We developed a genome‐wide SNP map for rhesus macaques containing 3869 validated markers with an average distance of 0.88 Mb and used the program VarLD to identify genomic areas with significant differences in linkage disequilibrium (LD) between Indian‐derived and Chinese rhesus macaques. Results Forty‐one statistically significant differences in LD between Chinese and Indian‐origin rhesus were detected on chromosomes 1, 4, 5 and 11. The region of greatest LD difference was located on the proximal end of chromosome one, which also contained the genes ELAVL4, MAST2 and HIVEP3. Conclusion These genomic areas provide entry to more detailed studies of gene function. This method is also applicable to the study of differences in biomarkers between regional populations of other species.     

Simian immunodeficiency virus SIVmac239 infection of major histocompatibility complex-identical cynomolgus macaques from Mauritius

Nonhuman primates are widely used to study correlates of protective immunity in AIDS research. Successful cellular immune responses have been difficult to identify because heterogeneity within macaque major histocompatibility complex (MHC) genes results in quantitative and qualitative differences in immune responses. Here we use microsatellite analysis to show that simian immunodeficiency virus (SIV)-susceptible cynomolgus macaques (Macaca fascicularis) from the Indian Ocean island of Mauritius have extremely simple MHC genetics, with six common haplotypes accounting for two-thirds of the MHC haplotypes in feral animals. Remarkably, 39% of Mauritian cynomolgus macaques carry at least one complete copy of the most frequent MHC haplotype, and 8% of these animals are homozygous. In stark contrast, entire MHC haplotypes are rarely conserved in unrelated Indian rhesus macaques. After intrarectal infection with highly pathogenic SIVmac239 virus, a pair of MHC-identical Mauritian cynomolgus macaques mounted concordant cellular immune responses comparable to those previously reported for a pair of monozygotic twins infected with the same strain of human immunodeficiency virus. Our identification of relatively abundant SIV-susceptible, MHC-identical macaques will facilitate research into protective cellular immunity.     

Mapping cynomolgus monkey MHC class I district on chromosome 6p13 using pooled cDNAs.

The cynomolgus monkey (Macaca fascicularis) is a frequently used animal model for studying human diseases, especially immune related ones. For a better understanding of its major histocompatibility complex (MHC) class I district chromosome location, we selected seven cDNA clones as probes for fluorescence in situ hybridization (FISH) from a lymphocyte cell line cDNA library. Expressed sequence tags (ESTs) from these clones were assembled into three clusters and annotated Mafa-A and Mafa-B genes. Further bioinformatics analysis shows that they had multiple duplications spanning approximately 2.8 Mb on the rhesus macaque MHC class I district. Using the FISH technique, we mapped the seven pooled cDNA clones to the short arm of the cynomolgus monkey chromosome 6 on 6p13. To our knowledge, this is the first report of the location of cynomolgus monkey MHC class I district. Using pooled adjacent cDNAs as probes also allows affordable, specific genome region mapping research.     

Interspecies pair housing of macaques in a research facility

The eighth edition of The Guide for the Care and Use of Laboratory Animals establishes social housing as the 'default' for social species including non-human primates. The advantages of social housing for primates have been well established, but small research facilities housing few primates in indoor cages have struggled with social housing as a result of limitations on appropriate housing and availability of compatible monkeys. Here, we report a novel approach to pair housing macaques - crossing species. We have successfully pair housed an intact male rhesus macaque with an intact male cynomolgus macaque, and an adult female rhesus macaque with numerous subadult female cynomolgus macaques. Monkeys in these pairs established dominant-subordinate relationships similar to same-species pairs. Rhesus and cynomolgus macaques can be successfully paired for the purpose of social housing in facilities with limited numbers of monkeys.     

Characterization of 47 MHC class I sequences in Filipino cynomolgus macaques

Cynomolgus macaques (Macaca fascicularis) provide increasingly common models for infectious disease research. Several geographically distinct populations of these macaques from Southeast Asia and the Indian Ocean island of Mauritius are available for pathogenesis studies. Though host genetics may profoundly impact results of such studies, similarities and differences between populations are often overlooked. In this study we identified 47 full-length MHC class I nucleotide sequences in 16 cynomolgus macaques of Filipino origin. The majority of MHC class I sequences characterized (39 of 47) were unique to this regional population. However, we discovered eight sequences with perfect identity and six sequences with close similarity to previously defined MHC class I sequences from other macaque populations. We identified two ancestral MHC haplotypes that appear to be shared between Filipino and Mauritian cynomolgus macaques, notably a Mafa-B haplotype that has previously been shown to protect Mauritian cynomolgus macaques against challenge with a simian/human immunodeficiency virus, SHIV_89.6P. We also identified a Filipino cynomolgus macaque MHC class I sequence for which the predicted protein sequence differs from Mamu-B*17 by a single amino acid. This is important because Mamu-B*17 is strongly associated with protection against simian immunodeficiency virus (SIV) challenge in Indian rhesus macaques. These findings have implications for the evolutionary history of Filipino cynomolgus macaques as well as for the use of this model in SIV/SHIV research protocols. Kevin J. Campbell and Ann M. Detmer contributed equally to this work.