Performance of genomic data sets on the estimation of the divergence time of New World and Old World anthropoids.

The origin of New World anthropoids has received renewed attention since the advent of molecular dating methods that relax the assumption of a strict molecular clock. However, the studies conducted to date have estimated the time of the separation of New World and Old World anthropoids at values as different as 70 and 22 Ma. With the aim of investigating the source of the discrepancies in the inferred ages, we have compared the performance of mitochondrial and nuclear markers in two pairs of datasets. We show that in the larger genomic samples, the dates of the separation of New and Old World anthropoids estimated from nuclear and mitochondrial data are significantly different. The precision of the estimates demonstrated that both markers rendered significantly different estimates. However, parametric estimates from the large nuclear dataset were highly cross-correlated. Cross-correlation of absolute divergence times and evolutionary rates was as great as -96%. Consequently, the age estimates from the large nuclear data were not reproducible, because Markov chains were unable to reach the same parametric values independently, even with the adoption of additional information from calibration priors. Thus, because branch length decomposition was not achieved, a comparison of the genomic age estimates from nuclear and mitochondrial datasets was statistically impractical. We demonstrate the importance of examining the output of Markov chain Monte Carlo analyses for correlation between rate and time in studies that use phylogenomic datasets to examine the chronological scales of primate evolution.

The precision of the hominid timescale estimated by relaxed clock methods.

The chronological scenario of the evolution of hominoid primates has been thoroughly investigated since the advent of the molecular clock hypothesis. With the availability of genomic sequences for all hominid genera and other anthropoids, we may have reached the point at which the information from sequence data alone will not provide further evidence for the inference of the hominid evolution timescale. To verify this conjecture, we have compiled a genomic data set for all of the anthropoid genera. Our estimate places the Homo/Pan divergence at approximately 7.4 Ma, the Gorilla lineage divergence at approximately 9.7 Ma, the basal Hominidae divergence at 18.1 Ma and the basal Hominoidea divergence at 20.6 Ma. By inferring the theoretical limit distribution of posterior densities under a Bayesian framework, we show that it is unlikely that lengthier alignments or the availability of new genomic sequences will provide additional information to reduce the uncertainty associated with the divergence time estimates of the four hominid genera. A reduction of this uncertainty will be achieved only by the inclusion of more informative calibration priors.

Combining fossil and molecular data to date the diversification of New World Primates.

Recent methodological advances in molecular dating associated with the growing availability of sequence data have prompted the study of the evolution of New World Anthropoidea in recent years. Motivated by questions regarding historical biogeography or the mode of evolution, these works aimed to obtain a clearer scenario of Platyrrhini origins and diversification. Although some consensus was found, disputed issues, especially those relating to the evolutionary affinities of fossil taxa, remain. The use of fossil taxa for divergence time analysis is traditionally restricted to the provision of calibration priors. However, new analytical approaches have been developed that incorporate fossils as terminals and, thus, directly assign ages to the fossil tips. In this study, we conducted a combined analysis of molecular and morphological data, including fossils, to derive the timescale of New World anthropoids. Differently from previous studies that conducted total-evidence analysis of molecules and morphology, our approach investigated the morphological clock alone. Our results corroborate the hypothesis that living platyrrhines diversified in the last 20 Ma and that Miocene Patagonian fossils compose an independent evolutionary radiation that diversified in the late Oligocene. When compared to the node ages inferred from the molecular timescale, the inclusion of fossils augmented the precision of the estimates for nodes constrained by the fossil tips. We show that morphological data can be analysed using the same methodological framework applied in relaxed molecular clock studies.

The evolution of South American endemic canids: a history of rapid diversification and morphological parallelism.

The origin of endemic South American canid fauna has been traditionally linked with the rise of the Isthmus of Panama, suggesting that diversification of the dog fauna on this continent occurred very rapidly. Nevertheless, despite its obvious biogeographic appeal, the tempo of Canid evolution in South America has never been studied thoroughly. This issue can be suitably tackled with the inference of a molecular timescale. In this study, using a relaxed molecular clock method, we estimated that the most recent common ancestor of South American canids lived around 4 Ma, whereas all other splits within the clade occurred after the rise of the Panamanian land bridge. We suggest that the early diversification of the ancestors of the two main lineages of South American canids may have occurred in North America, before the Great American Interchange. Moreover, a concatenated morphological and molecular analysis put some extinct canid species well within the South American radiation, and shows that the dental adaptations to hypercarnivory evolved only once in the South American clade.

