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.

Validation of a suite of ERP and QEEG biomarkers in a pre-competitive, industry-led study in subjects with schizophrenia and healthy volunteers

OBJECTIVE

Complexity and lack of standardization have mostly limited the use of event-related potentials (ERPs) and quantitative EEG (QEEG) biomarkers in drug development to small early phase trials. We present results from a clinical study on healthy volunteers (HV) and patients with schizophrenia (SZ) that assessed test-retest, group differences, variance, and correlation with functional assessments for ERP and QEEG measures collected at clinical and commercial trial sites with standardized instrumentation and methods, and analyzed through an automated data analysis pipeline.

METHODS

81 HV and 80 SZ were tested at one of four study sites. Subjects were administered two ERP/EEG testing sessions on separate visits. Sessions included a mismatch negativity paradigm, a 40 Hz auditory steady-state response paradigm, an eyes-closed resting state EEG, and an active auditory oddball paradigm. SZ subjects were also tested on the Brief Assessment of Cognition (BAC), Positive and Negative Syndrome Scale (PANSS), and Virtual Reality Functional Capacity Assessment Tool (VRFCAT).

RESULTS

Standardized ERP/EEG instrumentation and methods ensured few test failures. The automated data analysis pipeline allowed for near real-time analysis with no human intervention. Test-retest reliability was fair-to-excellent for most of the outcome measures. SZ subjects showed significant deficits in ERP and QEEG measures consistent with published academic literature. A subset of ERP and QEEG measures correlated with functional assessments administered to the SZ subjects.

CONCLUSIONS

With standardized instrumentation and methods, complex ERP/EEG testing sessions can be reliably performed at clinical and commercial trial sites to produce high-quality data in near real-time.

Contemplating Health Economics, Coding and Reimbursement in Orthotics, Prosthetics and Pedorthics

Reimbursement to U.S. healthcare service providers is largely transitioning from fee for service to fee for value for those clinicians who code using current procedural terminology and through their coding, describe their professional services. The Orthotic, Prosthetic and Pedorthic profession (O&P), currently codes using a system that describes the devices they evaluate for, fabricate, fit and maintain and their professional services are incorporated into their codes. These O&P codes, in contrast to those for other healthcare disciplines, are predominantly product based rather than service based, focusing on product features and function more than clinical service. This editorial manuscript provides a brief overview of the system the US O&P profession uses currently, particularly in the context of other healthcare professions transitioning to value based coding and reimbursement and culminates in a call to action for the profession to academically consider the strengths and weaknesses of the current system relative to alternative systems.

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 .

Effects of chronic caffeine exposure during adolescence and subsequent acute caffeine challenge during adulthood on rat brain serotonergic systems

Caffeine is the most commonly used drug in the world. However, animal studies suggest that chronic consumption of caffeine during adolescence can result in enhanced anxiety-like behavioral responses during adulthood. One mechanism through which chronic caffeine administration may influence subsequent anxiety-like responses is through actions on brainstem serotonergic systems. In order to explore potential effects of chronic caffeine consumption on brainstem serotonergic systems, we evaluated the effects of a 28-day exposure to chronic caffeine (0.3 g/L; postnatal day 28-56) or vehicle administration in the drinking water, followed by 24 h caffeine withdrawal, and subsequent challenge with caffeine (30 mg/kg; s.c.) or vehicle in adolescent male rats. In Experiment 1, acute caffeine challenge induced a widespread activation of serotonergic neurons throughout the dorsal raphe nucleus (DR); this effect was attenuated in rats that had been exposed to chronic caffeine consumption. In Experiment 2, acute caffeine administration profoundly decreased tph2 and slc22a3 mRNA expression throughout the DR, with no effects on htr1a or slc6a4 mRNA expression. Chronic caffeine exposure for four weeks during adolescence was sufficient to decrease tph2 mRNA expression in the DR measured 28 h after caffeine withdrawal. Chronic caffeine administration during adolescence did not impact the ability of acute caffeine to decrease tph2 or slc22a3 mRNA expression. Together, these data suggest that both chronic caffeine administration during adolescence and acute caffeine challenge during adulthood are important determinants of serotonergic function and serotonergic gene expression, effects that may contribute to chronic effects of caffeine on anxiety-like responses.

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.

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.

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.

