Primate genomes offer new view of human health and our past.

Sequencing efforts may also aid primate conservation.

Comparing 3D Genome Organization in Multiple Species Using Phylo-HMRF.

Recent whole-genome mapping approaches for the chromatin interactome have offered new insights into 3D genome organization. However, our knowledge of the evolutionary patterns of 3D genome in mammals remains limited. In particular, there are no existing phylogenetic-model-based methods to analyze chromatin interactions as continuous features. Here, we develop phylogenetic hidden Markov random field (Phylo-HMRF) to identify evolutionary patterns of 3D genome based on multi-species Hi-C data by jointly utilizing spatial constraints among genomic loci and continuous-trait evolutionary models. We used Phylo-HMRF to uncover cross-species 3D genome patterns based on Hi-C data from the same cell type in four primate species (human, chimpanzee, bonobo, and gorilla). The identified evolutionary patterns of 3D genome correlate with features of genome structure and function. This work provides a new framework to analyze multi-species continuous genomic features with spatial constraints and has the potential to help reveal the evolutionary principles of 3D genome organization.

Cost-effective scat-detection dogs: unleashing a powerful new tool for international mammalian conservation biology.

Recently, detection dogs have been utilized to collect fecal samples from cryptic and rare mammals. Despite the great promise of this technique for conservation biology, its broader application has been limited by the high cost (tens to hundreds of thousands of dollars) and logistical challenges of employing a scat-detection dog team while conducting international, collaborative research. Through an international collaboration of primatologists and the Chinese Ministry of Public Security, we trained and used a detection dog to find scat from three species of unhabituated, free-ranging primates, for less than $3,000. We collected 137 non-human primate fecal samples that we confirmed by sequencing taxonomically informative genetic markers. Our detection dog team had a 92% accuracy rate, significantly outperforming our human-only team. Our results demonstrate that detection dogs can locate fecal samples from unhabituated primates with variable diets, locomotion, and grouping patterns, despite challenging field conditions. We provide a model for in-country training, while also building local capacity for conservation and genetic monitoring. Unlike previous efforts, our approach will allow for the wide adoption of scat-detection dogs in international conservation biology.

Phylogenetic Analysis Supports a Link between DUF1220 Domain Number and Primate Brain Expansion.

The expansion of DUF1220 domain copy number during human evolution is a dramatic example of rapid and repeated domain duplication. Although patterns of expression, homology, and disease associations suggest a role in cortical development, this hypothesis has not been robustly tested using phylogenetic methods. Here, we estimate DUF1220 domain counts across 12 primate genomes using a nucleotide Hidden Markov Model. We then test a series of hypotheses designed to examine the potential evolutionary significance of DUF1220 copy number expansion. Our results suggest a robust association with brain size, and more specifically neocortex volume. In contradiction to previous hypotheses, we find a strong association with postnatal brain development but not with prenatal brain development. Our results provide further evidence of a conserved association between specific loci and brain size across primates, suggesting that human brain evolution may have occurred through a continuation of existing processes.

Neighboring Genes Show Correlated Evolution in Gene Expression.

When considering the evolution of a gene's expression profile, we commonly assume that this is unaffected by its genomic neighborhood. This is, however, in contrast to what we know about the lack of autonomy between neighboring genes in gene expression profiles in extant taxa. Indeed, in all eukaryotic genomes genes of similar expression-profile tend to cluster, reflecting chromatin level dynamics. Does it follow that if a gene increases expression in a particular lineage then the genomic neighbors will also increase in their expression or is gene expression evolution autonomous? To address this here we consider evolution of human gene expression since the human-chimp common ancestor, allowing for both variation in estimation of current expression level and error in Bayesian estimation of the ancestral state. We find that in all tissues and both sexes, the change in gene expression of a focal gene on average predicts the change in gene expression of neighbors. The effect is highly pronounced in the immediate vicinity (<100 kb) but extends much further. Sex-specific expression change is also genomically clustered. As genes increasing their expression in humans tend to avoid nuclear lamina domains and be enriched for the gene activator 5-hydroxymethylcytosine, we conclude that, most probably owing to chromatin level control of gene expression, a change in gene expression of one gene likely affects the expression evolution of neighbors, what we term expression piggybacking, an analog of hitchhiking.

