Moving primate genomics beyond the chimpanzee genome

Diagnosis and treatment of a poorly differentiated carcinoma in a male chimpanzee (Pan troglodytes)-A case study

This study reports the occurrence of a poorly differentiated carcinoma in a captive-born 28 year-old male chimpanzee (Pan troglodytes) who has a familial history of cancer. Pathological findings, surgical interventions, and experimental treatments are discussed.

High Satellite Repeat Turnover in Great Apes Studied with Short- and Long-Read Technologies

Satellite repeats are a structural component of centromeres and telomeres, and in some instances, their divergence is known to drive speciation. Due to their highly repetitive nature, satellite sequences have been understudied and underrepresented in genome assemblies. To investigate their turnover in great apes, we studied satellite repeats of unit sizes up to 50 bp in human, chimpanzee, bonobo, gorilla, and Sumatran and Bornean orangutans, using unassembled short and long sequencing reads. The density of satellite repeats, as identified from accurate short reads (Illumina), varied greatly among great ape genomes. These were dominated by a handful of abundant repeated motifs, frequently shared among species, which formed two groups: 1) the (AATGG)n repeat (critical for heat shock response) and its derivatives; and 2) subtelomeric 32-mers involved in telomeric metabolism. Using the densities of abundant repeats, individuals could be classified into species. However, clustering did not reproduce the accepted species phylogeny, suggesting rapid repeat evolution. Several abundant repeats were enriched in males versus females; using Y chromosome assemblies or Fluorescent In Situ Hybridization, we validated their location on the Y. Finally, applying a novel computational tool, we identified many satellite repeats completely embedded within long Oxford Nanopore and Pacific Biosciences reads. Such repeats were up to 59 kb in length and consisted of perfect repeats interspersed with other similar sequences. Our results based on sequencing reads generated with three different technologies provide the first detailed characterization of great ape satellite repeats, and open new avenues for exploring their functions.

Superoxide dismutase 1 is positively selected to minimize protein aggregation in great apes

Positive (adaptive) selection has recently been implied in human superoxide dismutase 1 (SOD1), a highly abundant antioxidant protein with energy signaling and antiaging functions, one of very few examples of direct selection on a human protein product (exon); the molecular drivers of this selection are unknown. We mapped 30 extant SOD1 sequences to the recently established mammalian species tree and inferred ancestors, key substitutions, and signatures of selection during the protein's evolution. We detected elevated substitution rates leading to great apes (Hominidae) at ~1 per 2 million years, significantly higher than in other primates and rodents, although these paradoxically generally evolve much faster. The high evolutionary rate was partly due to relaxation of some selection pressures and partly to distinct positive selection of SOD1 in great apes. We then show that higher stability and net charge and changes at the dimer interface were selectively introduced upon separation from old world monkeys and lesser apes (gibbons). Consequently, human, chimpanzee and gorilla SOD1s have a net charge of -6 at physiological pH, whereas the closely related gibbons and macaques have -3. These features consistently point towards selection against the malicious aggregation effects of elevated SOD1 levels in long-living great apes. The findings mirror the impact of human SOD1 mutations that reduce net charge and/or stability and cause ALS, a motor neuron disease characterized by oxidative stress and SOD1 aggregates and triggered by aging. Our study thus marks an example of direct selection for a particular chemical phenotype (high net charge and stability) in a single human protein with possible implications for the evolution of aging.

On the apparent rarity of epithelial cancers in captive chimpanzees

Malignant neoplasms arising from epithelial cells are called carcinomas. Such cancers are diagnosed in about one in three humans in 'developed' countries, with the most common sites affected being lung, breast, prostate, colon, ovary and pancreas. By contrast, carcinomas are said to be rare in captive chimpanzees, which share more than 99% protein sequence homology with humans (and possibly in other related 'great apes'-bonobos, gorillas and orangutans). Simple ascertainment bias is an unlikely explanation, as these nonhuman hominids are recipients of excellent veterinary care in research facilities and zoos, and are typically subjected to necropsies when they die. In keeping with this notion, benign tumours and cancers that are less common in humans are well documented in this population. In this brief overview, we discuss other possible explanations for the reported rarity of carcinomas in our closest evolutionary cousins, including inadequacy of numbers surveyed, differences in life expectancy, diet, genetic susceptibility, immune responses or their microbiomes, and other potential environmental factors. We conclude that while relative carcinoma risk is a likely difference between humans and chimpanzees (and possibly other 'great apes'), a more systematic survey of available data is required for validation of this claim.

Retrophylogenomics place tarsiers on the evolutionary branch of anthropoids

One of the most disputed issues in primate evolution and thus of our own primate roots, is the phylogenetic position of the Southeast Asian tarsier. While much molecular data indicate a basal place in the primate tree shared with strepsirrhines (prosimian monophyly hypothesis), data also exist supporting either an earlier divergence in primates (tarsier-first hypothesis) or a close relationship with anthropoid primates (Haplorrhini hypothesis). The use of retroposon insertions embedded in the Tarsius genome afforded us the unique opportunity to directly test all three hypotheses via three pairwise genome alignments. From millions of retroposons, we found 104 perfect orthologous insertions in both tarsiers and anthropoids to the exclusion of strepsirrhines, providing conflict-free evidence for the Haplorrhini hypothesis, and none supporting either of the other two positions. Thus, tarsiers are clearly the sister group to anthropoids in the clade Haplorrhini.

The evolution of embryo implantation

Embryo implantation varies widely in placental mammals. We review this variation in mammals with a special focus on two features: the depth of implantation and embryonic diapause. We discuss the two major types of implantation depth, superficial and interstitial, and map this character on a well-resolved molecular phylogenetic tree of placental mammals. We infer that relatively deep interstitial implantation has independently evolved at least eight times within placental mammals. Moreover, the superficial type of implantation represents the ancestral state for placental mammals. In addition, we review the genes involved in various phases of implantation, and suggest a future direction in investigating the molecular evolution of implantation-related genes.