ABO blood group phenotype frequency estimation using molecular phenotyping in rhesus and cynomolgus macaques

A much larger sample (N = 2369) was used to evaluate a previously reported distribution of the A, AB and B blood group phenotypes in rhesus and cynomolgus macaques from six different regional populations. These samples, acquired from 15 different breeding and research facilities in the United States, were analyzed using a real-time quantitative polymerase chain reaction (qPCR) assay that targets single nucleotide polymorphisms (SNPs) responsible for the macaque A, B and AB phenotypes. The frequency distributions of blood group phenotypes of the two species differ significantly from each other and significant regional differentiation within the geographic ranges of each species was also observed. The B blood group phenotype was prevalent in rhesus macaques, especially those from India, while the frequencies of the A, B and AB phenotypes varied significantly among cynomolgus macaques from different geographic regions. The Mauritian cynomolgus macaques, despite having originated in Indonesia, showed significant (P ≪ .01) divergence from the Indonesian animals at the ABO blood group locus. Most Mauritian animals belonged to the B blood group while the Indonesian animals were mostly A. The close similarity in blood group frequency distributions between the Chinese rhesus and Indochinese cynomolgus macaques demonstrates that the introgression between these two species extends beyond the zone of intergradation in Indochina. This study underscores the importance of ABO blood group phenotyping of the domestic supply of macaques and their biospecimens.

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.

Effects of naltrexone are influenced by childhood adversity during negative emotional processing in addiction recovery

Naltrexone is an opioid receptor antagonist used in the management of alcohol dependence. Although the endogenous opioid system has been implicated in emotion regulation, the effects of mu-opioid receptor blockade on brain systems underlying negative emotional processing are not clear in addiction. Individuals meeting criteria for alcohol dependence alone (n=18, alcohol) and in combination with cocaine and/or opioid dependence (n=21, alcohol/drugs) and healthy individuals without a history of alcohol or drug dependence (n=21) were recruited. Participants were alcohol and drug abstinent before entered into this double-blind, placebo-controlled, randomized, crossover study. Functional magnetic resonance imaging was used to investigate brain response while viewing aversive and neutral images relative to baseline on 50 mg of naltrexone and placebo. We found that naltrexone modulated task-related activation in the medial prefrontal cortex and functional connectivity between the anterior cingulate cortex and the hippocampus as a function of childhood adversity (for aversive versus neutral images) in all groups. Furthermore, there was a group-by-treatment-by-condition interaction in the right amygdala, which was mainly driven by a normalization of response for aversive relative to neutral images under naltrexone in the alcohol/drugs group. We conclude that early childhood adversity is one environmental factor that influences pharmacological response to naltrexone. Pharmacotherapy with naltrexone may also have some ameliorative effects on negative emotional processing in combined alcohol and drug dependence, possibly due to alterations in endogenous opioid transmission or the kappa-opioid receptor antagonist actions of naltrexone.

The genetic structure of native Americans in North America based on the Globalfiler® STRs

Current forensic STR databases, such as CODIS, lack population genetic data on Native American populations. Information from a geographically diverse array of tribes is necessary to provide improved statistical estimates of the strength of associations with DNA evidence. The Globalfiler® STR markers were used to characterize the genetic structure of ten tribal populations from seven geographic regions in North America, including those not presently represented in forensic databases. Samples from the Arctic region, Baja California, California/Great Basin, the Southeast, Mexico, the Midwest, and the Southwest were analyzed for allele frequencies, observed and expected heterozygosities, and F-statistics. The tribal samples exhibited an F_ST or θ value above the conservative 0.03 estimate recommended by the National Research Council (NRC) for calculating random match probabilities among Native Americans. The greater differentiation among tribal populations computed here (θ=0.04) warrants the inclusion of additional regional Native American samples into STR databases.

The Microbiome of the Built Environment and Human Behavior: Implications for Emotional Health and Well-Being in Postmodern Western Societies

