Primate Gene Research Articles

Retained duplications and deletions of CYP2C genes among primates

The human genome encodes about 60 functional enzymes of the cytochrome P450 superfamily, including four functional enzymes of the cytochrome P450 2C (CYP2C) subfamily. These enzymes have been shown to metabolize drugs and xenobiotic toxins, such as those in the diet, and are therefore of great importance for biomedical research and applications. While the pharmacology of P450 enzymes has been studied extensively, our understanding of molecular evolution of this gene family is incomplete, in part because a great variation exists in the number of CYP2C homologs across genomes. In humans, the enzymes encoded by these genes are responsible for the metabolism of more than 20% of clinical drugs, but this is not the naturalistic function of these enzymes, which is the metabolism of xenobiotics such as plant secondary metabolites. In this paper, we sought to correlate evolutionary relationships among primate CYP2C genes with known dietary profiles from these species, testing the hypothesis that these genes have evolved under the pressure of dietary toxins. Aside from a small number of deeply divergent genes, primate CYP2C paralogs form three separate clades: CYP2C18, CYP2C9/CYP2C19, and CYP2C8/CYP2C20. Our results showed that the CYP2C18 gene has been separately lost in Nomascus leucogenys and the Panini genomes, and there is no evidence that this gene has been under any positive selection among primates. While CYP2C20 has been retained in cercopithecoids, orthologous loci were separately lost in platyrrhines and hominoids. Notably, nine codons exhibited signature of positive selection. Finally, the CYP2C19 locus was duplicated in basal catarrhines, resulting in the birth of CYP2C9; but the ancestral locus was only retained in hominoid taxa. Overall, our findings support the hypothesis that primate CYP2C genes have evolved in response to selective pressures provided by dietary toxins, although not all gene clusters have evolved in the same manner. Our results may indicate an evolutionarily deep difference in ecology or physiology among higher-order primate taxa.

Similarity and variation in the insulin-like growth factor 2 - H19 locus in primates.

Insulin-like growth factor 2 (IGF2), a small, secreted protein, is critical for fetal and prenatal growth in humans and other mammals. The IGF2 gene and its mouse homolog comprise part of a conserved linkage group that is regulated by parental imprinting, with IGF2/ Igf2 being expressed from the paternal chromosome, and the adjacent H19 gene from the maternal chromosome. By using information extracted from public genomic and gene expression databases, I have now analyzed this locus in nine nonhuman primate species representing over 60 million years of evolutionary divergence from a common progenitor. Both IGF2 and H19 genes and the entire locus have been conserved among these primates. Each primate IGF2 gene except for gibbon and marmoset is composed of 10 exons and contains five potential promoters, each with distinctive 5'-untranslated exons. Similarly, except for marmoset and mouse lemur, H19 consists of six exons and has two promoters. DNA sequence conservation is high, not only in orthologous exons and promoters, but also in a putative imprinting control region located 5' to H19 and in multiple potential distal enhancer elements found 3' to H19. Collectively, these results support the hypothesis that common regulatory processes shaped the IGF2 - H19 locus before the onset of primate speciation more than 85 million years ago. This study also leads to the conclusion that inaccuracies in data presentation in genetic repositories could limit our ability to develop novel insights about roles of individual genes and multigene loci in mammalian physiology and disease.

The non-human primate kidney transcriptome in fetal development.

Little is known about the repertoire of non-human primate kidney genes expressed throughout development. The present work establishes an understanding of the primate renal transcriptome during fetal development in the context of renal maturation. The baboon kidney transcriptome was characterized at 60-day gestation (DG), 90 DG, 125 DG, 140 DG, 160 DG and adulthood (6-12years) using gene arrays and validated by QRT-PCR. Pathway and cluster analyses were used to characterize gene expression in the context of biological pathways. Pathway analysis indicated activation of pathways not previously reported as relevant to kidney development. Cluster analysis also revealed gene splice variants with discordant expression profiles during development. This study provides the first detailed genetic analysis of the developing primate kidney, and our findings of discordant expression of gene splice variants suggest that gene arrays likely provide a simplified view and demonstrate the need to study the fetal renal proteome.

