New Methods to Calculate Concordance Factors for Phylogenomic Datasets.
Ontology highlight
ABSTRACT:
SUBMITTER: Minh BQ
PROVIDER: S-EPMC7475031 | biostudies-literature | 2020 Sep
REPOSITORIES: biostudies-literature
Publications
New Methods to Calculate Concordance Factors for Phylogenomic Datasets.
Minh Bui Quang BQ Hahn Matthew W MW Lanfear Robert R
Molecular biology and evolution 20200901 9
We implement two measures for quantifying genealogical concordance in phylogenomic data sets: the gene concordance factor (gCF) and the novel site concordance factor (sCF). For every branch of a reference tree, gCF is defined as the percentage of "decisive" gene trees containing that branch. This measure is already in wide usage, but here we introduce a package that calculates it while accounting for variable taxon coverage among gene trees. sCF is a new measure defined as the percentage of deci ...[more]
Similar Datasets
Project description:BackgroundThe use of transcriptomic and genomic datasets for phylogenetic reconstruction has become increasingly common as researchers attempt to resolve recalcitrant nodes with increasing amounts of data. The large size and complexity of these datasets introduce significant phylogenetic noise and conflict into subsequent analyses. The sources of conflict may include hybridization, incomplete lineage sorting, or horizontal gene transfer, and may vary across the phylogeny. For phylogenetic analysis, this noise and conflict has been accommodated in one of several ways: by binning gene regions into subsets to isolate consistent phylogenetic signal; by using gene-tree methods for reconstruction, where conflict is presumed to be explained by incomplete lineage sorting (ILS); or through concatenation, where noise is presumed to be the dominant source of conflict. The results provided herein emphasize that analysis of individual homologous gene regions can greatly improve our understanding of the underlying conflict within these datasets.ResultsHere we examined two published transcriptomic datasets, the angiosperm group Caryophyllales and the aculeate Hymenoptera, for the presence of conflict, concordance, and gene duplications in individual homologs across the phylogeny. We found significant conflict throughout the phylogeny in both datasets and in particular along the backbone. While some nodes in each phylogeny showed patterns of conflict similar to what might be expected with ILS alone, the backbone nodes also exhibited low levels of phylogenetic signal. In addition, certain nodes, especially in the Caryophyllales, had highly elevated levels of strongly supported conflict that cannot be explained by ILS alone.ConclusionThis study demonstrates that phylogenetic signal is highly variable in phylogenomic data sampled across related species and poses challenges when conducting species tree analyses on large genomic and transcriptomic datasets. Further insight into the conflict and processes underlying these complex datasets is necessary to improve and develop adequate models for sequence analysis and downstream applications. To aid this effort, we developed the open source software phyparts ( https://bitbucket.org/blackrim/phyparts ), which calculates unique, conflicting, and concordant bipartitions, maps gene duplications, and outputs summary statistics such as internode certainy (ICA) scores and node-specific counts of gene duplications.
| S-EPMC4524127 | biostudies-literature
Project description:BACKGROUND:Partitioning involves estimating independent models of molecular evolution for different subsets of sites in a sequence alignment, and has been shown to improve phylogenetic inference. Current methods for estimating best-fit partitioning schemes, however, are only computationally feasible with datasets of fewer than 100 loci. This is a problem because datasets with thousands of loci are increasingly common in phylogenetics. METHODS:We develop two novel methods for estimating best-fit partitioning schemes on large phylogenomic datasets: strict and relaxed hierarchical clustering. These methods use information from the underlying data to cluster together similar subsets of sites in an alignment, and build on clustering approaches that have been proposed elsewhere. RESULTS:We compare the performance of our methods to each other, and to existing methods for selecting partitioning schemes. We demonstrate that while strict hierarchical clustering has the best computational efficiency on very large datasets, relaxed hierarchical clustering provides scalable efficiency and returns dramatically better partitioning schemes as assessed by common criteria such as AICc and BIC scores. CONCLUSIONS:These two methods provide the best current approaches to inferring partitioning schemes for very large datasets. We provide free open-source implementations of the methods in the PartitionFinder software. We hope that the use of these methods will help to improve the inferences made from large phylogenomic datasets.
| S-EPMC4012149 | biostudies-literature
Project description:In phylogenomics, incongruences between gene trees, resulting from both artifactual and biological reasons, can decrease the signal-to-noise ratio and complicate species tree inference. The amount of data handled today in classical phylogenomic analyses precludes manual error detection and removal. However, a simple and efficient way to automate the identification of outliers from a collection of gene trees is still missing. Here, we present PhylteR, a method that allows rapid and accurate detection of outlier sequences in phylogenomic datasets, i.e. species from individual gene trees that do not follow the general trend. PhylteR relies on DISTATIS, an extension of multidimensional scaling to 3 dimensions to compare multiple distance matrices at once. In PhylteR, these distance matrices extracted from individual gene phylogenies represent evolutionary distances between species according to each gene. On simulated datasets, we show that PhylteR identifies outliers with more sensitivity and precision than a comparable existing method. We also show that PhylteR is not sensitive to ILS-induced incongruences, which is a desirable feature. On a biological dataset of 14,463 genes for 53 species previously assembled for Carnivora phylogenomics, we show (i) that PhylteR identifies as outliers sequences that can be considered as such by other means, and (ii) that the removal of these sequences improves the concordance between the gene trees and the species tree. Thanks to the generation of numerous graphical outputs, PhylteR also allows for the rapid and easy visual characterization of the dataset at hand, thus aiding in the precise identification of errors. PhylteR is distributed as an R package on CRAN and as containerized versions (docker and singularity).
| S-EPMC10655845 | biostudies-literature
Project description:BackgroundCollective animal behavior, such as the flocking of birds or the shoaling of fish, has inspired a class of algorithms designed to optimize distance-based clusters in various applications, including document analysis and DNA microarrays. In a flocking model, individual agents respond only to their immediate environment and move according to a few simple rules. After several iterations the agents self-organize, and clusters emerge without the need for partitional seeds. In addition to its unsupervised nature, flocking offers several computational advantages, including the potential to reduce the number of required comparisons.FindingsIn the tool presented here, Clusterflock, we have implemented a flocking algorithm designed to locate groups (flocks) of orthologous gene families (OGFs) that share an evolutionary history. Pairwise distances that measure phylogenetic incongruence between OGFs guide flock formation. We tested this approach on several simulated datasets by varying the number of underlying topologies, the proportion of missing data, and evolutionary rates, and show that in datasets containing high levels of missing data and rate heterogeneity, Clusterflock outperforms other well-established clustering techniques. We also verified its utility on a known, large-scale recombination event in Staphylococcus aureus. By isolating sets of OGFs with divergent phylogenetic signals, we were able to pinpoint the recombined region without forcing a pre-determined number of groupings or defining a pre-determined incongruence threshold.ConclusionsClusterflock is an open-source tool that can be used to discover horizontally transferred genes, recombined areas of chromosomes, and the phylogenetic 'core' of a genome. Although we used it here in an evolutionary context, it is generalizable to any clustering problem. Users can write extensions to calculate any distance metric on the unit interval, and can use these distances to 'flock' any type of data.
| S-EPMC5078944 | biostudies-literature
Project description:For in vitro culture of plant and animal cells, one of the critical steps is to adjust the initial cell density. A typical example of this is isolated microspore culture, where specific cell densities have been determined for different species. Out of these ranges, microspore growth is not induced, or is severely reduced. A similar situation occurs in many other plant and animal cell culture systems. Traditionally, researchers have used counting chambers (hemacytometers) to calculate cell densities, but little is still known about their technical advantages. In addition, much less information is available about other, alternative methods. In this work, using isolated eggplant microspore cultures and fluorescent beads (fluorospheres) as experimental systems, we performed a comprehensive comparison of six methods to calculate cell densities: (1) a Neubauer improved hemacytometer, (2) an automated cell counter, (3) a manual-counting method, and three flow cytometry methods based on (4) autofluorescence, (5) propidium iodide staining, and (6) side scattered light (SSC).Our results show that from a technical perspective, hemacytometers are the most reasonable option for cell counting, which may explain their widely spread use. Automated cell counters represent a good compromise between precision and affordability, although with limited accuracy. Finally, the methods based on flow cytometry were, by far, the best in terms of reproducibility and agreement between them, but they showed deficient accuracy and precision.Together, our results show a thorough technical evaluation of each counting method, provide unambiguous arguments to decide which one is the most convenient for the particular case of each laboratory, and in general, shed light into the best way to determine cell densities for in vitro cell cultures. They may have an impact in such a practice not only in the context of microspore culture, but also in any other plant cell culture procedure, or in any process involving particle counting.
| S-EPMC5901878 | biostudies-other
Project description:MotivationWith the rapid growth rate of newly sequenced genomes, species tree inference from multiple genes has become a basic bioinformatics task in comparative and evolutionary biology. However, accurate species tree estimation is difficult in the presence of gene tree discordance, which is often due to incomplete lineage sorting (ILS), modelled by the multi-species coalescent. Several highly accurate coalescent-based species tree estimation methods have been developed over the last decade, including MP-EST. However, the running time for MP-EST increases rapidly as the number of species grows.ResultsWe present divide-and-conquer techniques that improve the scalability of MP-EST so that it can run efficiently on large datasets. Surprisingly, this technique also improves the accuracy of species trees estimated by MP-EST, as our study shows on a collection of simulated and biological datasets.
| S-EPMC4239662 | biostudies-literature
Project description:Background:As next generation sequence technology has advanced, there have been parallel advances in genome-scale analysis programs for determining evolutionary relationships as proxies for epidemiological relationship in public health. Most new programs skip traditional steps of ortholog determination and multi-gene alignment, instead identifying variants across a set of genomes, then summarizing results in a matrix of single-nucleotide polymorphisms or alleles for standard phylogenetic analysis. However, public health authorities need to document the performance of these methods with appropriate and comprehensive datasets so they can be validated for specific purposes, e.g., outbreak surveillance. Here we propose a set of benchmark datasets to be used for comparison and validation of phylogenomic pipelines. Methods:We identified four well-documented foodborne pathogen events in which the epidemiology was concordant with routine phylogenomic analyses (reference-based SNP and wgMLST approaches). These are ideal benchmark datasets, as the trees, WGS data, and epidemiological data for each are all in agreement. We have placed these sequence data, sample metadata, and "known" phylogenetic trees in publicly-accessible databases and developed a standard descriptive spreadsheet format describing each dataset. To facilitate easy downloading of these benchmarks, we developed an automated script that uses the standard descriptive spreadsheet format. Results:Our "outbreak" benchmark datasets represent the four major foodborne bacterial pathogens (Listeria monocytogenes, Salmonella enterica, Escherichia coli, and Campylobacter jejuni) and one simulated dataset where the "known tree" can be accurately called the "true tree". The downloading script and associated table files are available on GitHub: https://github.com/WGS-standards-and-analysis/datasets. Discussion:These five benchmark datasets will help standardize comparison of current and future phylogenomic pipelines, and facilitate important cross-institutional collaborations. Our work is part of a global effort to provide collaborative infrastructure for sequence data and analytic tools-we welcome additional benchmark datasets in our recommended format, and, if relevant, we will add these on our GitHub site. Together, these datasets, dataset format, and the underlying GitHub infrastructure present a recommended path for worldwide standardization of phylogenomic pipelines.
| S-EPMC5782805 | biostudies-literature
Project description:Rapid relaxed-clock dating methods are frequently applied to analyze phylogenomic data sets containing hundreds to thousands of sequences because of their accuracy and computational efficiency. However, the relative performance of different rapid dating methods is yet to be compared on the same data sets, and, thus, the power and pitfalls of selecting among these approaches remain unclear. We compared the accuracy, bias, and coverage probabilities of RelTime, treePL, and least-squares dating time estimates by applying them to analyze computer-simulated data sets in which evolutionary rates varied extensively among branches in the phylogeny. RelTime estimates were consistently more accurate than the other two, particularly when evolutionary rates were autocorrelated or shifted convergently among lineages. The 95% confidence intervals (CIs) around RelTime dates showed appropriate coverage probabilities (95% on average), but other methods produced rather low coverage probabilities because of overly narrow CIs of time estimates. Overall, RelTime appears to be a more efficient method for estimating divergence times for large phylogenies.
| S-EPMC8633771 | biostudies-literature
Project description:A single population of progenitor cells maintains many epithelial tissues. Transgenic mouse cell tracking has frequently been used to study the growth dynamics of competing clones in these tissues. A mathematical model (the 'single-progenitor model') has been argued to reproduce the observed progenitor dynamics accurately. This requires three parameters to describe the growth dynamics observed in transgenic mouse cell tracking-a division rate, a stratification rate and the probability of dividing symmetrically. Deriving these parameters is a time intensive and complex process. We compare the alternative strategies for analysing this source of experimental data, identifying an approximate Bayesian computation-based approach as the best in terms of efficiency and appropriate error estimation. We support our findings by explicitly modelling biological variation and consider the impact of different sampling regimes. All tested solutions are made available to allow new datasets to be analysed following our workflows. Based on our findings, we make recommendations for future experimental design.
| S-EPMC8097194 | biostudies-literature
Project description:BackgroundHospital readmissions are important patient outcomes that can be accurately captured with routinely collected administrative data. Hospital-specific readmission rates have been reported as a quality-of-care indicator. However, the extent to which these measures vary with different calculation methods is uncertain.MethodsWe identified all discharges from Ontario hospitals from 2005 to 2010 and determined whether patients died or were urgently readmitted within 30 days. For each hospital, we calculated 4 distinct observed-to-expected ratios, estimating the expected number of events using different adjustments for confounders (age and sex v. complete) and different units of analysis (all admissions v. single admission per patient).ResultsWe included 3 148 648 admissions to hospital for 1 802 704 patients in 162 hospitals. Ratios adjusted for age and sex alone had the greatest variation. Within hospitals, ranges of the 4 ratios averaged 31% of the overall estimate. Readmission ratios adjusted for age and sex showed the lowest correlation (Spearman correlation coefficient 0.48-0.68). Hospital rankings based on the different measures had an average range of 47.4 (standard deviation 32.2) out of 162.InterpretationWe found notable variation in rates of death or urgent readmission within 30 days based on the extent of adjustment for confounders and the unit of analysis. Slight changes in the methods used to calculate hospital-specific readmission rates influence their values and the consequent rankings of hospitals. Our results highlight the caution required when comparing hospital performance using rates of death or urgent readmission within 30 days.
| S-EPMC3478375 | biostudies-literature