Evolutionary Modeling of Genotype-Phenotype Associations, and Application to Primate Coding and Non-coding mtDNA Rate Variation

Variation in substitution rates across a phylogeny can be indicative of shifts in the evolutionary dynamics of a protein or non-protein coding regions. One way to understand these signals is to seek the phenotypic correlates of rate variation. Here, we extended a previously published likelihood method designed to detect evolutionary associations between genotypic evolutionary rate and phenotype over a phylogeny. In simulation with two discrete categories of phenotype, the method has a low false-positive rate and detects greater than 80% of true-positives with a tree length of three or greater and a three-fold or greater change in substitution rate given the phenotype. In addition, we successfully extend the test from two to four phenotype categories and evaluated its performance. We then applied the method to two major hypotheses for rate variation in the mitochondrial genome of primates-longevity and generation time as well as body mass which is correlated with many aspects of life history-using three categories of phenotype through discretization of continuous values. Similar to previous results for mammals, we find that the majority of mitochondrial protein-coding genes show associations consistent with the longevity and body mass predictions and that the predominant signal of association comes from the third codon position. We also found a significant association between maximum lifespan and the evolutionary rate of the control region of the mtDNA. In contrast, 24 protein-coding genes from the nuclear genome do not show a consistent pattern of association, which is inconsistent with the generation time hypothesis. These results show the extended method can robustly identify genotype-phenotype associations up to at least four phenotypic categories, and demonstrate the successful application of the method to study factors affecting neutral evolutionary rate in protein-coding and non-coding loci.

Adrenal androgen production in catarrhine primates and the evolution of adrenarche

Adrenarche is a developmental event involving differentiation of the adrenal gland and production of adrenal androgens, and has been hypothesized to play a role in the extension of the preadolescent phase of human ontogeny. It remains unclear whether any nonhuman primate species shows a similar suite of endocrine, biochemical, and morphological changes as are encompassed by human adrenarche. Here, we report serum concentrations of the adrenal androgens dehydroepiandrosterone (DHEA) and dehydroepiandrosterone sulfate (DHEAS) measured in 698 cross-sectional and mixed longitudinal serum samples from catarrhine primates ranging from 0.6 to 47 years of age. DHEAS in Pan is most similar to that of humans in both age-related pattern and absolute levels, and a transient early increase appears to be present in Gorilla. DHEA levels are highest in Cercocebus, Cercopithecus, and Macaca. We also tested for evidence of adaptive evolution in six genes that code for proteins involved in DHEA/S synthesis. Our genetic analyses demonstrate the protein-coding regions of these genes are highly conserved among sampled primates. We describe a tandem gene duplication event probably mediated by a retrotransposon that resulted in two 3-β-hydroxysteroid dehydrogenase/Delta 5-Delta 4 genes (HSD3B1 and HSD3B2) with tissue specific functions in catarrhines. In humans, HSD3B2 is expressed primarily in the adrenals, ovary, and testis, while HSD3B1 is expressed in the placenta. Taken together, our findings suggest that while adrenarche has been suggested to be unique to hominoids, the evolutionary roots for this developmental stage are more ancient.

Spontaneous abortion and preterm labor and delivery in nonhuman primates: evidence from a captive colony of chimpanzees (Pan troglodytes)

BACKGROUND

Preterm birth is a leading cause of perinatal mortality, yet the evolutionary history of this obstetrical syndrome is largely unknown in nonhuman primate species.

METHODOLOGY/PRINCIPAL FINDINGS

We examined the length of gestation during pregnancies that occurred in a captive chimpanzee colony by inspecting veterinary and behavioral records spanning a total of thirty years. Upon examination of these records we were able to confidently estimate gestation length for 93 of the 97 (96%) pregnancies recorded at the colony. In total, 78 singleton gestations resulted in live birth, and from these pregnancies we estimated the mean gestation length of normal chimpanzee pregnancies to be 228 days, a finding consistent with other published reports. We also calculated that the range of gestation in normal chimpanzee pregnancies is approximately forty days. Of the remaining fifteen pregnancies, only one of the offspring survived, suggesting viability for chimpanzees requires a gestation of approximately 200 days. These fifteen pregnancies constitute spontaneous abortions and preterm deliveries, for which the upper gestational age limit was defined as 2 SD from the mean length of gestation (208 days).

CONCLUSIONS/SIGNIFICANCE

The present study documents that preterm birth occurred within our study population of captive chimpanzees. As in humans, pregnancy loss is not uncommon in chimpanzees, In addition, our findings indicate that both humans and chimpanzees show a similar range of normal variation in gestation length, suggesting this was the case at the time of their last common ancestor (LCA). Nevertheless, our data suggest that whereas chimpanzees' normal gestation length is ∼20-30 days after reaching viability, humans' normal gestation length is approximately 50 days beyond the estimated date of viability without medical intervention. Future research using a comparative evolutionary framework should help to clarify the extent to which mechanisms at work in normal and preterm parturition are shared in these species.

Conversion events in gene clusters

BACKGROUND

Gene clusters containing multiple similar genomic regions in close proximity are of great interest for biomedical studies because of their associations with inherited diseases. However, such regions are difficult to analyze due to their structural complexity and their complicated evolutionary histories, reflecting a variety of large-scale mutational events. In particular, conversion events can mislead inferences about the relationships among these regions, as traced by traditional methods such as construction of phylogenetic trees or multi-species alignments.

RESULTS