Hormonal vitamin D up-regulates tissue-specific PD-L1 and PD-L2 surface glycoprotein expression in humans but not mice

PD-L1 (programmed death ligand 1) and PD-L2 are cell-surface glycoproteins that interact with programmed death 1 (PD-1) on T cells to attenuate inflammation. PD-1 signaling has attracted intense interest for its role in a pathophysiological context: suppression of anti-tumor immunity. Similarly, vitamin D signaling has been increasingly investigated for its non-classical actions in stimulation of innate immunity and suppression of inflammatory responses. Here, we show that hormonal 1,25-dihydroxyvitamin D (1,25D) is a direct transcriptional inducer of the human genes encoding PD-L1 and PD-L2 through the vitamin D receptor, a ligand-regulated transcription factor. 1,25D stimulated transcription of the gene encoding PD-L1 in epithelial and myeloid cells, whereas the gene encoding the more tissue-restricted PD-L2 was regulated only in myeloid cells. We identified and characterized vitamin D response elements (VDREs) located in both genes and showed that 1,25D treatment induces cell-surface expression of PD-L1 in epithelial and myeloid cells. In co-culture experiments with primary human T cells, epithelial cells pretreated with 1,25D suppressed activation of CD4+ and CD8+ cells and inhibited inflammatory cytokine production in a manner that was abrogated by anti-PD-L1 blocking antibody. Consistent with previous observations of species-specific regulation of immunity by vitamin D, the VDREs are present in primate genes, but neither the VDREs nor the regulation by 1,25D is present in mice. These findings reinforce the physiological role of 1,25D in controlling inflammatory immune responses but may represent a double-edged sword, as they suggest that elevated vitamin D signaling in humans could suppress anti-tumor immunity.

NLRP11 attenuates Toll-like receptor signalling by targeting TRAF6 for degradation via the ubiquitin ligase RNF19A

The adaptor protein TRAF6 has a central function in Toll-like receptor (TLR) signalling, yet the molecular mechanisms controlling its activity and stability are unclear. Here we show that NLRP11, a primate specific gene, inhibits TLR signalling by targeting TRAF6 for degradation. NLRP11 recruits the ubiquitin ligase RNF19A to catalyze K48-linked ubiquitination of TRAF6 at multiple sites, thereby leading to the degradation of TRAF6. Furthermore, deficiency in either NLRP11 or RNF19A abrogates K48-linked ubiquitination and degradation of TRAF6, which promotes activation of NF-κB and MAPK signalling and increases the production of proinflammatory cytokines. Therefore, our findings identify NLRP11 as a conserved negative regulator of TLR signalling in primate cells and reveal a mechanism by which the NLRP11-RNF19A axis targets TRAF6 for degradation.

Comparison of Transcriptomic Platforms for Analysis of Whole Blood from Ebola-Infected Cynomolgus Macaques

Ebola virus disease (EVD) is a serious illness with mortality rates of 20–90% in various outbreaks. EVD is characterized by robust virus replication and strong host inflammatory response. Analyzing host immune responses has increasingly involved multimodal approaches including transcriptomics to profile gene expression. We studied cynomolgus macaques exposed to Ebola virus Makona via different routes with the intent of comparing RNA-Seq to a NanoString nCounter codeset targeting 769 non-human primate (NHP) genes. RNA-Seq analysis of serial blood samples showed different routes led to the same overall transcriptional response seen in previously reported EBOV-exposed NHP studies. Both platforms displayed a strong correlation in gene expression patterns, including a strong induction of innate immune response genes at early times post-exposure, and neutrophil-associated genes at later time points. A 41-gene classifier was tested in both platforms for ability to cluster samples by infection status. Both NanoString and RNA-Seq could be used to predict relative abundances of circulating immune cell populations that matched traditional hematology. This demonstrates the complementarity of RNA-Seq and NanoString. Moreover, the development of an NHP-specific NanoString codeset should augment studies of filoviruses and other high containment infectious diseases without the infrastructure requirements of RNA-Seq technology.

