Genome sequencing and application of Taiwanese macaque Macaca cyclopis.

Formosan macaque (Macaca cyclopis) is the only non-human primate in Taiwan Island. We performed de novo hybrid assembly for M. cyclopis using Illumina paired-end short reads, mate-pair reads and Nanopore long reads and obtained 5065 contigs with a N50 of 2.66 megabases. M. cyclopis contigs > = 10 kb were assigned to chromosomes using Indian rhesus macaque (Macaca mulatta mulatta) genome assembly Mmul_10 as reference, resulting in a draft of M. cyclopis genome of 2,846,042,475 bases, distributed in 21 chromosomes. The draft genome contains 23,462 transcriptional origins (genes), capable of expressing 716,231 exons in 59,484 transcripts. Genome-based phylogenetic study using the assembled M. cyclopis genome together with genomes of four other macaque species, human, orangutan and chimpanzee showed similar result as previously reported. However, the M. cyclopis species was found to diverge from Chinese M. mulatta lasiota about 1.8 million years ago. Fossil gene analysis detected the presence of gap and pol endogenous viral elements of simian retrovirus in all macaques tested, including M. fascicularis, M. m. mulatta and M. cyclopis. However, M. cyclopis showed ~ 2 times less in number and more uniform in chromosomal locations. The constrain in foreign genome disturbance, presumably due to geographical isolation, should be able to simplify genomics-related investigations, making M. cyclopis an ideal primate species for medical research.

Genetic differentiation between and within Northern Native American language groups: an argument for the expansion of the Native American CODIS database.

The National Research Council recommends that genetic differentiation among subgroups of ethnic samples be lower than 3% of the total genetic differentiation within the ethnic sample to be used for estimating reliable random match probabilities for forensic use. Native American samples in the United States' Combined DNA Index System (CODIS) database represent four language families: Algonquian, Na-Dene, Eskimo-Aleut, and Salishan. However, a minimum of 27 Native American language families exists in the US, not including language isolates. Our goal was to ascertain whether genetic differences are correlated with language groupings and, if so, whether additional language families would provide a more accurate representation of current genetic diversity among tribal populations. The 21 short tandem repeat (STR) loci included in the Globalfiler® PCR Amplification Kit were used to characterize six indigenous language families, including three of the four represented in the CODIS database (i.e. Algonquian, Na-Dene, and Eskimo-Aleut), and two language isolates (Miwok and Seri) using major population genetic diversity metrics such as F statistics and Bayesian clustering analysis of genotype frequencies. Most of the genetic variation (97%) was found to be within language families instead of among them (3%). In contrast, when only the three of the four language families represented in both the CODIS database and the present study were considered, 4% of the genetic variation occurred among the language groups. Bayesian clustering resulted in a maximum posterior probability indicating three genetically distinct groups among the eight language families and isolates: (1) Eskimo, (2) Seri, and (3) all other language groups and isolates, thus confirming genetic subdivision among subgroups of the CODIS Native American database. This genetic structure indicates the need for an increased number of Native American populations based on language affiliation in the CODIS database as well as more robust sample sets for those language families. Supplemental data for this article is available online at https://doi.org/10.1080/20961790.2021.1963088 .

Population genetics of the ABO locus within the rhesus (Macaca mulatta) and cynomolgus (M. fascicularis) macaque hybrid zone.

Knowledge of the macaque ABO blood group system has been critical in the development of nonhuman primates (NHPs) as a translational model. Serving not only as a useful homologue of the disease-linked ABO system in humans, macaque ABO blood groups must be typed in colonies prior to performing experimental procedures requiring blood transfusion or transplantation. While the rates of blood type incompatibility and the distributions of A, B and AB blood groups are known in large samples of rhesus (Macaca mulatta) and cynomolgus (M. fascicularis) macaques, there is a dearth of blood type data from macaque populations occupying the rhesus-cynomolgus hybrid zone in Southeast Asia. Using molecular phenotyping, we profiled ABO blood group distributions of 232 macaques from 10 populations in the hybrid zone and compared them to pure blood populations of the two species. We found that while these distributions are significantly different in most populations, there was a lack of differentiation between the hybrid and cynomolgus macaques as well as between the Thai and neighbouring populations. This supports a more expansive model of hybridization between rhesus and cynomolgus macaques than often proposed and highlights the increased need for consideration of population genetic structure in biomedical studies that employ macaques as animal models. Additionally, we report an enrichment of indeterminate blood types in the hybrid populations.

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

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