A molecular study on the evolution of a subtype B variant frequently found in Brazil.

In spite of the remarkable diversity of HIV-1 env genes, several amino acids are extremely conserved, probably due to functional constraints. One example is the proline found at the second position of the GPGR motif. Several viruses, however, bear substitutions at this site, for instance, GWGR subtype B variant. GWGR viruses are described in Brazil since the beginning of the epidemics, but the extent of their dispersion or the geographical origin of the variant remains unknown. In the present study, phylogenetic trees were constructed in order to study the origin and spread of this variant. All GWGR sequences as well as a subset of subtype B sequences available were included in the analyses. Analyses of differential selection were also performed on GWGR and non-GWGR sequences in order to unveil evolutionary novelties due to the action of positive selection. Although the GWGR variant was found at least in 23 countries, its expansion probably has a single origin, and Brazil is the epicenter.

Evaluating DNA evidence in a genetically complex population.

In forensic genetics, the likelihood ratio (LR), measuring the value of DNA profile evidence, is computed from a database of allele frequencies. Here, we address the choice of database and adjustments for population structure and sample size in the context of Brazil. The Brazilian population underwent a complex process of colonization, migration and mating, which created an admixed genetic composition that makes it difficult to obtain an appropriate database for a given case. National databases are now available, as well as databases for many Brazilian states. However, those databases are not statistically random samples, and state boundaries may not accurately reflect the sub-structuring of genetic diversity. We compared the LR calculated using the relevant state-specific database with the statistics calculated when a national database and when international databases were used. We evaluated two methods of adjustment for population structure, due to Wright [13] and Balding and Nichols [14]. We also considered two adjustments for database sample size: the Balding size bias correction [15] and a minimum allele frequency [16]. Our results show that the use of a national database with the Balding and Nichols adjustment and θ = 0.002 generated lower LR values than did the state-specific database in more than 50% of the profiles simulated using the state-based allele frequencies, while θ = 0.01 produced lower LRs for more than 90% of the profiles. We conclude that the utilization of a national database for Brazilian cases can be justified in association with the appropriate adjustment for population structure.

An empirical examination of the standard errors of maximum likelihood phylogenetic parameters under the molecular clock via bootstrapping

The molecular clock theory has greatly enlightened our understanding of macroevolutionary events. Maximum likelihood (ML) estimation of divergence times involves the adoption of fixed calibration points, and the confidence intervals associated with the estimates are generally very narrow. The credibility intervals are inferred assuming that the estimates are normally distributed, which may not be the case. Moreover, calculation of standard errors is usually carried out by the curvature method and is complicated by the difficulty in approximating second derivatives of the likelihood function. In this study, a standard primate phylogeny was used to examine the standard errors of ML estimates via the bootstrap method. Confidence intervals were also assessed from the posterior distribution of divergence times inferred via Bayesian Markov Chain Monte Carlo. For the primate topology under evaluation, no significant differences were found between the bootstrap and the curvature methods. Also, Bayesian confidence intervals were always wider than those obtained by ML.

Analysis of differential selective forces acting on the coat protein (P3) of the plant virus family Luteoviridae

The coat protein (CP) of the family Luteoviridae is directly associated with the success of infection. It participates in various steps of the virus life cycle, such as virion assembly, stability, systemic infection, and transmission. Despite its importance, extensive studies on the molecular evolution of this protein are lacking. In the present study, we investigate the action of differential selective forces on the CP coding region using maximum likelihood methods. We found that the protein is subjected to heterogeneous selective pressures and some sites may be evolving near neutrality. Based on the proposed 3-D model of the CP S-domain, we showed that nearly neutral sites are predominantly located in the region of the protein that faces the interior of the capsid, in close contact with the viral RNA, while highly conserved sites are mainly part of beta-strands, in the protein's major framework.