PubMed Health⌕ Search

Biomedical subjects

D Posada

Publications and source records attributed to D Posada.

14 recordsLinked to original sources

Evaluation of methods for detecting recombination from DNA sequences: computer simulations.

Recombination is a key evolutionary process that shapes the architecture of genomes and the genetic structure of populations. Although many statistical methods are available for the detection of recombination from DNA sequences, their absolute and relative performance is still unknown. Here we evaluated the performance of 14 different recombination detection algorithms. We used the coalescent with recombination to simulate DNA sequences with different levels of recombination, genetic diversity, and rate variation among sites. Recombination detection methods were applied to these data sets, and whether they detected or not recombination was recorded. Different recombination methods showed distinct performance depending on the amount of recombination, genetic diversity, and rate variation among sites. The model of nucleotide substitution under which the data were generated did not seem to have a significant effect. Most methods increase power with more sequence divergence. In general, recombination detection methods seem to capture the presence of recombination, but they are not very powerful. Methods that use substitution patterns or incompatibility among sites were more powerful than methods based on phylogenetic incongruence. Most methods do not seem to infer more false positives than expected by chance. Especially depending on the amount of diversity in the data, different methods could be used to attain maximum power while minimizing false positives. Results shown here will provide some guidance in the selection of the most appropriate method/s for the analysis of the particular data at hand.

Computer Simulation↗

Intraspecific gene genealogies: trees grafting into networks.

Intraspecific gene evolution cannot always be represented by a bifurcating tree. Rather, population genealogies are often multifurcated, descendant genes coexist with persistent ancestors and recombination events produce reticulate relationships. Whereas traditional phylogenetic methods assume bifurcating trees, several networking approaches have recently been developed to estimate intraspecific genealogies that take into account these population-level phenomena.

Journal Article↗

The effect of branch length variation on the selection of models of molecular evolution.

Models of sequence evolution play an important role in molecular evolutionary studies. The use of inappropriate models of evolution may bias the results of the analysis and lead to erroneous conclusions. Several procedures for selecting the best-fit model of evolution for the data at hand have been proposed, like the likelihood ratio test (LRT) and the Akaike (AIC) and Bayesian (BIC) information criteria. The relative performance of these model-selecting algorithms has not yet been studied under a range of different model trees. In this study, the influence of branch length variation upon model selection is characterized. This is done by simulating sequence alignments under a known model of nucleotide substitution, and recording how often this true model is recovered by different model-fitting strategies. Results of this study agree with previous simulations and suggest that model selection is reasonably accurate. However, different model selection methods showed distinct levels of accuracy. Some LRT approaches showed better performance than the AIC or BIC information criteria. Within the LRTs, model selection is affected by the complexity of the initial model selected for the comparisons, and only slightly by the order in which different parameters are added to the model. A specific hierarchy of LRTs, which starts from a simple model of evolution, performed overall better than other possible LRT hierarchies, or than the AIC or BIC.

Bayes Theorem↗

Selecting models of nucleotide substitution: an application to human immunodeficiency virus 1 (HIV-1).

The blind use of models of nucleotide substitution in evolutionary analyses is a common practice in the viral community. Typically, a simple model of evolution like the Kimura two-parameter model is used for estimating genetic distances and phylogenies, either because other authors have used it or because it is the default in various phylogenetic packages. Using two statistical approaches to model fitting, hierarchical likelihood ratio tests and the Akaike information criterion, we show that different viral data sets are better explained by different models of evolution. We demonstrate our results with the analysis of HIV-1 sequences from a hierarchy of samples; sequences within individuals, individuals within subtypes, and subtypes within groups. We also examine results for three different gene regions: gag, pol, and env. The Kimura two-parameter model was not selected as the best-fit model for any of these data sets, despite its widespread use in phylogenetic analyses of HIV-1 sequences. Furthermore, the model complexity increased with increasing sequence divergence. Finally, the molecular-clock hypothesis was rejected in most of the data sets analyzed, throwing into question clock-based estimates of divergence times for HIV-1. The importance of models in evolutionary analyses and their repercussions on the derived conclusions are discussed.

Databases, Factual↗