Determination of major histocompatibility class I and class II genetic composition of the Caribbean Primate Center specific pathogen-free rhesus macaque (Macaca mulatta) colony based on massively parallel sequencing.

BACKGROUND: Knowledge of major histocompatibility complex (MHC) composition and distribution in rhesus macaque colonies is critical for management strategies that maximize the utility of this model for biomedical research. METHODS: Variation within the Mamu-A and Mamu-B (class I) and DRB, DQA/B, and DPA/B (class II) regions of 379 animals from the Caribbean Primate Research Center's (CPRC) specific pathogen free (SPF) colony was examined using massively parallel sequencing. RESULTS: Analyses of the 7 MHC loci revealed a background of Indian origin with high levels of variation despite past genetic bottlenecks. All loci exhibited mutual linkage disequilibria while conforming to Hardy-Weinberg expectations suggesting the achievement of mutation-selection balance. CONCLUSION: The CPRC's SPF colony is a significant resource for research on AIDS and other infectious agents. Characterizing colony-wide MHC variability facilitates the breeding and selection of animals bearing desired haplotypes and increases the investigator's ability to understand the immune responses mounted by these animals.

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

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

Mycobacterium kansasii Isolated from Tuberculinpositive Rhesus Macaques (Macaca mulatta) in the Absence of Disease.

Mycobacterial infections are of primary health concern in NHP colonies in biomedical research. NHP are constantly monitored and screened for Mycobacterium spp. We report 6 Chinese-origin rhesus macaques infected with Mycobacterium kansasii that exhibited positive tuberculin skin tests in the absence of disease. Two of these macaques were being used for research purposes; the remaining 4 macaques were residing at the contract quarantine company. Histopathology and acid-fast staining of fixed tissues from all macaques showed that all were free of disease. Thoracic radiographs were negative for any signs of disease or infection. Samples from bronchial lavage and tissues including lung, spleen, hilar and mesenteric lymph nodes tested negative by PCR assay for Mycobacterium spp. One of the research macaques tested culture-positive for M. kansasii and a poorly characterized M. avium complex organism. One macaque from the contract quarantine facility tested culture positive for M. kansasii. Genomic testing and target gene RNA expression analysis of the 2 M. kansasii isolates were performed to evaluate possible kinship and affected genes that might contribute to susceptibility to mycobacterial infection. Genotyping of the 2 isolates revealed 2 genetically distinct strains (strains 1 and 4). The presence of positive tuberculin skin tests in the absence of disease raises serious concerns regarding diagnostic methods used for infected NHP.

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

BACKGROUND: Most cynomolgus macaques (Macaca fascicularis) used in the United States as animal models are imported from Chinese breeding farms without documented ancestry. Cynomolgus macaques with varying rhesus macaque ancestry proportions may exhibit differences, such as susceptibility to malaria, that affect their suitability as a research model. METHODS: DNA of 400 cynomolgus macaques from 10 Chinese breeding farms was genotyped to characterize their regional origin and rhesus ancestry proportion. A nested PCR assay was used to detect Plasmodium cynomolgi infection in sampled individuals. RESULTS: All populations exhibited high levels of genetic heterogeneity and low levels of inbreeding and genetic subdivision. Almost all individuals exhibited an Indochinese origin and a rhesus ancestry proportion of 5%-48%. The incidence of P. cynomolgi infection in cynomolgus macaques is strongly associated with proportion of rhesus ancestry. CONCLUSIONS: The varying amount of rhesus ancestry in cynomolgus macaques underscores the importance of monitoring their genetic similarity in malaria research.