To correct the distorted information generated by such methods, we have developed an automated pipeline called CHAP (Cluster History Analysis Package) for detecting conversion events. We used this pipeline to analyze the conversion events that affected two well-studied gene clusters (α-globin and β-globin) and three gene clusters for which comparative sequence data were generated from seven primate species: CCL (chemokine ligand), IFN (interferon), and CYP2abf (part of cytochrome P450 family 2). CHAP is freely available at http://www.bx.psu.edu/miller_lab.

CONCLUSIONS

These studies reveal the value of characterizing conversion events in the context of studying gene clusters in complex genomes.

Complex transcriptional control of the AZFa gene DDX3Y in human testis

The human DEAD-box Y (DBY) RNA helicase (aka DDX3Y) gene is thought to be the major azoospermia factor a (AZFa) gene in proximal Yq11. Men with its deletion display no somatic pathologies, but suffer from complete absence of germ cells. Accordingly, DDX3Y protein is expressed only in the germline in spermatogonia, although the transcripts were found in many tissues. Here, we show the complex transcriptional control of a testis-specific DDX3Y transcript class with initiation at different sites upstream of the gene’s open reading frame (5′Untranslated Region; UTR) and with polyadenylation in their proximal 3′UTR. The most distal transcriptional start site (TSS; ∼1 kb upstream) was mapped in MSY2, a Y-specific minisatellite. As this testis-specific 5′UTR was subsequently processed by three alternative splicing events, it has been tentatively designated ‘exon-T’(estis). The MSY2 sequence unit was also found upstream of the mouse Ddx3y gene. However, only after its tandem amplification on the Y chromosome of Platyrrhini (new world monkeys) and Catarrhini (old world monkeys) did MSY2 become part of a novel distal promoter for DDX3Y expression in testis tissue and provides a second transcriptional start site (T-TSS-II) in Catarrhini. We therefore suggest that the development of a novel distal DDX3Y promoter in primates, which is activated only in testis tissue, is probably part of the gene’s germline translation control.

Human and great ape red blood cells differ in plasmalogen levels and composition

Background

Plasmalogens are ether phospholipids required for normal mammalian developmental, physiological, and cognitive functions. They have been proposed to act as membrane antioxidants and reservoirs of polyunsaturated fatty acids as well as influence intracellular signaling and membrane dynamics. Plasmalogens are particularly enriched in cells and tissues of the human nervous, immune, and cardiovascular systems. Humans with severely reduced plasmalogen levels have reduced life spans, abnormal neurological development, skeletal dysplasia, impaired respiration, and cataracts. Plasmalogen deficiency is also found in the brain tissue of individuals with Alzheimer disease.

Results

In a human and great ape cohort, we measured the red blood cell (RBC) levels of the most abundant types of plasmalogens. Total RBC plasmalogen levels were lower in humans than bonobos, chimpanzees, and gorillas, but higher than orangutans. There were especially pronounced cross-species differences in the levels of plasmalogens with a C16:0 moiety at the sn-1 position. Humans on Western or vegan diets had comparable total RBC plasmalogen levels, but the latter group showed moderately higher levels of plasmalogens with a C18:1 moiety at the sn-1 position. We did not find robust sex-specific differences in human or chimpanzee RBC plasmalogen levels or composition. Furthermore, human and great ape skin fibroblasts showed only modest differences in peroxisomal plasmalogen biosynthetic activity. Human and chimpanzee microarray data indicated that genes involved in plasmalogen biosynthesis show cross-species differential expression in multiple tissues.

Conclusion

We propose that the observed differences in human and great ape RBC plasmalogens are primarily caused by their rates of biosynthesis and/or turnover. Gene expression data raise the possibility that other human and great ape cells and tissues differ in plasmalogen levels. Based on the phenotypes of humans and rodents with plasmalogen disorders, we propose that cross-species differences in tissue plasmalogen levels could influence organ functions and processes ranging from cognition to reproduction to aging.

Adaptive evolution of four microcephaly genes and the evolution of brain size in anthropoid primates

The anatomical basis and adaptive function of the expansion in primate brain size have long been studied; however, we are only beginning to understand the genetic basis of these evolutionary changes. Genes linked to human primary microcephaly have received much attention as they have accelerated evolutionary rates along lineages leading to humans. However, these studies focus narrowly on apes, and the link between microcephaly gene evolution and brain evolution is disputed. We analyzed the molecular evolution of four genes associated with microcephaly (ASPM, CDK5RAP2, CENPJ, MCPH1) across 21 species representing all major clades of anthropoid primates. Contrary to prevailing assumptions, positive selection was not limited to or intensified along the lineage leading to humans. In fact we show that all four loci were subject to positive selection across the anthropoid primate phylogeny. We developed clearly defined hypotheses to explicitly test if selection on these loci was associated with the evolution of brain size. We found positive relationships between both CDK5RAP2 and ASPM and neonatal brain mass and somewhat weaker relationships between these genes and adult brain size. In contrast, there is no evidence linking CENPJ and MCPH1 to brain size evolution. The stronger association of ASPM and CDK5RAP2 evolution with neonatal brain size than with adult brain size is consistent with these loci having a direct effect on prenatal neuronal proliferation. These results suggest that primate brain size may have at least a partially conserved genetic basis. Our results contradict a previous study that linked adaptive evolution of ASPM to changes in relative cortex size; however, our analysis indicates that this conclusion is not robust. Our finding that the coding regions of two widely expressed loci has experienced pervasive positive selection in relation to a complex, quantitative developmental phenotype provides a notable counterexample to the commonly asserted hypothesis that cis-regulatory regions play a dominant role in phenotypic evolution.

Evolutionary genetics in wild primates: combining genetic approaches with field studies of natural populations

