BACKGROUND
animal research has provided the scientific community with extraordinary advances in medicine from the development of vaccines to the prevention and treatment of diseases.
unfortunately at present 85% of novel therapeutics fail in preclinical and early phase clinical trials and of the therapies that reach late phase trials an additional 50% fall short due to an inability to demonstrate efficacy and safety
<cit> .
reasons for these shortcomings include low patient recruitment, poor study design, and ineffective use of animal models
<cit> .
coupled with soaring drug development costs including both financial commitments and in years of labor, these shortfalls necessitate a biological and economic need for fundamental changes in the bench to bedside process.
furthermore, with advances in genome sequencing technologies there is a growing awareness that animal models fall short in terms of predictive power.
a recent study comparing the genomic responses of human inflammatory diseases to mouse models, for example, suggested that mice poorly mimic the human genetic response
<cit> .
continued progress in the understanding of human disease pathologies and the development of safe and effective therapies demands a more comprehensive understanding of animals in preclinical research.
although greater numbers of rodents are used in biomedical research, non-human primates are the gold standard of animal models in preclinical research offering advantages which include greater similarities in genome organization and sequence, behavior, and physiology
<cit> .
the rhesus  and cynomolgus  macaque are two of the most commonly used non-human primate species in research laboratories, sharing ~ <dig> % of their genome with humans
<cit> .
in academic research non-human primate use is most common in the fields of microbiology , biochemistry/pharmacology, and neuroscience
<cit> .
because of similarities in physiology and the central nervous system, non-human primates, for example, are crucial in stem cell-based regenerative medicine to ensure the efficacy and long-term safety of autologous cell therapies, which is not possible in rodents
<cit> .
in industry settings, non-human primates are important to drug development and are commonly found in drug metabolism and toxicology studies
<cit> .
despite these distinct advantages, drawbacks to non-human primates include greater genetic heterogeneity and higher costs which tend to lead, in turn, to small samples sizes
<cit> .
ultimately these disadvantages contribute to the limited use of non-human primates in biomedical research, particularly in academic settings.
this necessitates the need to optimize study design through careful animal selection, which can only be accomplished by gaining a more thorough understanding of the genetic variation inherent in non-human primates and more specifically the functional effects relative to similar variation in humans.
comparative genetic studies between non-human primates and humans have increased from early candidate gene studies through whole genomes, with limited but significant research now focusing on variation within species.
candidate polymorphism studies in non-human primates, for example, have revealed variation in the dopamine transporter
<cit> , tryptophan hydroxylase  <dig>
<cit> , the serotonin transporter
<cit> , monoamine oxidase a
<cit> , brain-derived neurotrophic factor
<cit> , neuropeptide y
<cit> , and corticotropin-releasing factor
<cit>  that parallel and functionally mimic variation found in humans.
in addition, not only are similar effects seen when these polymorphisms are compared in vitro but similar associations to organismal phenotypes also persist across human and non-human primate species.
g-protein coupled receptors  comprise the largest family of cell surface receptors.
though they share a similar seven transmembrane domain structural homology, they are extraordinarily diverse with the capacity to transduce messages triggered by ligands as varied as photons, organic odorants, nucleotides, nucleosides, peptides, lipids and proteins
<cit> .
consequently, excluding the olfactory subgenome, which represents a distinct class of gpcrs with targeted function
<cit> , this receptor superfamily represents the largest group of druggable targets
<cit>  comprising >50% of pharmacotherapies on the market today.
interestingly, only a third of these gpcrs have been explored for drug development portending a future active area of research for the discovery of novel therapeutics
<cit> .
polymorphisms in gpcrs however can affect drug efficacy through altered ligand binding, receptor activation/inactivation, and/or varied signaling cascades.
characterizing non-human primate variation in gpcrs can therefore complement the study of disease and pharmacotherapies whilst refining the translational capacity of non-human primates in preclinical research.
here the exonic sequence of non-olfactory gpcrs in  <dig> indian-origin rhesus,  <dig> chinese-origin rhesus, and  <dig> cynomolgus macaques was resequenced to gain a better understanding of the natural variation in gpcrs of common non-human primate models.
polymorphisms were then compared to fixed species differences and similar variation in humans.
predicted and known protein structural features were also used to better contextualize the changes and their likely functional effects.
comprehensive polymorphism data in non-human primates not only will facilitate characterization of functional variation at important drug targets and support a better understanding of disease but will also aid in informed a priori selection of animals in preclinical studies and increased translational validity of the non-human primate models ultimately leading to more safe and effective pharmacotherapies and treatments.
RESULTS
over  <dig> million reads were generated representing over  <dig> billion base pairs of sequence from  <dig> animals.
the number of reads per animal ranged from approximately  <dig> million to  <dig> million with a median of just over  <dig>  million.
these reads were aligned to the rhesus genome with the percentage of reads mapped confidently ranging from a minimum of  <dig> % to a maximum  <dig> %, with a median of  <dig> %.
of the  <dig> gpcrs targeted,  <dig> had complete coverage across the gene.
for the remainder, most had localized failures, often a single missing exon or portion of an exon, due to poor or inadequate annotation in the rhesus genome.
it is probable that rna-based approaches or improved annotation would ameliorate many of the failures.
while there were  <dig> animals for which more than 20% of regions were not called, presumably due to suboptimal dna quality or some other manual error in the processing stages, the median coverage for individual animals was  <dig> %.
over  <dig>  snps were identified across all regions and populations .
although the dna capture targeted exons, a large proportion of adjacent introns, upstream, and downstream flanking regions were also resequenced.
within exons, coding regions were the primary focus, though polymorphisms were also found in the 5′ and 3′ untranslated regions  in large numbers.
it is worth noting, however, that 3′ utrs, in particular, may be poorly annotated in the rhesus genome and difficult to comprehensively interrogate.
in coding sequence, nearly  <dig>  coding snps were identified including over  <dig>  non-synonymous and over  <dig>  synonymous snps.
as expected, regions showing the least evolutionary constraint show greater rates of polymorphism and greater numbers of higher frequency polymorphisms.
across non-coding regions, with the notable exception of the 5′ utr, singletons represent roughly 60% of all polymorphisms.
synonymous polymorphisms within coding regions are also at  <dig> %.
in comparison, non-synonymous polymorphisms show a much greater proportion of singletons,  <dig> %, consistent with a slightly deleterious genetic load.
the 5′ utr shows an intermediate proportion of singletons,  <dig> %, perhaps reflective of greater constraint due to a higher density of regulatory elements.
while much fewer, frameshift and nonsense  mutations in coding sequence were also observed.
for the most part these were rare events .
83%  of frameshift mutations were observed in a single individual and nearly 96%  of nonsense mutation were singletons.
among common mutations  private alleles predominated.
one note of caution, however, in that annotation difficulties within the rhesus genome may have overinflated these numbers.
because of the relative likelihood that these mutation will result in functional effects, often creating natural knockouts, particularly common mutations were further examined .
of note, is that five of the thirteen most common of these variants all occur in the celsr <dig> gene, notable for its extensive n-terminal domain.
this and other variation offers fertile ground for potential animal model development going forward.
common alleles defined as those observed in multiple individuals.
population demography
cynomolgus and rhesus macaques, despite being separate species, share polymorphisms
<cit>  and may show some evidence of natural admixture
<cit> .
both cynomolgus macaques and rhesus macaques are widely distributed across southeast asia and cryptic population substructure has been a pervasive problem in biomedical research.
in indian- and chinese-origin rhesus differences in susceptibility and progression of simian immunodeficiency virus  as a model of hiv/aids are the most recognized confounds in research laboratories
<cit>  though other behavioral and physiological differences also certainly exist
<cit> .
using structure
<cit> , rhesus and cynomolgus macaques were readily separated .
it is perhaps noteworthy that those animals that are less unambiguous are those for which fewer reads were generated and had lower levels of coverage across genes.
when only rhesus macaques were considered  the indian and chinese subpopulations readily separated, though three putative indian-origin animals showed significant proportions of chinese admixture, one a 50/ <dig> hybrid and two 75/ <dig> hybrids.
during retrospective investigation these animals were confirmed as known hybrids of the inferred proportions.
indian-origin rhesus macaques were sourced from three locations  but no genetic subdivision was observed.
with regards to the cynomolgus macaques, although all of the individuals used in this study were derived from mauritius stock, unexpected cryptic substructure was observed .
this substructure remains unexplained though recent published studies have indicated similar uncertainty as to the genetic homogeneity of the population
<cit> .
in any case, further study and consideration is warranted.
the demographic history of the subgroups can be confirmed by comparing the allele frequency spectra.
as predicted by population genetics theory, the vast majority of these snps
<cit>  are singletons.
in fact, singletons are overrepresented in all three populations  suggestive of recent population expansion .
again, however, cryptic population substructure in mauritian cynomolgus macaques is supported by an excess of high frequency alleles with a corresponding decline in mid-frequency alleles.
while the two populations of rhesus macaques behave similarly, the allele frequency spectrum of the chinese population appears more similar to that expected under neutrality while the indian population appears to have undergone a more recent population expansion.
these findings are contrary to conventional understandings of the population history of rhesus macaques and to previous genetic studies
<cit> .
it is possible that this discrepancy can be explained through greater artificial selection by humans as the indian rhesus macaques have been bred in biomedical research facilities under strong pressures to avoid inbreeding and to maximize genetic diversity, while chinese populations are more recently derived from wild caught animals.
it is also possible that cryptic differential natural selective regimes otherwise exist between the populations.
as expected, however, a greater percentage of higher frequency non-synonymous snps are lost in all populations, likely representing selection against deleterious alleles.
these findings extend when population specificity of snps is considered .
focusing exclusively on snps found in multiple individuals  the percentage of snps found in both indian and chinese rhesus populations is roughly one third with synonymous snps only slightly more likely to be found in both populations compared to non-synonymous snps .
but while synonymous snps are more likely to be private to indian-origin animals , non-synonymous snps are more often private to chinese-origin rhesus .
if non-synonymous snps are considered to be under greater selective constraint, then these findings are suggestive of either greater constraint in indian-origin animals  or a recent population expansion in these indian animals when compared to the chinese animals.
this latter finding is consistent with the allele frequency spectrum data though shares the same caveats with regard to human selective breeding.
previous studies have demonstrated that cynomolgus macaques share polymorphism with rhesus macaques
<cit> .
using control regions under selective neutrality or presumed constant selective pressures across the species, shared and private polymorphism was used to establish a divergence time of roughly  <dig>  mya and a consistent, if asymmetric, gene flow
<cit> .
studies focusing on the cytochrome p <dig> genes, important modulators of xenobiotic metabolism, have shown a relative increase in private polymorphism thought to perhaps represent the effects of differential selective regimes
<cit> .
interestingly, in gpcrs a greater percentage of non-synonymous snps  are shared between the species than synonymous snps .
this distinction is further muddied, however, when the two rhesus subpopulations are taken into account.
among synonymous snps the majority of shared polymorphisms  are shared among cynomolgus macaques and both rhesus subpopulations, compared to only  <dig> % of non-synonymous snps.
the preponderance of shared synonymous snps is consistent with previous, smaller-scale, findings on non-coding snps
<cit>  and is roughly consistent with expectations under neutrality.
the preponderance and distribution of non-synonymous snps, however, are perhaps indicative of balanced selection.
much of these findings have concentrated on general descriptions of the polymorphism profile of the macaque populations.
while these results have focused on protein-coding regions more likely under negative selective pressures than previous studies of presumably, or more likely, neutral variation, the results have by and large been the same.
to this point, the most notable finding is that non-synonymous polymorphisms seem more likely to be shared between populations than synonymous variation.
while informative, general demographic understandings are better approached through neutral variation and that was not the primary purpose here.
rather, the focus of this study was in identifying and understanding likely functionally relevant variation aimed at improving the usage of macaques as biomedical research models.
the focus on gpcrs, the most common of druggable targets, belies this goal.
distribution of variation
to understand the variation most likely to be functionally relevant in the gpcrs an initial focus was on polymorphism location with regards to secondary structure.
macaque sequences derived from existing annotation coupled with refinements from the consensus resequencing results were aligned with human sequences.
secondary structures for human proteins were pulled from the uniprot database.
the consensus macaque sequences were aligned and fixed divergent sites between macaque and humans were mapped onto secondary sequences.
in accordance with expectations, fixed synonymous mutations were distributed homogenously across the protein without regard for secondary structure.
non-synonymous differences, however, were non-randomly distributed across the secondary structure.
transmembrane domains were significantly more conserved than either intracellular or extracellular domains.
n-terminal and c-terminal domains were the most divergent between taxa and the first and second intracellular domains were the most conserved of the non-transmembrane domains.
these findings are consistent with understandings of gpcr structure and function given that transmembrane domains are expected to be under strong functional constraint to maintain secondary structure and hydrophobicity.
extracellular domains mediate ligand binding with functional residues largely spread across the three loops.
intracellular signaling domains are largely mediated through either the c-terminal domain or the third intracellular loop depending on the nature of the particular gpcr and, therefore, divergence in these domains suggests an evolutionary lability to these functions and drives a need for improved understanding.
as with fixed differences, synonymous snps in each of the populations are distributed evenly and consistently across the protein.
this distribution, driven by neutral mutation rate and largely unaffected by selection, is also seen in the distribution of singletons across the secondary structure .
in comparison, snps that are found in multiple individuals show distribution patterns across the proteins more similar to those seen in divergence with human .
this pattern also holds for human polymorphisms when the cutoff for common snps is arbitrarily placed at 1%.
again it is supposed that rare snps include many slightly deleterious mutations that are destined to be selected out of the population, while more common polymorphisms show patterns consistent with the effects of selective forces.
this can further be explored through the use of functional prediction algorithms.
three unique algorithms were used to classify each of the macaque non-synonymous changes: polyphen-2
<cit> , sift
<cit> , and evod
<cit> .
a consensus of these was used to classify non-synonymous snps as “ambiguous”, “deleterious”, “likely deleterious”, “likely neutral”, or “neutral” after established methods
<cit> .
regardless of the frequency of the snps, singletons or multiples, the percent identified as damaging was statistically the same .
there was also no difference in the proportion of damaging snps within the various populations and subpopulations.
this also did not significantly vary based on the secondary structure domain within the protein or on their distribution between subpopulations .
these findings run contrary to what is seen in humans.
in humans, as one would predict if these predicted deleterious snps are truly damaging, the more common the snp the less likely it is to be classified as deleterious
<cit> .
here not only is there not a correlation between frequency and likelihood of being damaging, but there also seems to be no correlation with secondary structure domain.
this is despite the fact that there does seem to be a correlation between non-synonymous snp frequency and domain as predicted by our conceptual understandings of gpcr structure and function.
there are several possible explanations for this observed phenomenon.
the first and more intriguing is that snps being classified as deleterious are perhaps more likely to change protein function but not necessarily in a selectively negative way.
some portion of these snps could thus be beneficial and driven to higher frequencies.
more likely, however, are much more mundane explanations that these algorithms simply are not designed to work well across species and do not or that the frequencies of alleles observed in these populations are the result of human selective breeding forces in biomedical research colonies and not representative of natural selective effects.
functional variation
regardless, the primary motivation for this study was to understand how functional variation in macaque gpcrs might be used to better understand evolutionary adaptation and the role of macaques as biomedical research models.
one question in particular is how variation in human gpcrs might compare to variation in their macaque orthologs and whether functional effects in humans could be better understood or possibly even modeled in macaques.
to investigate this, human polymorphisms with frequencies greater than one-half of one percent  were drawn from dbsnp.
while arbitrary, these criteria ensured the validity of the snp and at least a modicum of data.
it is important to note, however, that human snps were not chosen by frequencies in specific subpopulations and there are notable issues of ascertainment bias still present in the human data set.
human snps were then mapped to secondary structures following the same methodologies of the macaque polymorphisms and the two data sets were compared.
somewhat unexpectedly, though perhaps not in retrospect, nine recurrent mutations  were identified.
these mutations are present in both humans and macaques.
only snps present in multiple macaque animals were included and the animals sharing these “human” alleles were different so it is reasonably certain that they represent real macaque snps.
these polymorphisms do not represent true trans-species polymorphisms of a shared origin, but rather are recurrent mutations at the same position.
it remains unclear if this is due simply to chance or if there are similar underlying evolutionary pressures.
while there is neither functional information nor phenotypic associations with these snps in humans, it is perhaps interesting to note that consensus predications from polyphen- <dig>  sift, and evod show six of nine as “deleterious” or “likely deleterious”.
in comparison there are only five instances where the same ancestral amino acid was mutated to two different amino acids in human and macaques .
in these cases, the majority of changes are categorized as neutral, though in mrgprx <dig> both human, arg55leu, and macaque, arg55cys, polymorphisms are predicted to be deleterious.
aa: amino acid, maf: minor allele frequency.
aa: amino acid; maf: minor allele frequency.
in total,  <dig> instances were identified in which “common” human variation was found in the same gene and protein secondary structure domain as “common” macaque variation .
these spanned  <dig> distinct genes or roughly one-third of the gpcrs resequenced in this study.
although the majority of these were located in either the n-terminal  or c-terminal  domains, shared variation was found in every secondary structure domain.
the third intercellular domain, often associated with the signaling functions of the gpcrs, had the third greatest amount of shared variation .
further, more than half of all snps identified this way in macaques are predicted to be “deleterious” or “likely deleterious”.
of these, it is useful to highlight some specific examples.
the known parallel functional variation between human and rhesus macaques in oprm <dig> is recapitulated here.
in the n-terminal domain of the mu-opioid receptor, two human polymorphisms c17t  and a118g  show parallel in vitro functional effects with the indian rhesus macaque c77g  mutation
<cit>  as well as parallel phenotypic associations with alcohol consumption and response to naltrexone
<cit> .
this parallel function has already proven to be a useful tool in elucidating the role of the mu-opioid receptor in alcoholism.
prior to the rhesus macaque studies, human work had been inconclusive despite a relatively large number of studies
<cit> .
this variability across studies, inherent in human research due to genetic and environmental heterogeneity, could be quickly and simply teased apart using carefully selected and managed non-human primate models.
in another example, early studies have tentatively linked human variation in adra1a with complex pain and fibromyalgia
<cit>  and specific variation in the third intracellular domain, gly247arg, with receptor pharmacology
<cit> .
while not identical, one common polymorphism is found in the third intracellular loop in macaques, arg266leu, with predicted deleterious effects.
two polymorphisms are also found in the c-terminal domain, lys349arg and arg405his, where associations have also been seen in humans.
several other human variants with putative associations also have possible homologs in macaques.
in the oxytocin receptor , ala218thr has been associated with emotional empathy in humans
<cit> , while ser224cys, in the same receptor domain, is a common polymorphism in chinese-origin rhesus and cynomolgus macaques.
somatostatin receptor  <dig>  variation, phe327ser, has been associated with response to colorectal cancer treatment in humans
<cit> , and rhesus macaques and cynomolgus macaques harbor common polymorphisms ala357asp and met360val, respectively.
variation in follicle-stimulating hormone receptor  and histamine receptor h <dig>  have been associated with polycystic ovarian syndrome
<cit>  and breast cancer
<cit>  respectively and likewise similar polymorphisms may be observed in macaques.
these examples only scratch the surface with the focus here on common human variation, not pathogenic variation.
it is possible that there are additional examples of pathogenic variation that is modeled in macaques or human variation that simply has yet to be recognized as pathogenic due to the vagaries of human research.
common macaque polymorphism may illuminate the functional relevance of human variation even in the absence of known human associations.
variation found in the same genes and secondary structures in humans and macaques offer potentially informative targets for studies of functionally similar, though evolutionarily distinct, variation across species and for the improvement of understanding the molecular underpinnings of disease.
CONCLUSIONS
drug discovery and translational medicine benefit from strong animal models.
for too long poor animal models have led researchers down the wrong paths, leading, perhaps, to novel understandings and interesting results, but not to improved treatments in humans that have been promised.
in part, the scientific community has been playing the cards it was dealt, too quick to believe that shared phenotypes implied a shared molecular basis.
now, however, the revolution in sequencing technologies allows us to look closer at the molecular basis of disease than has ever been possible and, in doing so, we can more easily identify when shared phenotypes do share molecular bases and when they do not.
moreover, we can identify where similar molecular and genetic foundations exist, but do not lead to the same phenotypic effects.
non-human primates have long been known to share genetic and physiological similarities with humans.
this has made them the gold standard for preclinical research, though one for which it has not always been clear if the benefits outweighed the price.
by better understanding the genetics of non-human primates we lay clear the benefits, demonstrating where genetic similarities exist with humans and where non-human primates are most likely to be beneficial.
we also develop tools for maximizing the utility of non-human primates, ensuring that when they are used as biomedical research models they are used appropriately and result in the greatest power.
here we catalog the polymorphism in the gpcrs of rhesus macaques of indian and chinese origin and mauritian cynomolgus macaques.
together these species represent the most commonly used non-human primate biomedical research models and the genes represent the single largest family of drug targets.
this information can be used going forward to develop improved animal models and to better understand gene-phenotype associations.
by improving our animal models we improve the ability of our science to be translational and ultimately to bring basic research to bear on issues of human health.

 animal research has provided the scientific community with extraordinary advances in medicine from the development of vaccines to the prevention and treatment of diseases. unfortunately at present 85% of novel therapeutics fail in preclinical and early phase clinical trials and of the therapies that reach late phase trials an additional 50% fall short due to an inability to demonstrate efficacy and safety. reasons for these shortcomings include low patient recruitment, poor study design, and ineffective use of animal models. coupled with soaring drug development costs including both financial commitments and in years of labor, these shortfalls necessitate a biological and economic need for fundamental changes in the bench to bedside process. furthermore, with advances in genome sequencing technologies there is a growing awareness that animal models fall short in terms of predictive power. a recent study comparing the genomic responses of human inflammatory diseases to mouse models, for example, suggested that mice poorly mimic the human genetic response. continued progress in the understanding of human disease pathologies and the development of safe and effective therapies demands a more comprehensive understanding of animals in preclinical research. although greater numbers of rodents are used in biomedical research, non-human primates are the gold standard of animal models in preclinical research offering advantages which include greater similarities in genome organization and sequence, behavior, and physiology. the rhesus and cynomolgus macaque are two of the most commonly used non-human primate species in research laboratories, sharing ~ % of their genome with humans. in academic research non-human primate use is most common in the fields of microbiology, biochemistry/pharmacology, and neuroscience. because of similarities in physiology and the central nervous system, non-human primates, for example, are crucial in stem cell-based regenerative medicine to ensure the efficacy and long-term safety of autologous cell therapies, which is not possible in rodents. in industry settings, non-human primates are important to drug development and are commonly found in drug metabolism and toxicology studies. despite these distinct advantages, drawbacks to non-human primates include greater genetic heterogeneity and higher costs which tend to lead, in turn, to small samples sizes. ultimately these disadvantages contribute to the limited use of non-human primates in biomedical research, particularly in academic settings. this necessitates the need to optimize study design through careful animal selection, which can only be accomplished by gaining a more thorough understanding of the genetic variation inherent in non-human primates and more specifically the functional effects relative to similar variation in humans. comparative genetic studies between non-human primates and humans have increased from early candidate gene studies through whole genomes, with limited but significant research now focusing on variation within species. candidate polymorphism studies in non-human primates, for example, have revealed variation in the dopamine transporter, tryptophan hydroxylase, the serotonin transporter, monoamine oxidase a, brain-derived neurotrophic factor, neuropeptide y, and corticotropin-releasing factor that parallel and functionally mimic variation found in humans. in addition, not only are similar effects seen when these polymorphisms are compared in vitro but similar associations to organismal phenotypes also persist across human and non-human primate species. g-protein coupled receptors comprise the largest family of cell surface receptors. though they share a similar seven transmembrane domain structural homology, they are extraordinarily diverse with the capacity to transduce messages triggered by ligands as varied as photons, organic odorants, nucleotides, nucleosides, peptides, lipids and proteins. consequently, excluding the olfactory subgenome, which represents a distinct class of gpcrs with targeted function, this receptor superfamily represents the largest group of druggable targets comprising >50% of pharmacotherapies on the market today. interestingly, only a third of these gpcrs have been explored for drug development portending a future active area of research for the discovery of novel therapeutics. polymorphisms in gpcrs however can affect drug efficacy through altered ligand binding, receptor activation/inactivation, and/or varied signaling cascades. characterizing non-human primate variation in gpcrs can therefore complement the study of disease and pharmacotherapies whilst refining the translational capacity of non-human primates in preclinical research. here the exonic sequence of non-olfactory gpcrs in indian-origin rhesus, chinese-origin rhesus, and cynomolgus macaques was resequenced to gain a better understanding of the natural variation in gpcrs of common non-human primate models. polymorphisms were then compared to fixed species differences and similar variation in humans. predicted and known protein structural features were also used to better contextualize the changes and their likely functional effects. comprehensive polymorphism data in non-human primates not only will facilitate characterization of functional variation at important drug targets and support a better understanding of disease but will also aid in informed a priori selection of animals in preclinical studies and increased translational validity of the non-human primate models ultimately leading to more safe and effective pharmacotherapies and treatments.  over million reads were generated representing over billion base pairs of sequence from animals. the number of reads per animal ranged from approximately million to million with a median of just over million. these reads were aligned to the rhesus genome with the percentage of reads mapped confidently ranging from a minimum of % to a maximum %, with a median of %. of the gpcrs targeted, had complete coverage across the gene. for the remainder, most had localized failures, often a single missing exon or portion of an exon, due to poor or inadequate annotation in the rhesus genome. it is probable that rna-based approaches or improved annotation would ameliorate many of the failures. while there were animals for which more than 20% of regions were not called, presumably due to suboptimal dna quality or some other manual error in the processing stages, the median coverage for individual animals was %. over snps were identified across all regions and populations. although the dna capture targeted exons, a large proportion of adjacent introns, upstream, and downstream flanking regions were also resequenced. within exons, coding regions were the primary focus, though polymorphisms were also found in the 5′ and 3′ untranslated regions in large numbers. it is worth noting, however, that 3′ utrs, in particular, may be poorly annotated in the rhesus genome and difficult to comprehensively interrogate. in coding sequence, nearly coding snps were identified including over non-synonymous and over synonymous snps. as expected, regions showing the least evolutionary constraint show greater rates of polymorphism and greater numbers of higher frequency polymorphisms. across non-coding regions, with the notable exception of the 5′ utr, singletons represent roughly 60% of all polymorphisms. synonymous polymorphisms within coding regions are also at %. in comparison, non-synonymous polymorphisms show a much greater proportion of singletons, %, consistent with a slightly deleterious genetic load. the 5′ utr shows an intermediate proportion of singletons, %, perhaps reflective of greater constraint due to a higher density of regulatory elements. while much fewer, frameshift and nonsense mutations in coding sequence were also observed. for the most part these were rare events. 83% of frameshift mutations were observed in a single individual and nearly 96% of nonsense mutation were singletons. among common mutations private alleles predominated. one note of caution, however, in that annotation difficulties within the rhesus genome may have overinflated these numbers. because of the relative likelihood that these mutation will result in functional effects, often creating natural knockouts, particularly common mutations were further examined. of note, is that five of the thirteen most common of these variants all occur in the celsr gene, notable for its extensive n-terminal domain. this and other variation offers fertile ground for potential animal model development going forward. common alleles defined as those observed in multiple individuals. population demography cynomolgus and rhesus macaques, despite being separate species, share polymorphisms and may show some evidence of natural admixture. both cynomolgus macaques and rhesus macaques are widely distributed across southeast asia and cryptic population substructure has been a pervasive problem in biomedical research. in indian- and chinese-origin rhesus differences in susceptibility and progression of simian immunodeficiency virus as a model of hiv/aids are the most recognized confounds in research laboratories though other behavioral and physiological differences also certainly exist. using structure, rhesus and cynomolgus macaques were readily separated. it is perhaps noteworthy that those animals that are less unambiguous are those for which fewer reads were generated and had lower levels of coverage across genes. when only rhesus macaques were considered the indian and chinese subpopulations readily separated, though three putative indian-origin animals showed significant proportions of chinese admixture, one a 50/ hybrid and two 75/ hybrids. during retrospective investigation these animals were confirmed as known hybrids of the inferred proportions. indian-origin rhesus macaques were sourced from three locations but no genetic subdivision was observed. with regards to the cynomolgus macaques, although all of the individuals used in this study were derived from mauritius stock, unexpected cryptic substructure was observed. this substructure remains unexplained though recent published studies have indicated similar uncertainty as to the genetic homogeneity of the population. in any case, further study and consideration is warranted. the demographic history of the subgroups can be confirmed by comparing the allele frequency spectra. as predicted by population genetics theory, the vast majority of these snps are singletons. in fact, singletons are overrepresented in all three populations suggestive of recent population expansion. again, however, cryptic population substructure in mauritian cynomolgus macaques is supported by an excess of high frequency alleles with a corresponding decline in mid-frequency alleles. while the two populations of rhesus macaques behave similarly, the allele frequency spectrum of the chinese population appears more similar to that expected under neutrality while the indian population appears to have undergone a more recent population expansion. these findings are contrary to conventional understandings of the population history of rhesus macaques and to previous genetic studies. it is possible that this discrepancy can be explained through greater artificial selection by humans as the indian rhesus macaques have been bred in biomedical research facilities under strong pressures to avoid inbreeding and to maximize genetic diversity, while chinese populations are more recently derived from wild caught animals. it is also possible that cryptic differential natural selective regimes otherwise exist between the populations. as expected, however, a greater percentage of higher frequency non-synonymous snps are lost in all populations, likely representing selection against deleterious alleles. these findings extend when population specificity of snps is considered. focusing exclusively on snps found in multiple individuals the percentage of snps found in both indian and chinese rhesus populations is roughly one third with synonymous snps only slightly more likely to be found in both populations compared to non-synonymous snps. but while synonymous snps are more likely to be private to indian-origin animals, non-synonymous snps are more often private to chinese-origin rhesus. if non-synonymous snps are considered to be under greater selective constraint, then these findings are suggestive of either greater constraint in indian-origin animals or a recent population expansion in these indian animals when compared to the chinese animals. this latter finding is consistent with the allele frequency spectrum data though shares the same caveats with regard to human selective breeding. previous studies have demonstrated that cynomolgus macaques share polymorphism with rhesus macaques. using control regions under selective neutrality or presumed constant selective pressures across the species, shared and private polymorphism was used to establish a divergence time of roughly mya and a consistent, if asymmetric, gene flow. studies focusing on the cytochrome p genes, important modulators of xenobiotic metabolism, have shown a relative increase in private polymorphism thought to perhaps represent the effects of differential selective regimes. interestingly, in gpcrs a greater percentage of non-synonymous snps are shared between the species than synonymous snps. this distinction is further muddied, however, when the two rhesus subpopulations are taken into account. among synonymous snps the majority of shared polymorphisms are shared among cynomolgus macaques and both rhesus subpopulations, compared to only % of non-synonymous snps. the preponderance of shared synonymous snps is consistent with previous, smaller-scale, findings on non-coding snps and is roughly consistent with expectations under neutrality. the preponderance and distribution of non-synonymous snps, however, are perhaps indicative of balanced selection. much of these findings have concentrated on general descriptions of the polymorphism profile of the macaque populations. while these results have focused on protein-coding regions more likely under negative selective pressures than previous studies of presumably, or more likely, neutral variation, the results have by and large been the same. to this point, the most notable finding is that non-synonymous polymorphisms seem more likely to be shared between populations than synonymous variation. while informative, general demographic understandings are better approached through neutral variation and that was not the primary purpose here. rather, the focus of this study was in identifying and understanding likely functionally relevant variation aimed at improving the usage of macaques as biomedical research models. the focus on gpcrs, the most common of druggable targets, belies this goal. distribution of variation to understand the variation most likely to be functionally relevant in the gpcrs an initial focus was on polymorphism location with regards to secondary structure. macaque sequences derived from existing annotation coupled with refinements from the consensus resequencing results were aligned with human sequences. secondary structures for human proteins were pulled from the uniprot database. the consensus macaque sequences were aligned and fixed divergent sites between macaque and humans were mapped onto secondary sequences. in accordance with expectations, fixed synonymous mutations were distributed homogenously across the protein without regard for secondary structure. non-synonymous differences, however, were non-randomly distributed across the secondary structure. transmembrane domains were significantly more conserved than either intracellular or extracellular domains. n-terminal and c-terminal domains were the most divergent between taxa and the first and second intracellular domains were the most conserved of the non-transmembrane domains. these findings are consistent with understandings of gpcr structure and function given that transmembrane domains are expected to be under strong functional constraint to maintain secondary structure and hydrophobicity. extracellular domains mediate ligand binding with functional residues largely spread across the three loops. intracellular signaling domains are largely mediated through either the c-terminal domain or the third intracellular loop depending on the nature of the particular gpcr and, therefore, divergence in these domains suggests an evolutionary lability to these functions and drives a need for improved understanding. as with fixed differences, synonymous snps in each of the populations are distributed evenly and consistently across the protein. this distribution, driven by neutral mutation rate and largely unaffected by selection, is also seen in the distribution of singletons across the secondary structure. in comparison, snps that are found in multiple individuals show distribution patterns across the proteins more similar to those seen in divergence with human. this pattern also holds for human polymorphisms when the cutoff for common snps is arbitrarily placed at 1%. again it is supposed that rare snps include many slightly deleterious mutations that are destined to be selected out of the population, while more common polymorphisms show patterns consistent with the effects of selective forces. this can further be explored through the use of functional prediction algorithms. three unique algorithms were used to classify each of the macaque non-synonymous changes: polyphen-2, sift, and evod. a consensus of these was used to classify non-synonymous snps as “ambiguous”, “deleterious”, “likely deleterious”, “likely neutral”, or “neutral” after established methods. regardless of the frequency of the snps, singletons or multiples, the percent identified as damaging was statistically the same. there was also no difference in the proportion of damaging snps within the various populations and subpopulations. this also did not significantly vary based on the secondary structure domain within the protein or on their distribution between subpopulations. these findings run contrary to what is seen in humans. in humans, as one would predict if these predicted deleterious snps are truly damaging, the more common the snp the less likely it is to be classified as deleterious. here not only is there not a correlation between frequency and likelihood of being damaging, but there also seems to be no correlation with secondary structure domain. this is despite the fact that there does seem to be a correlation between non-synonymous snp frequency and domain as predicted by our conceptual understandings of gpcr structure and function. there are several possible explanations for this observed phenomenon. the first and more intriguing is that snps being classified as deleterious are perhaps more likely to change protein function but not necessarily in a selectively negative way. some portion of these snps could thus be beneficial and driven to higher frequencies. more likely, however, are much more mundane explanations that these algorithms simply are not designed to work well across species and do not or that the frequencies of alleles observed in these populations are the result of human selective breeding forces in biomedical research colonies and not representative of natural selective effects. functional variation regardless, the primary motivation for this study was to understand how functional variation in macaque gpcrs might be used to better understand evolutionary adaptation and the role of macaques as biomedical research models. one question in particular is how variation in human gpcrs might compare to variation in their macaque orthologs and whether functional effects in humans could be better understood or possibly even modeled in macaques. to investigate this, human polymorphisms with frequencies greater than one-half of one percent were drawn from dbsnp. while arbitrary, these criteria ensured the validity of the snp and at least a modicum of data. it is important to note, however, that human snps were not chosen by frequencies in specific subpopulations and there are notable issues of ascertainment bias still present in the human data set. human snps were then mapped to secondary structures following the same methodologies of the macaque polymorphisms and the two data sets were compared. somewhat unexpectedly, though perhaps not in retrospect, nine recurrent mutations were identified. these mutations are present in both humans and macaques. only snps present in multiple macaque animals were included and the animals sharing these “human” alleles were different so it is reasonably certain that they represent real macaque snps. these polymorphisms do not represent true trans-species polymorphisms of a shared origin, but rather are recurrent mutations at the same position. it remains unclear if this is due simply to chance or if there are similar underlying evolutionary pressures. while there is neither functional information nor phenotypic associations with these snps in humans, it is perhaps interesting to note that consensus predications from polyphen- sift, and evod show six of nine as “deleterious” or “likely deleterious”. in comparison there are only five instances where the same ancestral amino acid was mutated to two different amino acids in human and macaques. in these cases, the majority of changes are categorized as neutral, though in mrgprx both human, arg55leu, and macaque, arg55cys, polymorphisms are predicted to be deleterious. aa: amino acid, maf: minor allele frequency. aa: amino acid; maf: minor allele frequency. in total, instances were identified in which “common” human variation was found in the same gene and protein secondary structure domain as “common” macaque variation. these spanned distinct genes or roughly one-third of the gpcrs resequenced in this study. although the majority of these were located in either the n-terminal or c-terminal domains, shared variation was found in every secondary structure domain. the third intercellular domain, often associated with the signaling functions of the gpcrs, had the third greatest amount of shared variation. further, more than half of all snps identified this way in macaques are predicted to be “deleterious” or “likely deleterious”. of these, it is useful to highlight some specific examples. the known parallel functional variation between human and rhesus macaques in oprm is recapitulated here. in the n-terminal domain of the mu-opioid receptor, two human polymorphisms c17t and a118g show parallel in vitro functional effects with the indian rhesus macaque c77g mutation as well as parallel phenotypic associations with alcohol consumption and response to naltrexone. this parallel function has already proven to be a useful tool in elucidating the role of the mu-opioid receptor in alcoholism. prior to the rhesus macaque studies, human work had been inconclusive despite a relatively large number of studies. this variability across studies, inherent in human research due to genetic and environmental heterogeneity, could be quickly and simply teased apart using carefully selected and managed non-human primate models. in another example, early studies have tentatively linked human variation in adra1a with complex pain and fibromyalgia and specific variation in the third intracellular domain, gly247arg, with receptor pharmacology. while not identical, one common polymorphism is found in the third intracellular loop in macaques, arg266leu, with predicted deleterious effects. two polymorphisms are also found in the c-terminal domain, lys349arg and arg405his, where associations have also been seen in humans. several other human variants with putative associations also have possible homologs in macaques. in the oxytocin receptor, ala218thr has been associated with emotional empathy in humans, while ser224cys, in the same receptor domain, is a common polymorphism in chinese-origin rhesus and cynomolgus macaques. somatostatin receptor variation, phe327ser, has been associated with response to colorectal cancer treatment in humans, and rhesus macaques and cynomolgus macaques harbor common polymorphisms ala357asp and met360val, respectively. variation in follicle-stimulating hormone receptor and histamine receptor h have been associated with polycystic ovarian syndrome and breast cancer respectively and likewise similar polymorphisms may be observed in macaques. these examples only scratch the surface with the focus here on common human variation, not pathogenic variation. it is possible that there are additional examples of pathogenic variation that is modeled in macaques or human variation that simply has yet to be recognized as pathogenic due to the vagaries of human research. common macaque polymorphism may illuminate the functional relevance of human variation even in the absence of known human associations. variation found in the same genes and secondary structures in humans and macaques offer potentially informative targets for studies of functionally similar, though evolutionarily distinct, variation across species and for the improvement of understanding the molecular underpinnings of disease.  drug discovery and translational medicine benefit from strong animal models. for too long poor animal models have led researchers down the wrong paths, leading, perhaps, to novel understandings and interesting results, but not to improved treatments in humans that have been promised. in part, the scientific community has been playing the cards it was dealt, too quick to believe that shared phenotypes implied a shared molecular basis. now, however, the revolution in sequencing technologies allows us to look closer at the molecular basis of disease than has ever been possible and, in doing so, we can more easily identify when shared phenotypes do share molecular bases and when they do not. moreover, we can identify where similar molecular and genetic foundations exist, but do not lead to the same phenotypic effects. non-human primates have long been known to share genetic and physiological similarities with humans. this has made them the gold standard for preclinical research, though one for which it has not always been clear if the benefits outweighed the price. by better understanding the genetics of non-human primates we lay clear the benefits, demonstrating where genetic similarities exist with humans and where non-human primates are most likely to be beneficial. we also develop tools for maximizing the utility of non-human primates, ensuring that when they are used as biomedical research models they are used appropriately and result in the greatest power. here we catalog the polymorphism in the gpcrs of rhesus macaques of indian and chinese origin and mauritian cynomolgus macaques. together these species represent the most commonly used non-human primate biomedical research models and the genes represent the single largest family of drug targets. this information can be used going forward to develop improved animal models and to better understand gene-phenotype associations. by improving our animal models we improve the ability of our science to be translational and ultimately to bring basic research to bear on issues of human health. 
