The Influence of Taxon Sampling and Tree Shape on Molecular Dating: An Empirical Example from Mammalian Mitochondrial Genomes

Over the last decade, molecular dating methods have been among the most studied subjects in statistical phylogenetics. Although the evolutionary modelling of substitution rates and the handling of calibration information are the primary focus of species divergence time research, parameters that influence topological estimation, such as taxon sampling and tree shape, also have the potential to influence evolutionary age estimates. However, the impact of topological parameters on chronological estimates is rarely considered. In this study, we use mitochondrial genomes to evaluate the influence of tree shape and taxon sampling on the divergence times of selected nodes of the mammalian tree. Our results show that taxon sampling affects divergence time estimates; the credibility intervals for age estimates decrease as taxonomic sampling increases (i.e., estimates become more precise). The influence of taxonomic sampling was not observed on nodes that lay deep in the mammalian phylogeny, although the means of the posterior distributions tend to converge with increased taxon sampling, an effect that is independent of the location of the node. In the majority of cases, the effect of tree shape was negligible

Impact of the Partitioning Scheme on Divergence Times Inferred from Mammalian Genomic Data Sets

Data partitioning has long been regarded as an important parameter for phylogenetic inference. The division of heterogeneous multigene data sets into partitions with similar substitution patterns is known to increase the performance of probabilistic phylogenetic methods. However, the effect of the partitioning scheme on divergence time estimates has generally been ignored. To investigate the impact of data partitioning on the estimation of divergence times, we have constructed two genomic data sets. The first one with 15 nuclear genes comprising 50,928 bp were selected from the OrthoMam database; the second set was composed of complete mitochondrial genomes. We studied two partitioning schemes: concatenated supermatrices and partitioned gene analysis. We have also measured the impact of taxonomic sampling on the estimates. After drawing divergence time inferences using the uncorrelated relaxed clock in BEAST, we have compared the age estimates between the partitioning schemes. Our results show that, in general, both schemes resulted in similar chronological estimates, however the concatenated data sets were more efficient than the partitioned ones in attaining suitable effective sample sizes

The Use of Bioinformatics for Studying HIV Evolutionary and Epidemiological History in South America

The South American human immunodeficiency virus type 1 (HIV-1) epidemic is driven by several subtypes (B, C, and F1) and circulating and unique recombinant forms derived from those subtypes. Those variants are heterogeneously distributed around the continent in a country-specific manner. Despite some inconsistencies mainly derived from sampling biases and analytical constrains, most of studies carried out in the area agreed in pointing out specificities in the evolutionary dynamics of the circulating HIV-1 lineages. In this paper, we covered the theoretical basis, and the application of bioinformatics methods to reconstruct the HIV spatial-temporal dynamics, unveiling relevant information to understand the origin, geographical dissemination and the current molecular scenario of the HIV epidemic in the continent, particularly in the countries of Southern Cone

Phylogenetic nomenclature and evolution of mannose-binding lectin (MBL2) haplotypes

<p>Abstract</p> <p>Background</p> <p>Polymorphisms of the mannose-binding lectin gene (<it>MBL2</it>) affect the concentration and functional efficiency of the protein. We recently used haplotype-specific sequencing to identify 23 <it>MBL2 </it>haplotypes, associated with enhanced susceptibility to several diseases.</p> <p>Results</p> <p>In this work, we applied the same method in 288 and 470 chromosomes from Gabonese and European adults, respectively, and found three new haplotypes in the last group. We propose a phylogenetic nomenclature to standardize <it>MBL2 </it>studies and found two major phylogenetic branches due to six strongly linked polymorphisms associated with high MBL production. They presented high Fst values and were imbedded in regions with high nucleotide diversity and significant Tajima's D values. Compared to others using small sample sizes and unphased genotypic data, we found differences in haplotyping, frequency estimation, Fu and Li's D* and Fst results.</p> <p>Conclusion</p> <p>Using extensive testing for selective neutrality, we confirmed that stochastic evolutionary factors have had a major role in shaping this polymorphic gene worldwide.</p

Evolution of the B3 DNA Binding Superfamily: New Insights into REM Family Gene Diversification

Background: The B3 DNA binding domain includes five families: auxin response factor (ARF), abscisic acid-insensitive3 (ABI3), high level expression of sugar inducible (HSI), related to ABI3/VP1 (RAV) and reproductive meristem (REM). The release of the complete genomes of the angiosperm eudicots Arabidopsis thaliana and Populus trichocarpa, the monocot Orysa sativa, the bryophyte Physcomitrella patens,the green algae Chlamydomonas reinhardtii and Volvox carteri and the red algae Cyanidioschyzon melorae provided an exceptional opportunity to study the evolution of this superfamily. Methodology: In order to better understand the origin and the diversification of B3 domains in plants, we combined comparative phylogenetic analysis with exon/intron structure and duplication events. In addition, we investigated the conservation and divergence of the B3 domain during the origin and evolution of each family. Conclusions: Our data indicate that showed that the B3 containing genes have undergone extensive duplication events, and that the REM family B3 domain has a highly diverged DNA binding. Our results also indicate that the founding member of the B3 gene family is likely to be similar to the ABI3/HSI genes found in C. reinhardtii and V. carteri. Among the B3 families, ABI3, HSI, RAV and ARF are most structurally conserved, whereas the REM family has experienced a rapid divergence. Thes