Ecological and evolutionary studies of wild primates hold important keys to understanding both the shared characteristics of primate biology and the genetic and phenotypic differences that make specific lineages, including our own, unique. Although complementary genetic research on nonhuman primates has long been of interest, recent technological and methodological advances now enable functional and population genetic studies in an unprecedented manner. In the past several years, novel genetic data sets have revealed new information about the demographic history of primate populations and the genetics of adaptively important traits. In combination with the rich history of behavioral, ecological, and physiological work on natural primate populations, genetic approaches promise to provide a compelling picture of primate evolution in the past and in the present day.

Can chimpanzee biology highlight human origin and evolution?

The closest living relatives of humans are their chimpanzee/bonobo (Pan) sister species, members of the same subfamily "Homininae". This classification is supported by over 50 years of research in the fields of chimpanzee cultural diversity, language competency, genomics, anatomy, high cognition, psychology, society, self-consciousness and relation to others, tool use/production, as well as Homo level emotions, symbolic competency, memory recollection, complex multifaceted problem-solving capabilities, and interspecies communication. Language competence and symbolism can be continuously bridged from chimpanzee to man. Emotions, intercommunity aggression, body language, gestures, facial expressions, and vocalization of intonations seem to parallel between the sister taxa Homo and Pan. The shared suite of traits between Pan and Homo genus demonstrated in this article integrates old and new information on human-chimpanzee evolution, bilateral informational and cross-cultural exchange, promoting the urgent need for Pan cultures in the wild to be protected, as they are part of the cultural heritage of mankind. Also, we suggest that bonobos, Pan paniscus, based on shared traits with Australopithecus, need to be included in Australopithecine's subgenus, and may even represent living-fossil Australopithecines. Unfolding bonobo and chimpanzee biology highlights our common genetic and cultural evolutionary origins.

Comparative anatomy of the locus coeruleus in humans and nonhuman primates

The locus coeruleus (LC) is a dense cluster of neurons that projects axons throughout the neuroaxis and is located in the rostral pontine tegmentum extending from the level of the inferior colliculus to the motor nucleus of the trigeminal nerve. LC neurons are lost in the course of several neurodegenerative disorders, including Alzheimer's and Parkinson's diseases. In this study we used Nissl staining and tyrosine hydroxylase (TH) immunoreactivity to compare the human LC with that of closely related primate species, including great and lesser apes, and macaque monkeys. TH catalyzes the initial and rate-limiting step in catecholamine biosynthesis. The number of TH-immunoreactive (TH-ir) neurons was estimated in each species using stereologic methods. In the LC of humans the mean total number of TH-ir neurons was significantly higher compared to the other primates. Because the total number of TH-ir neurons in the LC was highly correlated with the species mean volume of the medulla oblongata, cerebellum, and neocortical gray matter, we conclude that much of the observed phylogenetic variation can be explained by anatomical scaling. Notably, the total number of LC neurons in humans was most closely predicted by the nonhuman allometric scaling relationship relative to medulla size, whereas the number of LC neurons in humans was considerably lower than predicted according to neocortex and cerebellum volume.

Colloquium paper: phylogenomic evidence of adaptive evolution in the ancestry of humans

In Charles Darwin's tree model for life's evolution, natural selection adaptively modifies newly arisen species as they branch apart from their common ancestor. In accord with this Darwinian concept, the phylogenomic approach to elucidating adaptive evolution in genes and genomes in the ancestry of modern humans requires a well supported and well sampled phylogeny that accurately places humans and other primates and mammals with respect to one another. For more than a century, first from the comparative immunological work of Nuttall on blood sera and now from comparative genomic studies, molecular findings have demonstrated the close kinship of humans to chimpanzees. The close genetic correspondence of chimpanzees to humans and the relative shortness of our evolutionary separation suggest that most distinctive features of the modern human phenotype had already evolved during our ancestry with chimpanzees. Thus, a phylogenomic assessment of being human should examine earlier stages of human ancestry as well as later stages. In addition, with the availability of a number of mammalian genomes, similarities in phenotype between distantly related taxa should be explored for evidence of convergent or parallel adaptive evolution. As an example, recent phylogenomic evidence has shown that adaptive evolution of aerobic energy metabolism genes may have helped shape such distinctive modern human features as long life spans and enlarged brains in the ancestries of both humans and elephants.

New perspectives on anthropoid origins

Adaptive shifts associated with human origins are brought to light as we examine the human fossil record and study our own genome and that of our closest ape relatives. However, the more ancient roots of many human characteristics are revealed through the study of a broader array of living anthropoids and the increasingly dense fossil record of the earliest anthropoid radiations. Genomic data and fossils of early primates in Asia and Africa clarify relationships among the major clades of primates. Progress in comparative anatomy, genomics, and molecular biology point to key changes in sensory ecology and brain organization that ultimately set the stage for the emergence of the human lineage.

Reconstructing the ups and downs of primate brain evolution: implications for adaptive hypotheses and Homo floresiensis

Background

Brain size is a key adaptive trait. It is often assumed that increasing brain size was a general evolutionary trend in primates, yet recent fossil discoveries have documented brain size decreases in some lineages, raising the question of how general a trend there was for brains to increase in mass over evolutionary time. We present the first systematic phylogenetic analysis designed to answer this question.

Results

We performed ancestral state reconstructions of three traits (absolute brain mass, absolute body mass, relative brain mass) using 37 extant and 23 extinct primate species and three approaches to ancestral state reconstruction: parsimony, maximum likelihood and Bayesian Markov-chain Monte Carlo. Both absolute and relative brain mass generally increased over evolutionary time, but body mass did not. Nevertheless both absolute and relative brain mass decreased along several branches. Applying these results to the contentious case of Homo floresiensis, we find a number of scenarios under which the proposed evolution of Homo floresiensis' small brain appears to be consistent with patterns observed along other lineages, dependent on body mass and phylogenetic position.

Conclusions