The apparent enhancement of CpG transversions in primate lineage is a consequence of multiple replacements.

We claim that the apparently enhanced CpG transversions in the form CpG to CpC/GpG or to ApG/CpT are caused by the hypermutable CpG to CpA/TpG transition. The nucleotide replacement counts obtained from the human/chimpanzee/gorilla/orangutan sequence alignments representing the replacements due to the evolutionary species divergence and the results of 1000 genomes project that provide us with the differences due to the intraspecies diversification were analyzed to estimate the ratio of CpG versus non-CpG transversion probabilities. The trinucleotide replacement counts were extracted from the regions that are free of functional constraints. The CpG transversion probabilities based upon the genomic comparisons were found to exceed more than twice the non-CpG transversions. The diversity data emerging from 14 population groups were partitioned in five classes as a function of the parameter quantifying the spread of the polymorphic allele among the group of individuals. The results based upon the human polymorphism exhibit a trend where CpG over non-CpG transversion probability ratio is less and less exceeding unity as the values of the derived allele frequency (DAF) of snps are diminishing. A computer simulation of a simplified model indicates that the phenomenon of the apparent enhancement of CpG transversions can have its source in the interference of the entropic effects with the maximum likelihood methodologies.

Unsupervised genome-wide recognition of local relationship patterns.

BACKGROUND: Phenomena such as incomplete lineage sorting, horizontal gene transfer, gene duplication and subsequent sub- and neo-functionalisation can result in distinct local phylogenetic relationships that are discordant with species phylogeny. In order to assess the possible biological roles for these subdivisions, they must first be identified and characterised, preferably on a large scale and in an automated fashion. RESULTS: We developed Saguaro, a combination of a Hidden Markov Model (HMM) and a Self Organising Map (SOM), to characterise local phylogenetic relationships among aligned sequences using cacti, matrices of pair-wise distance measures. While the HMM determines the genomic boundaries from aligned sequences, the SOM hypothesises new cacti in an unsupervised and iterative fashion based on the regions that were modelled least well by existing cacti. After testing the software on simulated data, we demonstrate the utility of Saguaro by testing two different data sets: (i) 181 Dengue virus strains, and (ii) 5 primate genomes. Saguaro identifies regions under lineage-specific constraint for the first set, and genomic segments that we attribute to incomplete lineage sorting in the second dataset. Intriguingly for the primate data, Saguaro also classified an additional ~3% of the genome as most incompatible with the expected species phylogeny. A substantial fraction of these regions was found to overlap genes associated with both the innate and adaptive immune systems. CONCLUSIONS: Saguaro detects distinct cacti describing local phylogenetic relationships without requiring any a priori hypotheses. We have successfully demonstrated Saguaro's utility with two contrasting data sets, one containing many members with short sequences (Dengue viral strains: n = 181, genome size = 10,700 nt), and the other with few members but complex genomes (related primate species: n = 5, genome size = 3 Gb), suggesting that the software is applicable to a wide variety of experimental populations. Saguaro is written in C++, runs on the Linux operating system, and can be downloaded from http://saguarogw.sourceforge.net/.

Phylogenetic patterns of GC-biased gene conversion in placental mammals and the evolutionary dynamics of recombination landscapes.

GC-biased gene conversion (gBGC) is a major evolutionary force shaping genomic nucleotide landscapes, distorting the estimation of the strength of selection, and having potentially deleterious effects on genome-wide fitness. Yet, a global quantitative picture, at large evolutionary scale, of the relative strength of gBGC compared with selection and random drift is still lacking. Furthermore, owing to its dependence on the local recombination rate, gBGC results in modulations of the substitution patterns along genomes and across time which, if correctly interpreted, may yield quantitative insights into the long-term evolutionary dynamics of recombination landscapes. Deriving a model of the substitution process at putatively neutral nucleotide positions from population-genetics arguments, and accounting for among-lineage and among-gene effects, we propose a reconstruction of the variation in gBGC intensity at the scale of placental mammals, and of its scaling with body-size and karyotypic traits. Our results are compatible with a simple population genetics model relating gBGC to effective population size and recombination rate. In addition, among-gene variation and phylogenetic patterns of exon-specific levels of gBGC reveal the presence of rugged recombination landscapes, and suggest that short-lived recombination hot-spots are a general feature of placentals. Across placental mammals, variation in gBGC strength spans two orders of magnitude, at its lowest in apes, strongest in lagomorphs, microbats or tenrecs, and near or above the nearly neutral threshold in most other lineages. Combined with among-gene variation, such high levels of biased gene conversion are likely to significantly impact midly selected positions, and to represent a substantial mutation load. Altogether, our analysis suggests a more important role of gBGC in placental genome evolution, compared with what could have been anticipated from studies conducted in anthropoid primates.