It is increasingly evident that inflammation is an important determinant of cognitive function and emotional behaviors that are dysregulated in stress-related psychiatric disorders, such as anxiety and affective disorders. Inflammatory responses to physical or psychological stressors are dependent on immunoregulation, which is indicated by a balanced expansion of effector T-cell populations and regulatory T cells. This balance is in part driven by microbial signals. The hygiene or "old friends" hypothesis posits that exposure to immunoregulation-inducing microorganisms is reduced in modern urban societies, leading to an epidemic of inflammatory disease and increased vulnerability to stress-related psychiatric disorders. With the global trend toward urbanization, humans are progressively spending more time in built environments, thereby, experiencing limited exposures to these immunoregulatory "old friends." Here, we evaluate the implications of the global trend toward urbanization, and how this transition may affect human microbial exposures and human behavior.

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.

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.

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.

POPULATION GENETICS. Genomic evidence for the Pleistocene and recent population history of Native Americans

How and when the Americas were populated remains contentious. Using ancient and modern genome-wide data, we found that the ancestors of all present-day Native Americans, including Athabascans and Amerindians, entered the Americas as a single migration wave from Siberia no earlier than 23 thousand years ago (ka) and after no more than an 8000-year isolation period in Beringia. After their arrival to the Americas, ancestral Native Americans diversified into two basal genetic branches around 13 ka, one that is now dispersed across North and South America and the other restricted to North America. Subsequent gene flow resulted in some Native Americans sharing ancestry with present-day East Asians (including Siberians) and, more distantly, Australo-Melanesians. Putative "Paleoamerican" relict populations, including the historical Mexican Pericúes and South American Fuego-Patagonians, are not directly related to modern Australo-Melanesians as suggested by the Paleoamerican Model.

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.

Population genetics of the California National Primate Research Center's (CNPRC) captive Callicebus cupreus colony

The California National Primate Research Center maintains a small colony of titi monkeys (Callicebus cupreus) for behavioral studies. While short tandem repeat (STR) markers are critical for the genetic management of the center's rhesus macaque (Macaca mulatta) breeding colony, STRs are not used for this purpose in the maintenance of the center's titi monkey colony. Consequently, the genetic structure of this titi monkey population has not been characterized. A lack of highly informative genetic markers in titi monkeys has also resulted in scant knowledge of the species' genetic variation in the wild. The purpose of this study was to develop a panel of highly polymorphic titi monkey STRs using a cross-species polymerase chain reaction (PCR) amplification protocol that could be used for the genetic management of the titi monkey colony. We screened 16 STR primer pairs and selected those that generated robust and reproducible polymorphic amplicons. Loci that were found to be highly polymorphic, very likely to be useful for parentage verification, pedigree assessment, and studying titi monkey population genetics, were validated using Hardy-Weinberg equilibrium and linkage disequilibrium analyses. The genetic data generated in this study were also used to assess directly the impact on the colony's genetic diversity of a recent adenovirus outbreak. While the adenovirus epizootic disease caused significant mortality (19 deaths among the 65 colony animals), our results suggest that the disease exhibited little or no influence on the overall genetic diversity of the colony.

A genetic comparison of two alleged subspecies of Philippine cynomolgus macaques

Two subspecies of cynomolgus macaques (Macaca fascicularis) are alleged to co-exist in the Philippines, M. f. philippensis in the north and M. f. fascicularis in the south. However, genetic differences between the cynomolgus macaques in the two regions have never been studied to document the propriety of their subspecies status. We genotyped samples of cynomolgus macaques from Batangas in southwestern Luzon and Zamboanga in southwestern Mindanao for 15 short tandem repeat (STR) loci and sequenced an 835 bp fragment of the mtDNA of these animals. The STR genotypes were compared with those of cynomolgus macaques from southern Sumatra, Singapore, Mauritius and Cambodia, and the mtDNA sequences of both Philippine populations were compared with those of cynomolgus macaques from southern Sumatra, Indonesia and Sarawak, Malaysia. We conducted STRUCTURE and PCA analyses based on the STRs and constructed a median joining network based on the mtDNA sequences. The Philippine population from Batangas exhibited much less genetic diversity and greater genetic divergence from all other populations, including the Philippine population from Zamboanga. Sequences from both Batangas and Zamboanga were most closely related to two different mtDNA haplotypes from Sarawak from which they are apparently derived. Those from Zamboanga were more recently derived than those from Batangas, consistent with their later arrival in the Philippines. However, clustering analyses do not support a sufficient genetic distinction of cynomolgus macaques from Batangas from other regional populations assigned to subspecies M. f. fascicularis to warrant the subspecies distinction M. f. philippensis.

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

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