Our results confirm that brain expansion began early in primate evolution and show that increases occurred in all major clades. Only in terms of an increase in absolute mass does the human lineage appear particularly striking, with both the rate of proportional change in mass and relative brain size having episodes of greater expansion elsewhere on the primate phylogeny. However, decreases in brain mass also occurred along branches in all major clades, and we conclude that, while selection has acted to enlarge primate brains, in some lineages this trend has been reversed. Further analyses of the phylogenetic position of Homo floresiensis and better body mass estimates are required to confirm the plausibility of the evolution of its small brain mass. We find that for our dataset the Bayesian analysis for ancestral state reconstruction is least affected by inclusion of fossil data suggesting that this approach might be preferable for future studies on other taxa with a poor fossil record.

A small, variable, and irregular killer cell Ig-like receptor locus accompanies the absence of MHC-C and MHC-G in gibbons

The killer cell Ig-like receptors (KIRs) of NK cells recognize MHC class I ligands and function in placental reproduction and immune defense against pathogens. During the evolution of monkeys, great apes, and humans, an ancestral KIR3DL gene expanded to become a diverse and rapidly evolving gene family of four KIR lineages. Characterizing the KIR locus are three framework regions, defining two intervals of variable gene content. By analysis of four KIR haplotypes from two species of gibbon, we find that the smaller apes do not conform to these rules. Although diverse and irregular in structure, the gibbon haplotypes are unusually small, containing only two to five functional genes. Comparison with the predicted ancestral hominoid KIR haplotype indicates that modern gibbon KIR haplotypes were formed by a series of deletion events, which created new hybrid genes as well as eliminating ancestral genes. Of the three framework regions, only KIR3DL3 (lineage V), defining the 5' end of the KIR locus, is present and intact on all gibbon KIR haplotypes. KIR2DL4 (lineage I) defining the central framework region has been a major target for elimination or inactivation, correlating with the absence of its putative ligand, MHC-G, in gibbons. Similarly, the MHC-C-driven expansion of lineage III KIR genes in great apes has not occurred in gibbons because they lack MHC-C. Our results indicate that the selective forces shaping the size and organization of the gibbon KIR locus differed from those acting upon the KIR of other hominoid species.

Inhibitory interneurons of the human prefrontal cortex display conserved evolution of the phenotype and related genes

Inhibitory interneurons participate in local processing circuits, playing a central role in executive cognitive functions of the prefrontal cortex. Although humans differ from other primates in a number of cognitive domains, it is not currently known whether the interneuron system has changed in the course of primate evolution leading to our species. In this study, we examined the distribution of different interneuron subtypes in the prefrontal cortex of anthropoid primates as revealed by immunohistochemistry against the calcium-binding proteins calbindin, calretinin and parvalbumin. In addition, we tested whether genes involved in the specification, differentiation and migration of interneurons show evidence of positive selection in the evolution of humans. Our findings demonstrate that cellular distributions of interneuron subtypes in human prefrontal cortex are similar to other anthropoid primates and can be explained by general scaling rules. Furthermore, genes underlying interneuron development are highly conserved at the amino acid level in primate evolution. Taken together, these results suggest that the prefrontal cortex in humans retains a similar inhibitory circuitry to that in closely related primates, even though it performs functional operations that are unique to our species. Thus, it is likely that other significant modifications to the connectivity and molecular biology of the prefrontal cortex were overlaid on this conserved interneuron architecture in the course of human evolution.

Phylogenomic analyses reveal convergent patterns of adaptive evolution in elephant and human ancestries

Specific sets of brain-expressed genes, such as aerobic energy metabolism genes, evolved adaptively in the ancestry of humans and may have evolved adaptively in the ancestry of other large-brained mammals. The recent addition of genomes from two afrotherians (elephant and tenrec) to the expanding set of publically available sequenced mammalian genomes provided an opportunity to test this hypothesis. Elephants resemble humans by having large brains and long life spans; tenrecs, in contrast, have small brains and short life spans. Thus, we investigated whether the phylogenomic patterns of adaptive evolution are more similar between elephant and human than between either elephant and tenrec lineages or human and mouse lineages, and whether aerobic energy metabolism genes are especially well represented in the elephant and human patterns. Our analyses encompassed approximately 6,000 genes in each of these lineages with each gene yielding extensive coding sequence matches in interordinal comparisons. Each gene's nonsynonymous and synonymous nucleotide substitution rates and dN/dS ratios were determined. Then, from gene ontology information on genes with the higher dN/dS ratios, we identified the more prevalent sets of genes that belong to specific functional categories and that evolved adaptively. Elephant and human lineages showed much slower nucleotide substitution rates than tenrec and mouse lineages but more adaptively evolved genes. In correlation with absolute brain size and brain oxygen consumption being largest in elephants and next largest in humans, adaptively evolved aerobic energy metabolism genes were most evident in the elephant lineage and next most evident in the human lineage.

Phylogenomics of primates and their ancestral populations

Genome assemblies are now available for nine primate species, and large-scale sequencing projects are underway or approved for six others. An explicitly evolutionary and phylogenetic approach to comparative genomics, called phylogenomics, will be essential in unlocking the valuable information about evolutionary history and genomic function that is contained within these genomes. However, most phylogenomic analyses so far have ignored the effects of variation in ancestral populations on patterns of sequence divergence. These effects can be pronounced in the primates, owing to large ancestral effective population sizes relative to the intervals between speciation events. In particular, local genealogies can vary considerably across loci, which can produce biases and diminished power in many phylogenomic analyses of interest, including phylogeny reconstruction, the identification of functional elements, and the detection of natural selection. At the same time, this variation in genealogies can be exploited to gain insight into the nature of ancestral populations. In this Perspective, I explore this area of intersection between phylogenetics and population genetics, and its implications for primate phylogenomics. I begin by "lifting the hood" on the conventional tree-like representation of the phylogenetic relationships between species, to expose the population-genetic processes that operate along its branches. Next, I briefly review an emerging literature that makes use of the complex relationships among coalescence, recombination, and speciation to produce inferences about evolutionary histories, ancestral populations, and natural selection. Finally, I discuss remaining challenges and future prospects at this nexus of phylogenetics, population genetics, and genomics.