The resurgence and genetic implications of New World primates in biomedical research.

There has been a recent resurgence of interest in New World monkeys within the biomedical research community, driven by both the sequencing of the common marmoset (Callithrix jacchus) genome and a growing demand for alternatives to Old World primates. New World monkeys offer attractive advantages over Old World species, including cheaper and simpler husbandry, while still maintaining a greater evolutionary proximity to humans compared with other animal models. Although numerous commonalities across primate species exist, there are also important genetic and reproductive differences that can and should play a critical role in selecting appropriate animal models. Common marmosets in particular have significantly reduced diversity at the major histocompatibility complex (MHC) loci and are born as hematopoietic chimeras. New World primates can make ideal translational models for research, but scientists must necessarily incorporate complete understandings of their genetic and phenotypic differences from humans and other model organisms.

Genome resource banking of biomedically important laboratory animals.

Genome resource banking is the systematic collection, storage, and redistribution of biomaterials in an organized, logistical, and secure manner. Genome cryobanks usually contain biomaterials and associated genomic information essential for progression of biomedicine, human health, and research. In that regard, appropriate genome cryobanks could provide essential biomaterials for both current and future research projects in the form of various cell types and tissues, including sperm, oocytes, embryos, embryonic or adult stem cells, induced pluripotent stem cells, and gonadal tissues. In addition to cryobanked germplasm, cryobanking of DNA, serum, blood products, and tissues from scientifically, economically, and ecologically important species has become a common practice. For revitalization of the whole organism, cryopreserved germplasm in conjunction with assisted reproductive technologies, offer a powerful approach for research model management, as well as assisting in animal production for agriculture, conservation, and human reproductive medicine. Recently, many developed and developing countries have allocated substantial resources to establish genome resources banks which are responsible for safeguarding scientifically, economically, and ecologically important wild type, mutant, and transgenic plants, fish, and local livestock breeds, as well as wildlife species. This review is dedicated to the memory of Dr. John K. Critser, who has made profound contributions to the science of cryobiology and establishment of genome research and resources centers for mice, rats, and swine. Emphasis will be given to application of genome resource banks to species with substantial contributions to the advancement of biomedicine and human health.

Codon substitution models based on residue similarity and their applications.

Codon models are now widely used to draw evolutionary inferences from alignments of homologous sequence data. Incorporating physicochemical properties of amino acids into codon models, two novel codon substitution models describing the evolution of protein-coding DNA sequences are presented based on the similarity scores of amino acids. To describe substitutions between codons a continue-time Markov process is used. Transition/transversion rate bias and nonsynonymous codon usage bias are allowed in the models. In our implementation, the parameters are estimated by maximum-likelihood (ML) method as in previous studies. Furthermore, instantaneous mutations involving more than one nucleotide position of a codon are considered in the second model. Then the two suggested models are applied to five real data sets. The analytic results indicate that the new codon models considering physicochemical properties of amino acids can provide a better fit to the data comparing with existing codon models, and then produce more reliable estimates of certain biologically important measures than existing methods.

Evaluating the phylogenetic position of Chinese tree shrew (Tupaia belangeri chinensis) based on complete mitochondrial genome: implication for using tree shrew as an alternative experimental animal to primates in biomedical research.