Expression levels of 10 candidate genes in lung tissue of vaccinated and TB-infected cynomolgus macaques

The expression of ten tuberculosis candidate genes in lung and lymph nodes of cynomolgus macaques vaccinated and experimentally infected with Mycobacterium tuberculosis (Mtb) was quantified. The expression of TNFα, IL10, IL1β, TLR4, IL17, IL6, IL12, and iNOS in the lungs of vaccinated animals was higher than that of non-vaccinated animals.

An empirical comparison of short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) for relatedness estimation in Chinese rhesus macaques (Macaca mulatta)

We compare the effectiveness of short tandem repeat (STR) and single nucleotide polymorphism (SNP) genotypes for estimating pairwise relatedness, using molecular data and pedigree records from a captive Chinese rhesus macaque population at the California National Primate Research Center. We find that a panel of 81 SNPs is as effective at estimating first-order kin relationships as a panel of 14 highly polymorphic STRs. We note, however, that the selected STRs provide more precise predictions of relatedness than the selected SNPs, and may be preferred in contexts that require the discrimination of kin related more distantly than first-order relatives. Additionally, we compare the performance of three commonly used relatedness estimation algorithms, and find that the Wang [2002] algorithm outperforms other algorithms when analyzing STR data, while the Queller & Goodnight [1989] algorithm outperforms other algorithms when analyzing SNP data. Future research is needed to address the number of SNPs required to reach the discriminatory power of a standard STR panel in relatedness estimation for primate colony management.

Familial aggregation of chronic diarrhea disease (CDD) in rhesus macaques (Macaca mulatta)

Chronic diarrheal disease (CDD) is a critical problem for breeders of captive rhesus macaque (Macaca mulatta), as it results in significant levels of morbidity and death annually. As with other inflammatory disorders, CDD is thought to be caused by environmental and/or genetic factors. Although correspondence between the characters defined as Mendelian by pedigree or segregation analysis and functional genes is difficult to establish, such analyses provide essential entry points into understanding CDD in captive bred rhesus macaques. To investigate the familial aggregation of CDD in captive rhesus macaque, we performed pedigree, segregation and heritability analyses on genealogical data from 55 severely affected individuals (probands) through whom relatives with a history of CDD were ascertained from routine computerized colony records comprising vital and demographic statistics of 10,814 rhesus macaques. We identified 175 rhesus macaques with CDD and estimated its incidence as approximately 2% in the colony. The disease strongly clustered in eight multi-generation pedigrees. Inspection of the pedigrees, segregation analysis and heritability estimate of CDD suggest that susceptibility to the disease is under strong genetic control. Identification of the locations of susceptibility genes in the rhesus macaque genome could facilitate the reduction of their frequency in captive breeding facilities.

The genetic composition of populations of cynomolgus macaques (Macaca fascicularis) used in biomedical research

BACKGROUND

The genetic composition of cynomolgus macaques used in biomedical research is not as well-characterized as that of rhesus macaques.

METHODS

Populations of cynomolgus macaques from Sumatra, Corregidor, Mauritius, Singapore, Cambodia, and Zamboanga were analyzed using 24 STRs.

RESULTS

The Sumatran and Cambodian populations exhibited the highest allelic diversity, while the Mauritian population exhibited the lowest. Sumatran cynomolgus macaques were the most genetically similar to all others, consistent with an Indonesian origin of the species. The high diversity among Cambodian animals may result from interbreeding with rhesus macaques. The Philippine and Mauritian samples were the most divergent from other populations, the former due to separation from the Sunda Shelf by deepwater and the latter due to anthropogenic translocation and extreme founder effects.