Alpha-fetoprotein related gene (ARG): a new member of the albumin gene family that is no longer functional in primates

The serum albumin gene family is comprised of albumin, alpha-fetoprotein, alpha-albumin (afamin), and the more distantly related Vitamin D binding protein. These genes arose from a common ancestor through a series of duplication events, are expressed primarily in the liver and tightly linked in all species where this has been investigated. Here, we describe a fifth member of the albumin gene family that we have named Alpha-fetoprotein Related Gene (ARG) since it exhibits greatest similarity to this family member. ARG is activated in the liver perinatally, but is expressed at very low levels. The ARG gene is present and intact in the mouse, rat, dog and horse genomes. In contrast, the ARG gene in human, chimpanzee, rhesus monkey, and marmoset contains a number of mutations common to all four species, indicating that this gene has been an inactive pseudogene in primates for at least 40 million years. Low expression and aberrant splicing of the ARG gene in the mouse liver suggests that ARG may have less functional significance than other members of the serum albumin gene family even in species where it is still intact.

Much Hype and Many Errors Uncovering Our Earliest Ancestor Little, Brown, New York, 2009. 304 pp. $25.99, C$28.99. ISBN 9780316070089

This account of Darwinius masillae , its discovery, and its importance was rushed into print as part of the hype surrounding the public announcement of the work.

Patterns of macroevolution among Primates inferred from a supermatrix of mitochondrial and nuclear DNA

Here, we present a new primate phylogeny inferred from molecular supermatrix analyses of size 42 kb containing 70% of missing data, and representing 75% of primate species diversity. The supermatrix was analysed using a gene-partitioned maximum likelihood approach to obtain an exhaustive molecular phylogenetic framework. All clades recovered from recent molecular works were upheld in our analyses demonstrating that the presence of missing data did not bias our supermatrix inference. The resulting phylogenetic tree was subsequently dated with a molecular dating method to provide a timescale for speciation events. Results obtained from our relaxed molecular clock analyses concurred with previous works based on the same fossil constraints. The resulting dated tree allowed to infer of macroevolutionary processes among the primates. Shifts in diversification rate and speciation rates were determined using the SymmeTREE method and a birthdeath process. No significant asymmetry was detected for the primate clade, but significant shifts in diversification rate were identified for seven clades: Anthropoidea, Lemuriformes, Lemuridae, Galagidae, Callithrix genus, the Cercopithecinae and Asian Macaca. Comparisons with previous primate supertree results reveal that (i) there was a diversification event at the root of the Lemuriformes, (ii) a higher diversification rate is detected for Cercopithecidae and Anthropoidea and (iii) a shift in diversification is always recovered for Macaca genus. Macroevolutionary inferences and primate divergence dates show that major primate diversification events occurred after the Paleogene, suggesting the extinction of ancient primate lineages.

A fully resolved genus level phylogeny of neotropical primates (Platyrrhini)

There are more than 125 species of extant New World monkeys (Primates: Platyrrhini) found in approximately 15 genera. The phylogenetic relationships of these neotropical primates have been extensively studied from a molecular perspective. While these studies have been successful at inferring many of the relationships within the platyrrhines, key questions remain. The current study provides a framework for using non-genic, non-coding markers in comparative primate phylogenomic studies in species whose genomes are not yet scheduled for complete sequencing. A random genomic shotgun library was generated from the nocturnal Owl monkey Aotus lemurinus. Eleven unlinked, non-coding, non-genic, non-repetitive, nuclear DNA markers derived from this library were sequenced in at least one representative species of every platyrrhine genus. The combined sequence from these markers yielded a 7.7 kb multiple sequence alignment of 22 taxa. We analyzed these markers independently and combined with a 10 kb dataset consisting of "traditional," previously published markers located within or directly adjacent to genes. Parsimony, maximum likelihood, and Bayesian analysis converged on a single topology for the platyrrhine generic relationships. Notably, we confidently inferred that Pitheciidae is the sister taxon to the other two platyrrhine families (Cebidae, Atelidae). This relationship is supported by high values of branch support as well as topology tests. Additionally, Aotus formed a sister taxon to a clade comprising Cebus and Saimiri. With a fully resolved platyrrhine phylogeny in place it is now possible to design and test hypotheses regarding the evolution and diversification of platyrrhine phenotypes and life histories.

Primate phylogenomics: developing numerous nuclear non-coding, non-repetitive markers for ecological and phylogenetic applications and analysis of evolutionary rate variation

Background

Genetic analyses are often limited by the availability of appropriate molecular markers. Markers from neutrally evolving genomic regions may be particularly useful for inferring evolutionary histories because they escape the constraints of natural selection. For the majority of taxa however, obtaining such markers is challenging. Advances in genomics have the potential to alleviate the shortage of neutral markers. Here we present a method to develop numerous markers from putatively neutral regions of primate genomes.

Results

We began with the available whole genome sequences of human, chimpanzee and macaque. Using computational methods, we identified a total of 280 potential amplicons from putatively neutral, non-coding, non-repetitive regions of these genomes. Subsequently we amplified, using experimental methods, many of these amplicons from diverse primate taxa, including a ring-tailed lemur, which is distantly related to the genomic resources. Using a subset of 10 markers, we demonstrate the utility of the developed markers in phylogenetic and evolutionary rate analyses. Particularly, we uncovered substantial evolutionary rate variation among lineages, some of which are previously not reported.

Conclusion