Tree shrew (Tupaia belangeri) is currently placed in Order Scandentia and has a wide distribution in Southeast Asia and Southwest China. Due to its unique characteristics, such as small body size, high brain-to-body mass ratio, short reproductive cycle and life span, and low-cost of maintenance, tree shrew has been proposed to be an alternative experimental animal to primates in biomedical research. However, there are some debates regarding the exact phylogenetic affinity of tree shrew to primates. In this study, we determined the mtDNA entire genomes of three Chinese tree shrews (T. belangeri chinensis) and one Malayan flying lemur (Galeopterus variegatus). Combined with the published data for species in Euarchonta, we intended to discern the phylogenetic relationship among representative species of Dermoptera, Scandentia and Primates. The mtDNA genomes of Chinese tree shrews and Malayan flying lemur shared similar gene organization and structure with those of other mammals. Phylogenetic analysis based on 12 concatenated mitochondrial protein-encoding genes revealed a closer relationship between species of Scandentia and Glires, whereas species of Dermoptera were clustered with Primates. This pattern was consistent with previously reported phylogeny based on mtDNA data, but differed from the one reconstructed on the basis of nuclear genes. Our result suggested that the matrilineal affinity of tree shrew to primates may not be as close as we had thought. The ongoing project for sequencing the entire genome of Chinese tree shrew will provide more information to clarify this important issue.

Evidence for a convergent slowdown in primate molecular rates and its implications for the timing of early primate evolution.

A long-standing problem in primate evolution is the discord between paleontological and molecular clock estimates for the time of crown primate origins: the earliest crown primate fossils are ~56 million y (Ma) old, whereas molecular estimates for the haplorhine-strepsirrhine split are often deep in the Late Cretaceous. One explanation for this phenomenon is that crown primates existed in the Cretaceous but that their fossil remains have not yet been found. Here we provide strong evidence that this discordance is better-explained by a convergent molecular rate slowdown in early primate evolution. We show that molecular rates in primates are strongly and inversely related to three life-history correlates: body size (BS), absolute endocranial volume (EV), and relative endocranial volume (REV). Critically, these traits can be reconstructed from fossils, allowing molecular rates to be predicted for extinct primates. To this end, we modeled the evolutionary history of BS, EV, and REV using data from both extinct and extant primates. We show that the primate last common ancestor had a very small BS, EV, and REV. There has been a subsequent convergent increase in BS, EV, and REV, indicating that there has also been a convergent molecular rate slowdown over primate evolution. We generated a unique timescale for primates by predicting molecular rates from the reconstructed phenotypic values for a large phylogeny of living and extinct primates. This analysis suggests that crown primates originated close to the K-Pg boundary and possibly in the Paleocene, largely reconciling the molecular and fossil timescales of primate evolution.

Dating primate divergences through an integrated analysis of palaeontological and molecular data.

Estimation of divergence times is usually done using either the fossil record or sequence data from modern species. We provide an integrated analysis of palaeontological and molecular data to give estimates of primate divergence times that utilize both sources of information. The number of preserved primate species discovered in the fossil record, along with their geological age distribution, is combined with the number of extant primate species to provide initial estimates of the primate and anthropoid divergence times. This is done by using a stochastic forwards-modeling approach where speciation and fossil preservation and discovery are simulated forward in time. We use the posterior distribution from the fossil analysis as a prior distribution on node ages in a molecular analysis. Sequence data from two genomic regions (CFTR on human chromosome 7 and the CYP7A1 region on chromosome 8) from 15 primate species are used with the birth-death model implemented in mcmctree in PAML to infer the posterior distribution of the ages of 14 nodes in the primate tree. We find that these age estimates are older than previously reported dates for all but one of these nodes. To perform the inference, a new approximate Bayesian computation (ABC) algorithm is introduced, where the structure of the model can be exploited in an ABC-within-Gibbs algorithm to provide a more efficient analysis.

Modelling the ancestral sequence distribution and model frequencies in context-dependent models for primate non-coding sequences.