CONCLUSIONS

Investigators should verify their research subjects' origin, ancestry, and pedigree to minimize risks to biomedical experimentation from genetic variance stemming from close kinship and mixed ancestry as these can obscure treatment effects.

Population and landscape genetics of an introduced species (M. fascicularis) on the island of Mauritius

The cynomolgus macaque, Macaca fascicularis, was introduced onto the island of Mauritius in the early 17(th) century. The species experienced explosive population growth, and currently exists at high population densities. Anecdotes collected from nonhuman primate trappers on the island of Mauritius allege that animals from the northern portion of the island are larger in body size than and superior in condition to their conspecifics in the south. Although previous genetic studies have reported Mauritian cynomolgus macaques to be panmictic, the individuals included in these studies were either from the southern/central or an unknown portion of the island. In this study, we sampled individuals broadly throughout the entire island of Mauritius and used spatial principle component analysis to measure the fine-scale correlation between geographic and genetic distance in this population. A stronger correlation between geographic and genetic distance was found among animals in the north than in those in the southern and central portions of the island. We found no difference in body weight between the two groups, despite anecdotal evidence to the contrary. We hypothesize that the increased genetic structure among populations in the north is related to a reduction in dispersal distance brought about by human habitation and tourist infrastructure, but too recent to have produced true genetic differentiation.

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.

Identifying human-rhesus macaque gene orthologs using heterospecific SNP probes

We genotyped a Chinese and an Indian-origin rhesus macaque using the Affymetrix Genome-Wide Human SNP Array 6.0 and cataloged 85,473 uniquely mapping heterospecific SNPs. These SNPs were assigned to rhesus chromosomes according to their probe sequence alignments as displayed in the human and rhesus reference sequences. The conserved gene order (synteny) revealed by heterospecific SNP maps is in concordance with that of the published human and rhesus macaque genomes. Using these SNPs' original human rs numbers, we identified 12,328 genes annotated in humans that are associated with these SNPs, 3674 of which were found in at least one of the two rhesus macaques studied. Due to their density, the heterospecific SNPs allow fine-grained comparisons, including approximate boundaries of intra- and extra-chromosomal rearrangements involving gene orthologs, which can be used to distinguish rhesus macaque chromosomes from human chromosomes.

Molecular ABO phenotyping in cynomolgus macaques using real-time quantitative PCR

Macaques are commonly used in biomedical research as animal models of human disease. The ABO phenotype of donors and recipients plays an important role in the success of transplantation and stem cell research of both human and macaque tissue. Traditional serological methods for ABO phenotyping can be time consuming, provide ambiguous results and/or require tissue that is unavailable or unsuitable. We developed a novel method to detect the A, B, and AB phenotypes of macaques using real-time quantitative polymerase chain reaction. This method enables the simple and rapid screening of these phenotypes in macaques without the need for fresh blood or saliva. This study reports the distribution of the A, B, and AB phenotypes of captive cynomolgus macaques that, while regionally variable, closely resembles that of rhesus macaques. Blood group B, as in rhesus macaques, predominates in cynomolgus macaques and its frequency distribution leads to a probability of major incompatibility of 41%. No silencing mutations have been identified in exon 6 or 7 in macaques that could be responsible for the O phenotype, that, although rare, have been reported. The excess homozygosity of rhesus and cynomolgus macaque genotypes in this study, that assumes the absence of the O allele, suggests the possibility of some mechanism preventing the expression of the A and B transferases.

Population genetics of the Washington National Primate Research Center's (WaNPRC) captive pigtailed macaque (Macaca nemestrina) population