The Role of bZIP Transcription Factors in Green Plant Evolution: Adaptive Features Emerging from Four Founder Genes

BACKGROUND: Transcription factors of the basic leucine zipper (bZIP) family control important processes in all eukaryotes. In plants, bZIPs are regulators of many central developmental and physiological processes including photomorphogenesis, leaf and seed formation, energy homeostasis, and abiotic and biotic stress responses. Here we performed a comprehensive phylogenetic analysis of bZIP genes from algae, mosses, ferns, gymnosperms and angiosperms. METHODOLOGY/PRINCIPAL FINDINGS: We identified 13 groups of bZIP homologues in angiosperms, three more than known before, that represent 34 Possible Groups of Orthologues (PoGOs). The 34 PoGOs may correspond to the complete set of ancestral angiosperm bZIP genes that participated in the diversification of flowering plants. Homologous genes dedicated to seed-related processes and ABA-mediated stress responses originated in the common ancestor of seed plants, and three groups of homologues emerged in the angiosperm lineage, of which one group plays a role in optimizing the use of energy. CONCLUSIONS/SIGNIFICANCE: Our data suggest that the ancestor of green plants possessed four bZIP genes functionally involved in oxidative stress and unfolded protein responses that are bZIP-mediated processes in all eukaryotes, but also in light-dependent regulations. The four founder genes amplified and diverged significantly, generating traits that benefited the colonization of new environments

Assignment of Calibration Information to Deeper Phylogenetic Nodes is More Effective in Obtaining Precise and Accurate Divergence Time Estimates

Divergence time estimation has become an essential tool for understanding macroevolutionary events. Molecular dating aims to obtain reliable inferences, which, within a statistical framework, means jointly increasing the accuracy and precision of estimates. Bayesian dating methods exhibit the propriety of a linear relationship between uncertainty and estimated divergence dates. This relationship occurs even if the number of sites approaches infinity and places a limit on the maximum precision of node ages. However, how the placement of calibration information may affect the precision of divergence time estimates remains an open question. In this study, relying on simulated and empirical data, we investigated how the location of calibration within a phylogeny affects the accuracy and precision of time estimates. We found that calibration priors set at median and deep phylogenetic nodes were associated with higher precision values compared to analyses involving calibration at the shallowest node. The results were independent of the tree symmetry. An empirical mammalian dataset produced results that were consistent with those generated by the simulated sequences. Assigning time information to the deeper nodes of a tree is crucial to guarantee the accuracy and precision of divergence times. This finding highlights the importance of the appropriate choice of outgroups in molecular dating

Arrival and diversification of mabuyine skinks (Squamata: Scincidae) in the Neotropics based on a fossil-calibrated timetree

Background The evolution of South American Mabuyinae skinks holds significant biogeographic interest because its sister lineage is distributed across the African continent and adjacent islands. Moreover, at least one insular species, Trachylepis atlantica, has independently reached the New World through transoceanic dispersal. To clarify the evolutionary history of both Neotropical lineages, this study aimed to infer an updated timescale using the largest species and gene sampling dataset ever assembled for this group. By extending the analysis to the Scincidae family, we could employ fossil information to estimate mabuyinae divergence times and carried out a formal statistical biogeography analysis. To unveil macroevolutionary patterns, we also inferred diversification rates for this lineage and evaluated whether the colonization of South American continent significantly altered the mode of Mabuyinae evolution. Methods A time-calibrated phylogeny was inferred under the Bayesian framework employing fossil information. This timetree was used to (i) evaluate the historical biogeography of mabuiyines using the statistical approach implemented in BioGeoBEARS; (ii) estimate macroevolutionary diversification rates of the South American Mabuyinae lineages and the patterns of evolution of selected traits, namely, the mode of reproduction, body mass and snout–vent length; (iii) test the hypothesis of differential macroevolutionary patterns in South American lineages in BAMM and GeoSSE; and (iv) re-evaluate the ancestral state of the mode of reproduction of mabuyines. Results Our results corroborated the hypothesis that the occupation of the South American continent by Mabuyinae consisted of two independent dispersion events that occurred between the Oligocene and the Miocene. We found significant differences in speciation rates between the New World and the remaining Mabuyinae clades only in GeoSSE. The influence of phenotypic traits on diversification rates was not supported by any method. Ancestral state reconstruction suggested that the ancestor of South American mabuyine was likely viviparous. Discussion Our analyses further corroborated the existence of a transoceanic connection between Africa and South America in the Eocene/Oligocene period (Atlantogea). Following colonization of the isolated South America and subsequent dispersal through the continent by the ancestral mabuyine stock, we detected no difference in macroevolutionary regimes of New World clades. This finding argued against the ecological opportunity model as an explanation for the diversity of living mabuyines