We successfully developed numerous markers from putatively neutral regions of primate genomes using a strategy combining computational and experimental methods. Applying these markers to phylogenetic and evolutionary rate variation analyses exemplifies the utility of these markers. Diverse ecological and evolutionary analyses will benefit from these markers. Importantly, methods similar to those presented here can be applied to other taxa in the near future.

Ancient mtDNA genetic variants modulate mtDNA transcription and replication

Although the functional consequences of mitochondrial DNA (mtDNA) genetic backgrounds (haplotypes, haplogroups) have been demonstrated by both disease association studies and cell culture experiments, it is not clear which of the mutations within the haplogroup carry functional implications and which are "evolutionary silent hitchhikers". We set forth to study the functionality of haplogroup-defining mutations within the mtDNA transcription/replication regulatory region by in vitro transcription, hypothesizing that haplogroup-defining mutations occurring within regulatory motifs of mtDNA could affect these processes. We thus screened >2500 complete human mtDNAs representing all major populations worldwide for natural variation in experimentally established protein binding sites and regulatory regions comprising a total of 241 bp in each mtDNA. Our screen revealed 77/241 sites showing point mutations that could be divided into non-fixed (57/77, 74%) and haplogroup/sub-haplogroup-defining changes (i.e., population fixed changes, 20/77, 26%). The variant defining Caucasian haplogroup J (C295T) increased the binding of TFAM (Electro Mobility Shift Assay) and the capacity of in vitro L-strand transcription, especially of a shorter transcript that maps immediately upstream of conserved sequence block 1 (CSB1), a region associated with RNA priming of mtDNA replication. Consistent with this finding, cybrids (i.e., cells sharing the same nuclear genetic background but differing in their mtDNA backgrounds) harboring haplogroup J mtDNA had a >2 fold increase in mtDNA copy number, as compared to cybrids containing haplogroup H, with no apparent differences in steady state levels of mtDNA-encoded transcripts. Hence, a haplogroup J regulatory region mutation affects mtDNA replication or stability, which may partially account for the phenotypic impact of this haplogroup. Our analysis thus demonstrates, for the first time, the functional impact of particular mtDNA haplogroup-defining control region mutations, paving the path towards assessing the functionality of both fixed and un-fixed genetic variants in the mitochondrial genome.

Evolutionary analysis of the highly dynamic CHEK2 duplicon in anthropoids

BACKGROUND

Segmental duplications (SDs) are euchromatic portions of genomic DNA (> or = 1 kb) that occur at more than one site within the genome, and typically share a high level of sequence identity (>90%). Approximately 5% of the human genome is composed of such duplicated sequences. Here we report the detailed investigation of CHEK2 duplications. CHEK2 is a multiorgan cancer susceptibility gene encoding a cell cycle checkpoint kinase acting in the DNA-damage response signalling pathway. The continuous presence of the CHEK2 gene in all eukaryotes and its important role in maintaining genome stability prompted us to investigate the duplicative evolution and phylogeny of CHEK2 and its paralogs during anthropoid evolution.

RESULTS

To study CHEK2 duplicon evolution in anthropoids we applied a combination of comparative FISH and in silico analyses. Our comparative FISH results with a CHEK2 fosmid probe revealed the single-copy status of CHEK2 in New World monkeys, Old World monkeys and gibbons. Whereas a single CHEK2 duplication was detected in orangutan, a multi-site signal pattern indicated a burst of duplication in African great apes and human. Phylogenetic analysis of paralogous and ancestral CHEK2 sequences in human, chimpanzee and rhesus macaque confirmed this burst of duplication, which occurred after the radiation of orangutan and African great apes. In addition, we used inter-species quantitative PCR to determine CHEK2 copy numbers. An amplification of CHEK2 was detected in African great apes and the highest CHEK2 copy number of all analysed species was observed in the human genome. Furthermore, we detected variation in CHEK2 copy numbers within the analysed set of human samples.

CONCLUSION

Our detailed analysis revealed the highly dynamic nature of CHEK2 duplication during anthropoid evolution. We determined a burst of CHEK2 duplication after the radiation of orangutan and African great apes and identified the highest CHEK2 copy number in human. In conclusion, our analysis of CHEK2 duplicon evolution revealed that SDs contribute to inter-species variation. Furthermore, our qPCR analysis led us to presume CHEK2 copy number variation in human, and molecular diagnostics of the cancer susceptibility gene CHEK2 inside the duplicated region might be hampered by the individual-specific set of duplicons.

Evolution of a behavior-linked microsatellite-containing element in the 5' flanking region of the primate AVPR1A gene

Background

The arginine vasopressin V1a receptor (V1aR) modulates social cognition and behavior in a wide variety of species. Variation in a repetitive microsatellite element in the 5' flanking region of the V1aR gene (AVPR1A) in rodents has been associated with variation in brain V1aR expression and in social behavior. In humans, the 5' flanking region of AVPR1A contains a tandem duplication of two ~350 bp, microsatellite-containing elements located approximately 3.5 kb upstream of the transcription start site. The first block, referred to as DupA, contains a polymorphic (GT)₂₅ microsatellite; the second block, DupB, has a complex (CT)₄-(TT)-(CT)₈-(GT)₂₄ polymorphic motif, known as RS3. Polymorphisms in RS3 have been associated with variation in sociobehavioral traits in humans, including autism spectrum disorders. Thus, evolution of these regions may have contributed to variation in social behavior in primates. We examined the structure of these regions in six ape, six monkey, and one prosimian species.

Results

Both tandem repeat blocks are present upstream of the AVPR1A coding region in five of the ape species we investigated, while monkeys have only one copy of this region. As in humans, the microsatellites within DupA and DupB are polymorphic in many primate species. Furthermore, both single (lacking DupB) and duplicated alleles (containing both DupA and DupB) are present in chimpanzee (Pan troglodytes) populations with allele frequencies of 0.795 and 0.205 for the single and duplicated alleles, respectively, based on the analysis of 47 wild-caught individuals. Finally, a phylogenetic reconstruction suggests two alternate evolutionary histories for this locus.