Pigtailed macaques (Macaca nemestrina) provide an important model for biomedical research on human disease and for studying the evolution of primate behavior. The genetic structure of captive populations of pigtailed macaques is not as well described as that of captive rhesus (M. mulatta) or cynomolgus (M. fascicularis) macaques. The Washington National Primate Research Center houses the largest captive colony of pigtailed macaques located in several different housing facilities. Based on genotypes of 18 microsatellite (short tandem repeat [STR]) loci, these pigtailed macaques are more genetically diverse than captive rhesus macaques and exhibit relatively low levels of inbreeding. Colony genetic management facilitates the maintenance of genetic variability without compromising production goals of a breeding facility. The periodic introduction of new founders from specific sources to separate housing facilities at different times influenced the colony's genetic structure over time and space markedly but did not alter its genetic diversity significantly. Changes in genetic structure over time were predominantly due to the inclusion of animals from the Yerkes National Primate Research Center in the original colony and after 2005. Strategies to equalize founder representation in the colony have maximized the representation of the founders' genomes in the extant population. Were exchange of animals among the facilities increased, further differentiation could be avoided. The use of highly differentiated animals may confound interpretations of phenotypic differences due to the inflation of the genetic contribution to phenotypic variance of heritable traits.

4040 SNPs for genomic analysis in the rhesus macaque (Macaca mulatta)

Although the rhesus macaque (Macaca mulatta) is commonly used for biomedical research and becoming a preferred model for translational medicine, quantification of genome-wide variation has been slow to follow the publication of the genome in 2007. Here we report the properties of 4040 single nucleotide polymorphisms discovered and validated in Chinese and Indian rhesus macaques from captive breeding colonies in the United States. Frequency-matched measures of linkage disequilibrium were much greater in the Indian sample. Although the majority of polymorphisms were shared between the two populations, rare alleles were over twice as common in the Chinese sample. Indian rhesus had higher rates of heterozygosity, as well as previously undetected substructure, potentially due to admixture from Burma in wild populations and demographic events post-captivity.

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.

A rapid quantitative real-time PCR-based DNA quantification assay coupled with species--assignment capabilities for two hybridizing Macaca species

Regional populations of rhesus and long-tailed macaques exhibit fundamental differences in mitochondrial DNA, short tandem repeat and single nucleotide polymorphism variation between mainland and insular Southeast Asian populations. Some studies have revealed genetic admixture between these species due to natural hybridization and human-assisted intercrosses. A quantitative real-time PCR (qPCR) assay was developed to efficiently determine the species of origin of a macaque biological sample, and to quantify the species-specific template DNA. Prior knowledge of species identity and DNA concentrations are crucial for maintaining cost-effective methods and accurate DNA analysis. DNA from 109 regionally representative rhesus and long-tailed macaques was qPCR amplified to determine the species and template quantities. Of the 19 Vietnamese long-tailed macaques, 3 samples were discovered to be hybrids.

Genotyping single nucleotide polymorphisms (SNPs) across species in Old World Monkeys

The development of DNA markers is becoming increasingly useful in the field of primatology for studies on paternity, population history, and biomedical research. In this study, we determine the efficacy of using cross-species amplification to identify single nucleotide polymorphisms (SNPs) in closely related species. The DNA of 93 individuals representing seven Old World Monkey species was analyzed to identify SNPs using cross-species amplification and genotyping. The loci genotyped were 653 SNPs identified and validated in rhesus macaques. Of the 653 loci analyzed, 27% were estimated to be polymorphic in the samples studied. SNPs identified at the same locus among different species (coincident SNPs) were found in six of the seven species studied with longtail macaques exhibiting the highest number of coincident SNPs (84). The distribution of coincident SNPs among species is not biased based on proximity to genes in the samples studied. In addition, the frequency of coincident SNPs is not consistent with expectations based on their phylogenetic relationships. This study demonstrates that cross-species amplification and genotyping using the Illumina Golden Gate Array is a useful method to identify a large number of SNPs in closely related species, although issues with ascertainment bias may limit the type of studies where this method can be applied.

Association between hybrid status and reproductive success of captive male and female rhesus macaques (Macaca mulatta) at the California National Primate Research Center (CNPRC)