Faster Evolving Primate Genes Are More Likely to Duplicate.

An attractive and long-standing hypothesis regarding the evolution of genes after duplication posits that the duplication event creates new evolutionary possibilities by releasing a copy of the gene from constraint. Apparent support was found in numerous analyses, particularly, the observation of higher rates of evolution in duplicated as compared with singleton genes. Could it, instead, be that more duplicable genes (owing to mutation, fixation, or retention biases) are intrinsically faster evolving? To uncouple the measurement of rates of evolution from the determination of duplicate or singleton status, we measure the rates of evolution in singleton genes in outgroup primate lineages but classify these genes as to whether they have duplicated or not in a crown group of great apes. We find that rates of evolution are higher in duplicable genes prior to the duplication event. In part this is owing to a negative correlation between coding sequence length and rate of evolution, coupled with a bias toward smaller genes being more duplicable. The effect is masked by difference in expression rate between duplicable genes and singletons. Additionally, in contradiction to the classical assumption, we find no convincing evidence for an increase in after duplication, nor for rate asymmetry between duplicates. We conclude that high rates of evolution of duplicated genes are not solely a consequence of the duplication event, but are rather a predictor of duplicability. These results are consistent with a model in which successful gene duplication events in mammals are skewed toward events of minimal phenotypic impact.

Transposable elements are the primary source of novelty in primate gene regulation.

Gene regulation shapes the evolution of phenotypic diversity. We investigated the evolution of liver promoters and enhancers in six primate species using ChIP-seq (H3K27ac and H3K4me1) to profile cis-regulatory elements (CREs) and using RNA-seq to characterize gene expression in the same individuals. To quantify regulatory divergence, we compared CRE activity across species by testing differential ChIP-seq read depths directly measured for orthologous sequences. We show that the primate regulatory landscape is largely conserved across the lineage, with 63% of the tested human liver CREs showing similar activity across species. Conserved CRE function is associated with sequence conservation, proximity to coding genes, cell-type specificity, and transcription factor binding. Newly evolved CREs are enriched in immune response and neurodevelopmental functions. We further demonstrate that conserved CREs bind master regulators, suggesting that while CREs contribute to species adaptation to the environment, core functions remain intact. Newly evolved CREs are enriched in young transposable elements (TEs), including Long-Terminal-Repeats (LTRs) and SINE-VNTR-Alus (SVAs), that significantly affect gene expression. Conversely, only 16% of conserved CREs overlap TEs. We tested the cis-regulatory activity of 69 TE subfamilies by luciferase reporter assays, spanning all major TE classes, and showed that 95.6% of tested TEs can function as either transcriptional activators or repressors. In conclusion, we demonstrated the critical role of TEs in primate gene regulation and illustrated potential mechanisms underlying evolutionary divergence among the primate species through the noncoding genome.

Identification and characterization of the tyrosinase gene (TYR) and its transcript variants (TYR_1 and TYR_2) in the crab-eating macaque (Macaca fascicularis)

Tyrosinase is a copper-containing enzyme that regulates melanin biosynthesis and is encoded by the tyrosinase (TYR) gene. Previous studies demonstrated that mutations in TYR could lead to oculocutaneous albinism type 1 (OCA1) owing to the failure of melanin formation. Although a previous study found that albinism in the rhesus monkey was derived from a mutation in TYR, the identification and characterization of this gene in non-human primates has not been achieved thus far. Thus, using the rapid amplification of cDNA ends (RACE) and internal reverse transcription PCR (RT-PCR) we identified the full-length sequence of TYR in the crab-eating macaque, and two different transcript variants (TYR_1 and TYR_2). While TYR_1 comprised five exons and its coding sequence was highly similar to that of humans, TYR_2 comprised four exons and was generated by a third-exon-skipping event. Interestingly, these two transcripts were also present in the African green monkey (Old World monkey) and the common marmoset (New World monkey). Deduced amino acid sequence analyses revealed that TYR_2 had a shorter C-terminal region than TYR_1 owing to the exon-skipping event. Thus, the present study is the first to identify and characterize a full-length TYR gene in a non-human primate, while the further validation of the third-exon-skipping in TYR indicates that this event is well conserved in the primate lineage. Therefore, this study provides useful and important information for the study of albinism using non-human primate models.