Conclusion

There is no obvious relationship between the presence of the RS3 duplication and social organization in primates. However, polymorphisms identified in some species may be useful in future genetic association studies. In particular, the presence of both single and duplicated alleles in chimpanzees provides a unique opportunity to assess the functional role of this duplication in contributing to variation in social behavior in primates. While our initial studies show no signs of directional selection on this locus in chimps, pharmacological and genetic association studies support a potential role for this region in influencing V1aR expression and social behavior.

Evolutionary dynamics of segmental duplications from human Y-chromosomal euchromatin/heterochromatin transition regions

Human chromosomal regions enriched in segmental duplications are subject to extensive genomic reorganization. Such regions are particularly informative for illuminating the evolutionary history of a given chromosome. We have analyzed 866 kb of Y-chromosomal non-palindromic segmental duplications delineating four euchromatin/heterochromatin transition regions (Yp11.2/Yp11.1, Yq11.1/Yq11.21, Yq11.23/Yq12, and Yq12/PAR2). Several computational methods were applied to decipher the segmental duplication architecture and identify the ancestral origin of the 41 different duplicons. Combining computational and comparative FISH analysis, we reconstruct the evolutionary history of these regions. Our analysis indicates a continuous process of transposition of duplicated sequences onto the evolving higher primate Y chromosome, providing unique insights into the development of species-specific Y-chromosomal and autosomal duplicons. Phylogenetic sequence comparisons show that duplicons of the human Yp11.2/Yp11.1 region were already present in the macaque-human ancestor as multiple paralogs located predominantly in subtelomeric regions. In contrast, duplicons from the Yq11.1/Yq11.21, Yq11.23/Yq12, and Yq12/PAR2 regions show no evidence of duplication in rhesus macaque, but map to the pericentromeric regions in chimpanzee and human. This suggests an evolutionary shift in the direction of duplicative transposition events from subtelomeric in Old World monkeys to pericentromeric in the human/ape lineage. Extensive chromosomal relocation of autosomal-duplicated sequences from euchromatin/heterochromatin transition regions to interstitial regions as demonstrated on the pygmy chimpanzee Y chromosome support a model in which substantial reorganization and amplification of duplicated sequences may contribute to speciation.

The human progesterone receptor shows evidence of adaptive evolution associated with its ability to act as a transcription factor

The gene encoding the progesterone receptor (PGR) acts as a transcription factor, and participates in the regulation of reproductive processes including menstruation, implantation, pregnancy maintenance, parturition, mammary development, and lactation. Unlike other mammals, primates do not exhibit progesterone withdrawal at the time of parturition. Because progesterone-mediated reproductive features vary among mammals, PGR is an attractive candidate gene for studies of adaptive evolution. Thus, we sequenced the progesterone receptor coding regions in a diverse range of species including apes, Old World monkeys, New World monkeys, prosimian primates, and other mammals. Adaptive evolution occurred on the human and chimpanzee lineages as evidenced by statistically significant increases in nonsynonymous substitution rates compared to synonymous substitution rates. Positive selection was rarely observed in other lineages. In humans, amino acid replacements occurred mostly in a region of the gene that has been shown to have an inhibitory function (IF) on the ability of the progesterone receptor to act as a transcription factor. Moreover, many of the nonsynonymous substitutions in primates occurred in the N-terminus. This suggests that cofactor interaction surfaces might have been altered, resulting in altered progesterone-regulated gene transcriptional effects. Further evidence that the changes conferred an adaptive advantage comes from SNP analysis indicating only one of the IF changes is polymorphic in humans. In chimpanzees, amino acid changes occurred in both the inhibitory and transactivation domains. Positive selection provides the basis for the hypothesis that changes in structure and function of the progesterone receptor during evolution contribute to the diversity of primate reproductive biology, especially in parturition.

The molecular signature of selection underlying human adaptations

In the last decade, advances in human population genetics and comparative genomics have resulted in important contributions to our understanding of human genetic diversity and genetic adaptation. For the first time, we are able to reliably detect the signature of natural selection from patterns of DNA polymorphism. Identifying the effects of natural selection in this way provides a crucial piece of evidence needed to support hypotheses of human adaptation. This review provides a detailed description of the theory and analytical approaches used to detect signatures of natural selection in the human genome. We discuss these methods in relation to four classic human traits--skin color, the Duffy blood group, bitter-taste sensation, and lactase persistence. By highlighting these four traits we are able to discuss the ways in which analyses of DNA polymorphism can lead to inferences regarding past histories of selection. Specifically, we can infer the importance of specific regimes of selection (i.e. directional selection, balancing selection, and purifying selection) in the evolution of a trait because these different types of selection leave different patterns of DNA polymorphism. In addition, we demonstrate how these types of data can be used to estimate the time frame in which selection operated on a trait. As the field has advanced, a general issue that has come to the forefront is how specific demographic events in human history, such as population expansions, bottlenecks, and subdivision of populations, have also left a signature across the genome that can interfere with our detection of the footprint of selection at particular genes. Therefore, we discuss this general problem with respect to the four traits reviewed here, and describe the ways in which the signature of selection can be teased from a background signature of demographic history. Finally, we move from a discussion of analyses of selection motivated by a "candidate-gene" approach, in which a priori information led to the analysis of specific gene, to discussion of "genome-scanning" approaches that are directed at discovering new genes that have been under positive selection. Such scans can be designed to detect those genes that have been positively selected in our divergence from chimpanzees, as well as those genes that have been under selection as human populations have migrated, differentiated, and adapted to specific geographic environments. We predict that both approaches will be applied in the future, enabling a greater insight into human species-wide adaptations, as well as the specific adaptations of human populations.