Flanking region variation of ForenSeq™ DNA Signature Prep Kit STR and SNP loci in Yavapai Native Americans.

Massively parallel sequencing (MPS) offers advantages over current capillary electrophoresis-based analysis of short tandem repeat (STR) loci for human identification testing. In particular STR repeat motif sequence information can be obtained, thereby increasing the discrimination power of some loci. While sequence variation within the repeat region is observed relatively frequently in some of the commonly used STRs, there is an additional degree of variation found in the flanking regions adjacent to the repeat motif. Repeat motif and flanking region sequence variation have been described for major population groups, however, not for more isolated populations. Flanking region sequence variation in STR and single nucleotide polymorphism (SNP) loci in the Yavapai population was analyzed using the ForenSeq™ DNA Signature Prep Kit and STRait Razor v2s. Seven and 14 autosomal STRs and identity-informative single nucleotide polymorphisms (iiSNPs), respectively, had some degree of flanking region variation. Three and four of these identity-informative loci, respectively, showed ≥5% increase in expected heterozygosity. The combined length- and sequence-based random match probabilities (RMPs) for 27 autosomal STRs were 6.11×10-26 and 2.79×10-29, respectively. When combined with 94 iiSNPs (a subset of which became microhaplotypes) the combined RMP was 5.49×10-63. Analysis of length-based and sequence-based autosomal STRs in STRUCTURE indicated that the Yavapai are most similar to the Hispanic population. While producing minimal increase in X- and Y-STR discrimination potential, access to flanking region data enabled identification of one novel X-STR and three Y-STR alleles relative to previous reports. Five ancestry-informative SNPs (aiSNPs) and two phenotype-informative SNPs (piSNPs) exhibited notable flanking region variation.

Mitochondrial DNA Analysis of Mazahua and Otomi Indigenous Populations from Estado de México Suggests a Distant Common Ancestry.

The indigenous Mazahua and Otomi have inhabited the same localities in Estado de México since pre-Columbian times. Their languages, Mazahua and Otomi, belong to the Oto-Manguean linguistic family, and although they share cultural traditions and a regional history that suggest close genetic relationships and common ancestry, the historical records concerning their origin are confusing. To understand the biological relationships between Mazahua and Otomi, we analyzed mitochondrial DNA (mtDNA) genetic variation. We identified the mtDNA haplogroups by restriction fragment length polymorphism typing and sequenced hypervariable region 1 of the mtDNA control region in 141 Mazahua and 100 Otomi. These results showed that Otomi exhibit a higher frequency of haplogroup A than B, whereas Mazahua exhibit the opposite pattern. In the Otomi EM population the most frequent subhaplogroups are, in order of frequency, A2, B2, and C1, whereas in the Mazahua 1 population they are B2, D1, and A2. The most frequent haplotypes (Ht) of haplogroups A and B are Ht2 (A) and Ht58 (B2g1) in Mazahua 1 and Ht8 (A2), Ht22 (A2ao1), and Ht53 (B2c2b) in Otomi EM. The genetic differences between the Mazahua 1 and Otomi EM suggest a distant shared ancestry and a moderate degree of maternal admixture that has not obscured the difference of their mtDNA patterns. These unexpected results suggest the Mazahua and Otomi probably descend from the same group but separated very early and admixed with other Mesoamerican populations before their arrival in Central Mexico. The historical evidence of conflicting relations between the Mazahua and Otomi and the almost nonexistence of marriage between them could be responsible for maintaining only a moderate degree of maternal admixture.

Distribution and prevalence of malaria parasites among long-tailed macaques (Macaca fascicularis) in regional populations across Southeast Asia.