BACKGROUND: Recent approaches for context-dependent evolutionary modelling assume that the evolution of a given site depends upon its ancestor and that ancestor's immediate flanking sites. Because such dependency pattern cannot be imposed on the root sequence, we consider the use of different orders of Markov chains to model dependence at the ancestral root sequence. Root distributions which are coupled to the context-dependent model across the underlying phylogenetic tree are deemed more realistic than decoupled Markov chains models, as the evolutionary process is responsible for shaping the composition of the ancestral root sequence. RESULTS: We find strong support, in terms of Bayes Factors, for using a second-order Markov chain at the ancestral root sequence along with a context-dependent model throughout the remainder of the phylogenetic tree in an ancestral repeats dataset, and for using a first-order Markov chain at the ancestral root sequence in a pseudogene dataset. Relaxing the assumption of a single context-independent set of independent model frequencies as presented in previous work, yields a further drastic increase in model fit. We show that the substitution rates associated with the CpG-methylation-deamination process can be modelled through context-dependent model frequencies and that their accuracy depends on the (order of the) Markov chain imposed at the ancestral root sequence. In addition, we provide evidence that this approach (which assumes that root distribution and evolutionary model are decoupled) outperforms an approach inspired by the work of Arndt et al., where the root distribution is coupled to the evolutionary model. We show that the continuous-time approximation of Hwang and Green has stronger support in terms of Bayes Factors, but the parameter estimates show minimal differences. CONCLUSIONS: We show that the combination of a dependency scheme at the ancestral root sequence and a context-dependent evolutionary model across the remainder of the tree allows for accurate estimation of the model's parameters. The different assumptions tested in this manuscript clearly show that designing accurate context-dependent models is a complex process, with many different assumptions that require validation. Further, these assumptions are shown to change across different datasets, making the search for an adequate model for a given dataset quite challenging.

Divergence, demography and gene loss along the human lineage.

Genomic DNA sequences are an irreplaceable source for reconstructing the vanished past of living organisms. Based on updated sequence data, this paper summarizes our studies on species divergence time, ancient population size and functional loss of genes in the primate lineage leading to modern humans (Homo sapiens sapiens). The inter- and intraspecific comparisons of DNA sequences suggest that the human lineage experienced a rather severe bottleneck in the Middle Pleistocene, throughout which period the subdivided African population played a predominant role in shaping the genetic architecture of modern humans. Also, published and newly identified human-specific pseudogenes (HSPs) are enumerated in order to infer their significance for human evolution. Of the 121 candidate genes obtained, authentic HSPs turn out to comprise only 25 olfactory receptor genes, four T cell receptor genes and nine other genes. The fixation of HSPs has been too rare over the past 6-7 Myr to account for species differences between humans and chimpanzees.

Genome assembly quality: assessment and improvement using the neutral indel model.

We describe a statistical and comparative-genomic approach for quantifying error rates of genome sequence assemblies. The method exploits not substitutions but the pattern of insertions and deletions (indels) in genome-scale alignments for closely related species. Using two- or three-way alignments, the approach estimates the amount of aligned sequence containing clusters of nucleotides that were wrongly inserted or deleted during sequencing or assembly. Thus, the method is well-suited to assessing fine-scale sequence quality within single assemblies, between different assemblies of a single set of reads, and between genome assemblies for different species. When applying this approach to four primate genome assemblies, we found that average gap error rates per base varied considerably, by up to sixfold. As expected, bacterial artificial chromosome (BAC) sequences contained lower, but still substantial, predicted numbers of errors, arguing for caution in regarding BACs as the epitome of genome fidelity. We then mapped short reads, at approximately 10-fold statistical coverage, from a Bornean orangutan onto the Sumatran orangutan genome assembly originally constructed from capillary reads. This resulted in a reduced gap error rate and a separation of error-prone from high-fidelity sequence. Over 5000 predicted indel errors in protein-coding sequence were corrected in a hybrid assembly. Our approach contributes a new fine-scale quality metric for assemblies that should facilitate development of improved genome sequencing and assembly strategies.

Nepotistic cooperation in non-human primate groups.