Alzheimer’s disease related genes during primate evolution

During primate evolution, the neuronal and cognition-related genes have evolved rapidly. These genes seem to induce neurological illnesses such as Alzheimer’s disease (AD). In this study, we analyzed genes APOE, TOMM40, and PICALM known as the risk factors of AD. We performed bioinformatics analyses in relation to evolution, phylogeny, and protein structure for those genes in humans, Neanderthals, chimpanzees, bonobos, gorillas, orangutans, crab-eating monkeys, and rhesus monkeys. Cholesterol-related genes showed relatively rapid evolution toward a lower risk of AD. Neanderthals showed relatively higher polymorphism in genes APOE, TOMM40, and PICALM than humans did. Phylogeny indicated different topologies in the trichotomy of humans, chimpanzees, and gorillas in terms of genes APOE, TOMM40, and PICALM. These results provide to hominin-specific patterns in three genes, and give clues to the modern human-specific traits of AD and shed light on further functional research helping to understand AD.

The Binding Sites of miR-619-5p in the mRNAs of Human and Orthologous Genes

BackgroundNormally, one miRNA interacts with the mRNA of one gene. However, there are miRNAs that can bind to many mRNAs, and one mRNA can be the target of many miRNAs. This significantly complicates the study of the properties of miRNAs and their diagnostic and medical applications.ResultsThe search of 2,750 human microRNAs (miRNAs) binding sites in 12,175 mRNAs of human genes using the MirTarget program has been completed. For the binding sites of the miR-619-5p the hybridization free energy of the bonds was equal to 100% of the maximum potential free energy. The mRNAs of 201 human genes have complete complementary binding sites of miR-619-5p in the 3’UTR (214 sites), CDS (3 sites), and 5’UTR (4 sites). The mRNAs of CATAD1, ICA1L, GK5, POLH, and PRR11 genes have six miR-619-5p binding sites, and the mRNAs of OPA3 and CYP20A1 genes have eight and ten binding sites, respectively. All of these miR-619-5p binding sites are located in the 3’UTRs. The miR-619-5p binding site in the 5’UTR of mRNA of human USP29 gene is found in the mRNAs of orthologous genes of primates. Binding sites of miR-619-5p in the coding regions of mRNAs of C8H8orf44, C8orf44, and ISY1 genes encode the WLMPVIP oligopeptide, which is present in the orthologous proteins. Binding sites of miR-619-5p in the mRNAs of transcription factor genes ZNF429 and ZNF429 encode the AHACNP oligopeptide in another reading frame. Binding sites of miR-619-5p in the 3’UTRs of all human target genes are also present in the 3’UTRs of orthologous genes of mammals. The completely complementary binding sites for miR-619-5p are conservative in the orthologous mammalian genes.ConclusionsThe majority of miR-619-5p binding sites are located in the 3’UTRs but some genes have miRNA binding sites in the 5’UTRs of mRNAs. Several genes have binding sites for miRNAs in the CDSs that are read in different open reading frames. Identical nucleotide sequences of binding sites encode different amino acids in different proteins. The binding sites of miR-619-5p in 3’UTRs, 5’UTRs and CDSs are conservative in the orthologous mammalian genes.

Segmental duplications and evolutionary acquisition of UV damage response in the SPATA31 gene family of primates and humans