BACKGROUND: Plasmodium knowlesi and Plasmodium cynomolgi are two malaria parasites naturally transmissible between humans and wild macaque through mosquito vectors, while Plasmodium inui can be experimentally transmitted from macaques to humans. One of their major natural hosts, the long-tailed macaque (Macaca fascicularis), is host to two other species of Plasmodium (Plasmodium fieldi and Plasmodium coatneyi) and is widely distributed in Southeast Asia. This study aims to determine the distribution of wild macaques infected with malarial parasites by examining samples derived from seven populations in five countries across Southeast Asia. METHODS: Plasmodium knowlesi, P. cynomolgi, P. coatneyi, P. inui and P. fieldi, were detected using nested PCR assays in DNA samples from 276 wild-caught long-tailed macaques. These samples had been derived from macaques captured at seven locations, two each in the Philippines (n = 68) and Indonesia (n = 70), and one each in Cambodia (n = 54), Singapore (n = 40) and Laos (n = 44). The results were compared with previous studies of malaria parasites in long-tailed macaques from other locations in Southeast Asia. Fisher exact test and Chi square test were used to examine the geographic bias of the distribution of Plasmodium species in the macaque populations. RESULTS: Out of 276 samples tested, 177 were Plasmodium-positive, with P. cynomolgi being the most common and widely distributed among all long-tailed macaque populations (53.3 %) and occurring in all populations examined, followed by P. coatneyi (20.4 %), P. inui (12.3 %), P. fieldi (3.4 %) and P. knowlesi (0.4 %). One P. knowlesi infection was detected in a macaque from Laos, representing the first documented case of P. knowlesi in wildlife in Laos. Chi square test showed three of the five parasites (P. knowlesi, P. coatneyi, P. cynomolgi) with significant bias in prevalence towards macaques from Malaysian Borneo, Cambodia, and Southern Sumatra, respectively. CONCLUSIONS: The prevalence of malaria parasites, including those that are transmissible to humans, varied among all sampled regional populations of long-tailed macaques in Southeast Asia. The new discovery of P. knowlesi infection in Laos, and the high prevalence of P. cynomolgi infections in wild macaques in general, indicate the strong need of public advocacy in related countries.

Native American population data based on the Globalfiler(®) autosomal STR loci.

Native American population data are limited and thus impact computing accurate statistical parameters for forensic investigations. Thus, additional information should be generated from geographically representative tribes in North America, particularly from those that are not included in existing population databases for forensic use. The Globafiler(®) PCR Amplification kit was used to produce STR genotypic data for 533 individuals who represent 31 Native American tribal populations derived from eight geographically diverse regions in North America. Population genetic estimates from 21 autosomal STRs are reported.

Genetic analysis of the Yavapai Native Americans from West-Central Arizona using the Illumina MiSeq FGx™ forensic genomics system.

Forensically-relevant genetic markers were typed for sixty-two Yavapai Native Americans using the ForenSeq™ DNA Signature Prep Kit.These data are invaluable to the human identity community due to the greater genetic differentiation among Native American tribes than among other subdivisions within major populations of the United States. Autosomal, X-chromosomal, and Y-chromosomal short tandem repeat (STR) and identity-informative (iSNPs), ancestry-informative (aSNPs), and phenotype-informative (pSNPs) single nucleotide polymorphism (SNP) allele frequencies are reported. Sequence-based allelic variants were observed in 13 autosomal, 3 X, and 3 Y STRs. These observations increased observed and expected heterozygosities for autosomal STRs by 0.081±0.068 and 0.073±0.063, respectively, and decreased single-locus random match probabilities by 0.051±0.043 for 13 autosomal STRs. The autosomal random match probabilities (RMPs) were 2.37×10-26 and 2.81×10-29 for length-based and sequence-based alleles, respectively. There were 22 and 25 unique Y-STR haplotypes among 26 males, generating haplotype diversities of 0.95 and 0.96, for length-based and sequencebased alleles, respectively. Of the 26 haplotypes generated, 17 were assigned to haplogroup Q, three to haplogroup R1b, two each to haplogroups E1b1b and L, and one each to haplogroups R1a and I1. Male and female sequence-based X-STR random match probabilities were 3.28×10-7 and 1.22×10-6, respectively. The average observed and expected heterozygosities for 94 iSNPs were 0.39±0.12 and 0.39±0.13, respectively, and the combined iSNP RMP was 1.08×10-32. The combined STR and iSNP RMPs were 2.55×10-58 and 3.02×10-61 for length-based and sequence-based STR alleles, respectively. Ancestry and phenotypic SNP information, performed using the ForenSeq™ Universal Analysis Software, predicted black hair, brown eyes, and some probability of East Asian ancestry for all but one sample that clustered between European and Admixed American ancestry on a principal components analysis. These data serve as the first population assessment using the ForenSeq™ panel and highlight the value of employing sequence-based alleles for forensic DNA typing to increase heterozygosity, which is beneficial for identity testing in populations with reduced genetic diversity.