Darwin was struck by the many similarities between humans and other primates and believed that these similarities were the product of common ancestry. He would be even more impressed by the similarities if he had known what we have learned about primates over the last 50 years. Genetic kinship has emerged as the primary organizing force in the evolution of primate social organization and the patterning of social behaviour in non-human primate groups. There are pronounced nepotistic biases across the primate order, from tiny grey mouse lemurs (Microcebus murinus) that forage alone at night but cluster with relatives to sleep during the day, to cooperatively breeding marmosets that rely on closely related helpers to rear their young, rhesus macaque (Macaca mulatta) females who acquire their mother's rank and form strict matrilineal dominance hierarchies, male howler monkeys that help their sons maintain access to groups of females and male chimpanzees (Pan troglodytes) that form lasting relationships with their brothers. As more evidence of nepotism has accumulated, important questions about the evolutionary processes underlying these kin biases have been raised. Although kin selection predicts that altruism will be biased in favour of relatives, it is difficult to assess whether primates actually conform to predictions derived from Hamilton's rule: br > c. In addition, other mechanisms, including contingent reciprocity and mutualism, could contribute to the nepotistic biases observed in non-human primate groups. There are good reasons to suspect that these processes may complement the effects of kin selection and amplify the extent of nepotistic biases in behaviour.

Widespread genomic signatures of natural selection in hominid evolution.

Selection acting on genomic functional elements can be detected by its indirect effects on population diversity at linked neutral sites. To illuminate the selective forces that shaped hominid evolution, we analyzed the genomic distributions of human polymorphisms and sequence differences among five primate species relative to the locations of conserved sequence features. Neutral sequence diversity in human and ancestral hominid populations is substantially reduced near such features, resulting in a surprisingly large genome average diversity reduction due to selection of 19-26% on the autosomes and 12-40% on the X chromosome. The overall trends are broadly consistent with "background selection" or hitchhiking in ancestral populations acting to remove deleterious variants. Average selection is much stronger on exonic (both protein-coding and untranslated) conserved features than non-exonic features. Long term selection, rather than complex speciation scenarios, explains the large intragenomic variation in human/chimpanzee divergence. Our analyses reveal a dominant role for selection in shaping genomic diversity and divergence patterns, clarify hominid evolution, and provide a baseline for investigating specific selective events.

Patterns of co-speciation and host switching in primate malaria parasites.

BACKGROUND: The evolutionary history of many parasites is dependent on the evolution of their hosts, leading to an association between host and parasite phylogenies. However, frequent host switches across broad phylogenetic distances may weaken this close evolutionary link, especially when vectors are involved in parasites transmission, as is the case for malaria pathogens. Several studies suggested that the evolution of the primate-infective malaria lineages may be constrained by the phylogenetic relationships of their hosts, and that lateral switches between distantly related hosts may have been occurred. However, no systematic analysis has been quantified the degree of phylogenetic association between primates and their malaria parasites. METHODS: Here phylogenetic approaches have been used to discriminate statistically between events due to co-divergence, duplication, extinction and host switches that can potentially cause historical association between Plasmodium parasites and their primate hosts. A Bayesian reconstruction of parasite phylogeny based on genetic information for six genes served as basis for the analyses, which could account for uncertainties about the evolutionary hypotheses of malaria parasites. RESULTS: Related lineages of primate-infective Plasmodium tend to infect hosts within the same taxonomic family. Different analyses testing for congruence between host and parasite phylogenies unanimously revealed a significant association between the corresponding evolutionary trees. The most important factor that resulted in this association was host switching, but depending on the parasite phylogeny considered, co-speciation and duplication may have also played some additional role. Sorting seemed to be a relatively infrequent event, and can occur only under extreme co-evolutionary scenarios. The concordance between host and parasite phylogenies is heterogeneous: while the evolution of some malaria pathogens is strongly dependent on the phylogenetic history of their primate hosts, the congruent evolution is less emphasized for other parasite lineages (e.g. for human malaria parasites). Estimation of ancestral states of host use along the phylogenetic tree of parasites revealed that lateral transfers across distantly related hosts were likely to occur in several cases. Parasites cannot infect all available hosts, and they should preferentially infect hosts that provide a similar environment for reproduction. Marginally significant evidence suggested that there might be a consistent variation within host ranges in terms of physiology. CONCLUSION: The evolution of primate malarias is constrained by the phylogenetic associations of their hosts. Some parasites can preserve a great flexibility to infect hosts across a large phylogenetic distance, thus host switching can be an important factor in mediating host ranges observed in nature. Due to this inherent flexibility and the potential exposure to various vectors, the emergence of new malaria disease in primates including humans cannot be predicted from the phylogeny of parasites.