The California National Primate Research Center (CNPRC) houses more than 1,000 rhesus macaques (Macaca mulatta) of mixed Chinese-Indian ancestry. Most of these animals are kept in outdoor field cages, the colony's long term breeding resource. Since 2001, hybrids comprised between 4 and 49% of the field cage populations, but in most cases have represented a maximum of 10% of those populations. The increasing prevalence of hybrids is partly due to management efforts to distribute genetic diversity effectively and minimize genetic subdivisions. However, other factors may also contribute to the spread of hybrids within the colony, most notably variance in socio-sexual behaviors and physical attributes. It is known that hybrids of some species exhibit heterosis, such as early maturation, that can enhance reproductive success, and anecdotal observations of mixed groups of hybrid, Indian and Chinese animals at the CNPRC suggest that hybrids are more sexually active. To determine whether hybrids experienced a reproductive advantage, a study was conducted using birth records of 5,611 offspring born in the CNPRC colony between 2003 and 2009. We found that while the degree of Chinese ancestry (DCA) appeared to influence the maturational schedule of both males and females (maturation was inversely related to proportion of Chinese ancestry), DCA had no independent effect on either male or female RS or rank. Therefore, we have found no evidence that a hybrid phenotype confers an absolute reproductive advantage in our colony.

A simple multiplex polymerase chain reaction to determine ABO blood types of rhesus macaques (Macaca mulatta)

Rhesus macaques are the most common nonhuman primate model organism used in biomedical research. Their increasingly frequent use as subjects in studies involving transplantation requires that blood and other tissue antigens of donors and recipients be compatible. We report here an easy and rapid multiplex polymerase chain reaction (PCR) to determine the ABO blood group phenotypes of rhesus macaques that can be performed with only small amounts of DNA. We phenotyped 78 individuals and found this species to exhibit the A, B and AB phenotypes in frequencies that vary by geographic region. The probability of randomly pairing rhesus macaque donors and recipients that exhibit major ABO phenotype incompatibility is approximately 0.35 and 0.45 for Indian and Chinese rhesus macaques, respectively.

The effect of SNP discovery method and sample size on estimation of population genetic data for Chinese and Indian rhesus macaques (Macaca mulatta)

This study was designed to address issues regarding sample size and marker location that have arisen from the discovery of SNPs in the genomes of poorly characterized primate species and the application of these markers to the study of primate population genetics. We predict the effect of discovery sample size on the probability of discovering both rare and common SNPs and then compare this prediction with the proportion of common and rare SNPs discovered when different numbers of individuals are sequenced. Second, we examine the effect of genomic region on estimates of common population genetic data, comparing markers from both coding and non-coding regions of the rhesus macaque genome and the population genetic data calculated from these markers, to measure the degree and direction of bias introduced by SNPs located in coding versus non-coding regions of the genome. We found that both discovery sample size and genomic region surveyed affect SNP marker attributes and population genetic estimates, even when these are calculated from an expanded data set containing more individuals than the original discovery data set. Although none of the SNP detection methods or genomic regions tested in this study was completely uninformative, these results show that each has a different kind of genetic variation that is suitable for different purposes, and each introduces specific types of bias. Given that each SNP marker has an individual evolutionary history, we calculated that the most complete and unbiased representation of the genetic diversity present in the individual can be obtained by incorporating at least 10 individuals into the discovery sample set, to ensure the discovery of both common and rare polymorphisms.

Detecting signatures of inter-regional and inter-specific hybridization among the Chinese rhesus macaque specific pathogen-free (SPF) population using single nucleotide polymorphic (SNP) markers

BACKGROUND

While rates of gene flow between rhesus and longtail macaque populations near their hybrid zone in Indochina have been quantified elsewhere, this study demonstrates that the inter-specific introgression is not limited to the Indochinese hybrid zone but is more geographically widespread.

METHODS

Twelve rhesus and longtail macaque populations were analyzed using single nucleotide polymorphic (SNP) loci.

RESULTS

There is evidence for inter-specific admixture between Chinese rhesus and mainland longtails, with implications for genetic diversity both in the Chinese super-SPF population at the California National Primate Research Center and in other primate facilities. Eastern Chinese rhesus appeared more highly derived than western Chinese rhesus, and allele sharing between longtails and Chinese rhesus was not random with regard to geographic distance, but no significant nuclear genetic differences between eastern and western Chinese rhesus were detected among the 245 genic SNPs assayed.

CONCLUSION

The implications of this inter-specific admixture for the use of Chinese rhesus and mainland longtail in biomedical research should be considered.