Performing monkeys of Bangladesh: characterizing their source and genetic variation.

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

Genetic Structure of First Nation Communities in the Pacific Northwest.

This study presents genetic data for nine Native American populations from northern North America. Analyses of genetic variation focus on the Pacific Northwest (PNW). Using mitochondrial, Y chromosomal, and autosomal DNA variants, we aimed to more closely address the relationships of geography and language with present genetic diversity among the regional PNW Native American populations. Patterns of genetic diversity exhibited by the three genetic systems were consistent with our hypotheses: genetic variation was more strongly explained by geographic proximity than by linguistic structure. Our findings were corroborated through a variety on analytic approaches, with the unrooted trees for the three genetic systems consistently separating inland from coastal PNW populations. Furthermore, analyses of molecular variance support the trends exhibited by the unrooted trees, with geographic partitioning of PNW populations (FCT = 19.43%, p = 0.010 ± 0.009) accounting for over twice as much of the observed genetic variation as linguistic partitioning of the same populations (FCT = 9.15%, p = 0.193 ± 0.013). These findings demonstrate a consensus with previous PNW population studies examining the relationships of genome-wide variation, mitochondrial haplogroup frequencies, and skeletal morphology with geography and language.

Identifying rhesus macaque gene orthologs using heterospecific human CNV probes.

We used the Affymetrix(®) Genome-Wide Human SNP Array 6.0 to identify heterospecific markers and compare copy number and structural genomic variation between humans and rhesus macaques. Over 200,000 human copy number variation (CNV) probes were mapped to a Chinese and an Indian rhesus macaque sample. Observed genomic rearrangements and synteny were in agreement with the results of a previously published genomic comparison between humans and rhesus macaques. Comparisons between each of the two rhesus macaques and humans yielded 206 regions with copy numbers that differed by at least two fold in the Indian rhesus macaque and human, 32 in the Chinese rhesus macaque and human, and 147 in both rhesus macaques. The detailed genomic map and preliminary CNV data are useful for better understanding genetic variation in rhesus macaques, identifying derived changes in human CNVs that may have evolved by selection, and determining the suitability of rhesus macaques as human models for particular biomedical studies.

Gene expression and TB pathogenesis in rhesus macaques: TR4, CD40, CD40L, FAS (CD95), and TNF are host genetic markers in peripheral blood mononuclear cells that are associated with severity of TB lesions.

Tuberculosis (TB) pathologic lesions in rhesus macaques resemble those in humans. The expression levels of several host TB candidate genes in the peripheral blood mononuclear cells (PBMCs) of six rhesus macaques experimentally infected with Mycobacterium tuberculosis were quantified pre-infection and at several dates post-infection. Quantitative measures of TB histopathology in the lungs including: granuloma count, granuloma size, volume of granulomatous and non-granulomatous lesions, and direct bacterial load, were used as the outcomes of a multi-level Bayesian regression model in which expression levels of host genes at various dates were used as predictors. The results indicate that the expression levels of TR4, CD40, CD40L, FAS (CD95) and TNF in PBMC were associated with quantitative measures of the severity of TB histopathologic lesions in the lungs of the study animals. Moreover, no reliable association between the expression levels of IFNE in PBMCs and the severity of TB lesions in the lungs of the study animals was found. In conclusion, PBMC expression profiles derived from the above-listed host genes might be appropriate biomarkers for probabilistic diagnosis and/or prognosis of TB severity in rhesus macaques.