BackgroundSegmental duplications are an abundant source for novel gene functions and evolutionary adaptations. This mechanism of generating novelty was very active during the evolution of primates particularly in the human lineage. Here, we characterize the evolution and function of the SPATA31 gene family (former designation FAM75A), which was previously shown to be among the gene families with the strongest signal of positive selection in hominoids. The mouse homologue for this gene family is a single copy gene expressed during spermatogenesis.ResultsWe show that in primates, the SPATA31 gene duplicated into SPATA31A and SPATA31C types and broadened the expression into many tissues. Each type became further segmentally duplicated in the line towards humans with the largest number of full-length copies found for SPATA31A in humans. Copy number estimates of SPATA31A based on digital PCR show an average of 7.5 with a range of 5–11 copies per diploid genome among human individuals. The primate SPATA31 genes also acquired new protein domains that suggest an involvement in UV response and DNA repair. We generated antibodies and show that the protein is re-localized from the nucleolus to the whole nucleus upon UV-irradiation suggesting a UV damage response. We used CRISPR/Cas mediated mutagenesis to knockout copies of the gene in human primary fibroblast cells. We find that cell lines with reduced functional copies as well as naturally occurring low copy number HFF cells show enhanced sensitivity towards UV-irradiation.ConclusionThe acquisition of new SPATA31 protein functions and its broadening of expression may be related to the evolution of the diurnal life style in primates that required a higher UV tolerance. The increased segmental duplications in hominoids as well as its fast evolution suggest the acquisition of further specific functions particularly in humans.

A quick and robust MHC typing method for free-ranging and captive primate species

Gene products of the major histocompatibility complex (MHC) of human and non-human primates play a crucial role in adaptive immunity, and most of the relevant genes not only show a high degree of variability (polymorphism) but also copy number variation (CNV) is observed. Due to this diversity, MHC proteins influence the capability of individuals to cope with various pathogens. MHC and/or MHC-linked gene products such as odorant receptor genes are thought to influence mate choice and reproductive success. Therefore, MHC typing of wild and captive primate populations is considered to be useful in conservation biology, which is, however, often hampered by the need of invasive and time-consuming methods. All intact Mhc-DRB genes in primates appear to possess a complex and highly divergent microsatellite, DRB-STR. A panel of 154 pedigreed olive baboons (Papio anubis) was examined for their DRB content by DRB-STR analysis of genomic DNA. Using the same methodology on DNA of feces samples, DRB variability of a silvery gibbon population (Hylobates moloch) (N = 24), an endangered species, could successfully be studied. In both species, length determination of the DRB-STR resulted in the definition of unique genotyping patterns that appeared to be specific for a certain chromosome. Moreover, the different STR lengths were shown to segregate with the allelic variation of the respective gene. The results obtained expand data gained previously on DRB-STR typing in macaques, great apes, and humans and strengthen the conclusion that this protocol is applicable in molecular ecology, conservation biology, and colony management, especially of endangered primate species.

Pharmacoinformatics, Adaptive Evolution, and Elucidation of Six Novel Compounds for Schizophrenia Treatment by Targeting DAOA (G72) Isoforms.

Studies on Schizophrenia so far reveal a complex picture of neurological malfunctioning reported to be strongly associated with DAOA. Detailed sequence analyses proved DAOA as a primate specific gene having conserved gene desert region on both upstream and downstream region. The analyses of 10 MB chromosomal region of primates, birds, rodents, and reptiles having DAOA evidenced the conserved part in primates and in the rest of species, while DAOA is only present in primates. DAOA has four isoforms having one interaction partner DAO. Protein-protein analyses of four DAOA isoforms with DAO were performed individually and find potential interacting residues computationally. It was observed that molecular docking of approved FDA drugs revealed efficient results but there was no common drug with effective binding to all DAOA isoforms. Library of compounds was constructed by virtual screening of 2D similarity search against recommended SZ drugs in conjunction with their physiochemical properties. Molecular docking resulted in six novel compounds exhibiting maximum binding affinity with selected four DAOA isoforms. However not the entire schizophrenic population responds to the single drug and interestingly in this study six novel compounds having promising results and same binding site to that DAOA that may be used to interact with DAO against four DAOA isoforms were observed.