Molecular systematics of European Hyalodaphnia: the role of contemporary hybridization in ancient species.

We examined phylogenetic relationships among Daphnia using mitochondrial DNA (mtDNA) sequences from the small subunit ribosomal RNA (12S), cytochrome c oxidase subunit I and nuclear DNA sequences from the first and second internal transcribed spacer representing 1612 base positions. Phylogenetic analyses using several species of the three main Daphnia subgenera, Ctenodaphnia, Hyalodaphnia and Daphnia, revealed that the Hyalodaphnia are a monophyletic sister group of the Daphnia. Most Hyalodaphnia species occur on one continent, whereas only three are found in North America and Europe. Endemicity of species is associated with variation in thermal tolerance and habitat differentiation. Although many species of the Hyalodaphnia are known to hybridize in nature, mtDNA divergence is relatively high ca. 9%) compared to other hybridizing arthropods (ca. 3%). Reproductive isolation in Daphnia seems to evolve significantly slower than genetic isolation. We related these findings to what is known about the ecology and genetics of Daphnia in order to better understand the evolutionary diversification of lineages. The relationship of these data to phylogenetic patterns is discussed in the context of speciation processes in Daphnia.

Animals↗

Population genetics of the porB gene of Neisseria gonorrhoeae: different dynamics in different homology groups.

The porB locus codes for the major outer membrane protein of Neisseria gonorrhoeae. Alleles of this locus have been assigned to two homology groups based on close sequence and immunological relationships and are designated as either PIA or PIB. Several population parameters were estimated and compared among these two groups using a data set of 22 PIA sequences and 91 PIB sequences obtained from diverse geographic localities and from time periods spanning approximately 50 years. Recombination appears to be extensive in the porB gene. While the recombination rates are similar for the PIA and PIB sequences, the relative contribution of recombination to genetic diversity is higher for the PIA sequences. Alleles belonging to the PIB group show greater genetic diversity than do those in the PIA group. Although phylogenetic analysis did not reveal temporal or geographic clustering of sequences, estimates of gene flow and the fixation index suggested that PIB sequences exhibit population substructure based on geographic locality. Selection acts in these homology groups in a different way. While positive Darwinian selection is the dominant force driving the evolution of the PIA sequences, purifying selection operates also on the PIB sequences. These differences may be attributable to the greater propensity of PIA strains, as compared with PIB strains, to cause disseminated gonococcal infection, which would expose the former to intense selection pressure from the host immune system. The molecular evolution of Neisseria gonorrhoeae seems to be driven by the simultaneous action of selection and recombination, but under different rates and selection pressures for the PIA and PIB homology groups.

Alleles↗

Mitochondrial DNA phylogeography and population history of the grey wolf canis lupus

The grey wolf (Canis lupus) and coyote (C. latrans) are highly mobile carnivores that disperse over great distances in search of territories and mates. Previous genetic studies have shown little geographical structure in either species. However, population genetic structure is also influenced by past isolation events and population fluctuations during glacial periods. In this study, control region sequence data from a worldwide sample of grey wolves and a more limited sample of coyotes were analysed. The results suggest that fluctuating population sizes during the late Pleistocene have left a genetic signature on levels of variation in both species. Genealogical measures of nucleotide diversity suggest that historical population sizes were much larger in both species and grey wolves were more numerous than coyotes. Currently, about 300 000 wolves and 7 million coyotes exist. In grey wolves, genetic diversity is greater than that predicted from census population size, reflecting recent historical population declines. By contrast, nucleotide diversity in coyotes is smaller than that predicted by census population size, reflecting a recent population expansion following the extirpation of wolves from much of North America. Both species show little partitioning of haplotypes on continental or regional scales. However, a statistical parsimony analysis indicates local genetic structure that suggests recent restricted gene flow.

Journal Article↗

MODELTEST: testing the model of DNA substitution.

SUMMARY: The program MODELTEST uses log likelihood scores to establish the model of DNA evolution that best fits the data. AVAILABILITY: The MODELTEST package, including the source code and some documentation is available at http://bioag.byu. edu/zoology/crandall_lab/modeltest.html.

Algorithms↗