Complete Taiwanese Macaque (Macaca cyclopis) Mitochondrial Genome: Reference-Assisted de novo Assembly with Multiple k-mer Strategy.

The Taiwanese (Formosan) macaque (Macaca cyclopis) is the only nonhuman primate endemic to Taiwan. This primate species is valuable for evolutionary studies and as subjects in medical research. However, only partial fragments of the mitochondrial genome (mitogenome) of this primate species have been sequenced, not mentioning its nuclear genome. We employed next-generation sequencing to generate 2 x 90 bp paired-end reads, followed by reference-assisted de novo assembly with multiple k-mer strategy to characterize the M. cyclopis mitogenome. We compared the assembled mitogenome with that of other macaque species for phylogenetic analysis. Our results show that, the M. cyclopis mitogenome consists of 16,563 nucleotides encoding for 13 protein-coding genes, 2 ribosomal RNAs and 22 transfer RNAs. Phylogenetic analysis indicates that M. cyclopis is most closely related to M. mulatta lasiota (Chinese rhesus macaque), supporting the notion of Asia-continental origin of M. cyclopis proposed in previous studies based on partial mitochondrial sequences. Our work presents a novel approach for assembling a mitogenome that utilizes the capabilities of de novo genome assembly with assistance of a reference genome. The availability of the complete Taiwanese macaque mitogenome will facilitate the study of primate evolution and the characterization of genetic variations for the potential usage of this species as a non-human primate model for medical research.

Heterospecific SNP diversity in humans and rhesus macaque (Macaca mulatta).

BACKGROUND: Conservation of single nucleotide polymorphisms (SNPs) between human and other primates (i.e., heterospecific SNPs) in candidate genes can be used to assess the utility of those organisms as models for human biomedical research. METHODS: A total of 59,691 heterospecific SNPs in 22 rhesus macaques and 20 humans were analyzed for human trait associations and 4207 heterospecific SNPs biallelic in both taxa were compared for genetic variation. RESULTS: Variation comparisons at the 4207 SNPs showed that humans were more genetically diverse than rhesus macaques with observed and expected heterozygosities of 0.337 and 0.323 vs. 0.119 and 0.102, and minor allele frequencies of 0.239 and 0.063, respectively. In total, 431 of the 59,691 heterospecific SNPs are reportedly associated with human-specific traits. CONCLUSION: While comparisons between human and rhesus macaque genomes are plausible, functional studies of heterospecific SNPs are necessary to determine whether rhesus macaque alleles are associated with the same phenotypes as their corresponding human alleles.

Evolutionary distance of amino acid sequence orthologs across macaque subspecies: identifying candidate genes for SIV resistance in Chinese rhesus macaques.

We use the Reciprocal Smallest Distance (RSD) algorithm to identify amino acid sequence orthologs in the Chinese and Indian rhesus macaque draft sequences and estimate the evolutionary distance between such orthologs. We then use GOanna to map gene function annotations and human gene identifiers to the rhesus macaque amino acid sequences. We conclude methodologically by cross-tabulating a list of amino acid orthologs with large divergence scores with a list of genes known to be involved in SIV or HIV pathogenesis. We find that many of the amino acid sequences with large evolutionary divergence scores, as calculated by the RSD algorithm, have been shown to be related to HIV pathogenesis in previous laboratory studies. Four of the strongest candidate genes for SIVmac resistance in Chinese rhesus macaques identified in this study are CDK9, CXCL12, TRIM21, and TRIM32. Additionally, ANKRD30A, CTSZ, GORASP2, GTF2H1, IL13RA1, MUC16, NMDAR1, Notch1, NT5M, PDCD5, RAD50, and TM9SF2 were identified as possible candidates, among others. We failed to find many laboratory experiments contrasting the effects of Indian and Chinese orthologs at these sites on SIVmac pathogenesis, but future comparative studies might hold fertile ground for research into the biological mechanisms underlying innate resistance to SIVmac in Chinese rhesus macaques.