Evolution: Genomic remodelling in the primate brain.

In many mammals, the gene Ostn is expressed in muscles and bones. The discovery that the primate OSTN gene has been repurposed to also act in neurons provides clues to how humans evolved their cognitive abilities. See Article p.242 Much of the research on the gene expression networks that drive brain development has been performed in mice. Relatively little is known about how expression networks in other animal groups—particularly primates, in which the cerebral cortex is expanded—might differ from the mouse model. Here, Michael Greenberg and colleagues identify a non-neuronal secreted factor in mice, Osteocrin, that may have been re-purposed evolutionarily as a neuronal development gene in primates. Osteocrin is specifically expressed in the neocortex of the humans and macaques. In mice it is enriched in bone and muscle tissues, but not in the brain.

Differences in Toxicological and Pharmacological Responses Mediated by Polymorphic Cytochromes P450 and Related Drug-Metabolizing Enzymes.

Research over the past 30 years has elucidated the roles of polymorphic human liver cytochrome P450 (P450) enzymes associated with toxicological and/or pharmacological actions. Thalidomide exerts its various pharmacological and toxic actions in primates through multiple mechanisms, including nonspecific modification of many protein networks after bioactivation by autoinduced human P450 enzymes. To overcome species differences between rodents, currently, nonhuman primates and/or mouse models with transplanted human hepatocytes are used. Interindividual variability of P450-dependent drug clearances in cynomolgus monkeys and common marmosets is partly accounted for by polymorphic P450 variants and/or aging, just as it is in humans with increased prevalence of polypharmacy. Genotyping of P450 genes in nonhuman primates would be beneficial before and/or after drug metabolism and toxicity testing and evaluation as well in humans. Genome-wide association studies in humans have been rapidly advanced; however, unique whole-gene deletion of P450 2A6 was subsequently developed to cover nicotine-related lung cancer risk study. Regarding polypharmacy, toxicological research should generally be aimed at identifying the risk of adverse drug events following specific potential drug exposures by examining single or multiple metabolic pathways involving single or multiple drug-metabolizing enzymes. Current and next-generation research of drug metabolism and disposition resulting in drug toxicity would be addressed under advanced knowledge of polymorphic and age-related intra- and/or interspecies differences of drug-metabolizing enzymes. In the near future, humanized animal models combining transplanted hepatocytes and a humanized immune system may be available to study human immune reactions caused by human-type drug metabolites. Such sophisticated models should provide preclinical predictions of human drug metabolism and potential toxicity.

Natural mutations in IFITM3 modulate post-translational regulation and toggle antiviral specificity.

The interferon‐induced transmembrane (IFITM) proteins protect host cells from diverse virus infections. IFITM proteins also incorporate into HIV‐1 virions and inhibit virus fusion and cell‐to‐cell spread, with IFITM3 showing the greatest potency. Here, we report that amino‐terminal mutants of IFITM3 preventing ubiquitination and endocytosis are more abundantly incorporated into virions and exhibit enhanced inhibition of HIV‐1 fusion. An analysis of primate genomes revealed that IFITM3 is the most ancient antiviral family member of the IFITM locus and has undergone a repeated duplication in independent host lineages. Some IFITM3 genes in nonhuman primates, including those that arose following gene duplication, carry amino‐terminal mutations that modify protein localization and function. This suggests that “runaway” IFITM3 variants could be selected for altered antiviral activity. Furthermore, we show that adaptations in IFITM3 result in a trade‐off in antiviral specificity, as variants exhibiting enhanced activity against HIV‐1 poorly restrict influenza A virus. Overall, we provide the first experimental evidence that diversification of IFITM3 genes may boost the antiviral coverage of host cells and provide selective functional advantages.

Species specific exome probes reveal new insights in positively selected genes in nonhuman primates.