while the vast majority of these snps are singletons, roughly non-synonymous and synonymous snps were found in multiple individuals. gpcrs also represent the single largest class of drug targets and associated pharmacogenetic effects are modulated, in part, by polymorphisms. recently, non-human primate models have been developed focusing on naturally-occurring, functionally-parallel polymorphisms in candidate genes. snp frequencies in macaques follow a similar pattern of divergence from humans and new polymorphisms in primates have been identified that may parallel those seen in humans, helping to establish better non-human primate models of disease. as expected, regions showing the least evolutionary constraint show greater rates of polymorphism and greater numbers of higher frequency polymorphisms. in all three populations, polymorphism and divergence is highly concentrated in n-terminal and c-terminal domains and the third intracellular loop region of gpcrs, regions critical to ligand-binding and signaling. g-protein coupled receptors play an inordinately large role in human health.
 g-protein coupled receptors play an inordinately large role in human health. variation in the genes that encode these receptors is associated with numerous disorders across the entire spectrum of disease. gpcrs also represent the single largest class of drug targets and associated pharmacogenetic effects are modulated, in part, by polymorphisms. recently, non-human primate models have been developed focusing on naturally-occurring, functionally-parallel polymorphisms in candidate genes. this work aims to extend those studies broadly across the roughly non-olfactory gpcrs. initial efforts include resequencing indian-origin rhesus macaques, chinese-origin rhesus macaques, and cynomolgus macaques.  using the agilent target enrichment system, capture baits were designed for gpcrs off the human and rhesus exonic sequence. using next generation sequencing technologies, nearly snps were identified in coding sequences including over non-synonymous and more than synonymous protein-coding snps. as expected, regions showing the least evolutionary constraint show greater rates of polymorphism and greater numbers of higher frequency polymorphisms. while the vast majority of these snps are singletons, roughly non-synonymous and synonymous snps were found in multiple individuals.  in all three populations, polymorphism and divergence is highly concentrated in n-terminal and c-terminal domains and the third intracellular loop region of gpcrs, regions critical to ligand-binding and signaling. snp frequencies in macaques follow a similar pattern of divergence from humans and new polymorphisms in primates have been identified that may parallel those seen in humans, helping to establish better non-human primate models of disease. resequencingsingle-nucleotide polymorphismindian-origin rhesus macaqueschinese-origin rhesus macaquescynomolgus macaques 