Nonhuman primates (NHP) are important biomedical animal models for the study of human disease. Of these, the most widely used models in biomedical research currently are from the genus Macaca. However, evolutionary genetic divergence between human and NHP species makes human-based probes inefficient for the capture of genomic regions of NHP for sequencing and study. Here we introduce a new method to resequence the exome of NHP species by a designed capture approach specifically targeted to the NHP, and demonstrate its superior performance on four NHP species or subspecies. Detailed investigation on biomedically relevant genes demonstrated superior capture by the new approach. We identified 28 genes that appeared to be pseudogenized and inactivated in macaque. Finally, we identified 187 genes showing strong evidence for positive selection across all branches of the primate phylogeny including many novel findings.

Primate genotyping via high resolution melt analysis: rapid and reliable identification of color vision status in wild lemurs

Analyses of genetic polymorphisms can aid our understanding of intra- and interspecific variation in primate sociality, ecology, and behavior. Studies of primate opsin genes are prime examples of this, as single nucleotide variants (SNVs) in the X-linked opsin gene underlie variation in color vision. For primate species with polymorphic trichromacy, genotyping opsin SNVs can generally indicate whether individual primates are red-green color-blind (denoted homozygous M or homozygous L) or have full trichromatic color vision (heterozygous ML). Given the potential influence of color vision on behavior and fitness, characterizing the color vision status of study subjects is becoming commonplace for many primate field projects. Such studies traditionally involve a multi-step sequencing-based method that can be costly and time-consuming. Here we present a new reliable, rapid, and relatively inexpensive method for characterizing color vision in primate populations using high resolution melt analysis (HRMA). Using lemurs as a case study, we characterized variation at exons 3 and/or 5 of the X-linked opsin gene for 87 individuals representing nine species. We scored opsin genotypes and color vision status using both traditional sequencing-based methods as well as our novel melting-curve based HRMA protocol. For each species, the melting curves of varying genotypes (homozygous M, homozygous L, heterozygous ML) differed in melting temperature and/or shape. Melting curves for each sample were consistent across replicates, and genotype-specific melting curves were consistent across DNA sources (blood vs. feces). We show that opsin genotypes can be quickly and reliably scored using HRMA once lab-specific reference curves have been developed based on known genotypes. Although the protocol presented here focuses on genotyping lemur opsin loci, we also consider the larger potential for applying this approach to various types of genetic studies of primate populations.

Conservation and Innovation of APOBEC3A Restriction Functions during Primate Evolution.

LINE-1 (long interspersed element-1) retroelements are the only active autonomous endogenous retroelements in human genomes. Their retrotransposition activity has created close to 50% of the current human genome. Due to the apparent costs of this proliferation, host genomes have evolved multiple mechanisms to curb LINE-1 retrotransposition. Here, we investigate the evolution and function of the LINE-1 restriction factor APOBEC3A, a member of the APOBEC3 cytidine deaminase gene family. We find that APOBEC3A genes have evolved rapidly under diversifying selection in primates, suggesting changes in APOBEC3A have been recurrently selected in a host-pathogen "arms race." Nonetheless, in contrast to previous reports, we find that the LINE-1 restriction activity of APOBEC3A proteins has been strictly conserved throughout simian primate evolution in spite of its pervasive diversifying selection. Based on these results, we conclude that LINE-1s have not driven the rapid evolution of APOBEC3A in primates. In contrast to this conserved LINE-1 restriction, we find that a subset of primate APOBEC3A genes have enhanced antiviral restriction. We trace this gain of antiviral restriction in APOBEC3A to the common ancestor of a subset of Old World monkeys. Thus, APOBEC3A has not only maintained its LINE-1 restriction ability, but also evolved a gain of antiviral specificity against other pathogens. Our findings suggest that while APOBEC3A has evolved to restrict additional pathogens, only those adaptive amino acid changes that leave LINE-1 restriction unperturbed have been tolerated.