A new resolution function to evaluate tree shape statistics.
Ontology highlight
ABSTRACT:
SUBMITTER: Hayati M
PROVIDER: S-EPMC6874070 | biostudies-literature | 2019
REPOSITORIES: biostudies-literature
Publications
A new resolution function to evaluate tree shape statistics.
Hayati Maryam M Shadgar Bita B Chindelevitch Leonid L
PloS one 20191121 11
Phylogenetic trees are frequently used in biology to study the relationships between a number of species or organisms. The shape of a phylogenetic tree contains useful information about patterns of speciation and extinction, so powerful tools are needed to investigate the shape of a phylogenetic tree. Tree shape statistics are a common approach to quantifying the shape of a phylogenetic tree by encoding it with a single number. In this article, we propose a new resolution function to evaluate th ...[more]
Similar Datasets
Project description:The shape of phylogenetic trees can be used to gain evolutionary insights. A tree's shape specifies the connectivity of a tree, while its branch lengths reflect either the time or genetic distance between branching events; well-known measures of tree shape include the Colless and Sackin imbalance, which describe the asymmetry of a tree. In other contexts, network science has become an important paradigm for describing structural features of networks and using them to understand complex systems, ranging from protein interactions to social systems. Network science is thus a potential source of many novel ways to characterize tree shape, as trees are also networks. Here, we tailor tools from network science, including diameter, average path length, and betweenness, closeness, and eigenvector centrality, to summarize phylogenetic tree shapes. We thereby propose tree shape summaries that are complementary to both asymmetry and the frequencies of small configurations. These new statistics can be computed in linear time and scale well to describe the shapes of large trees. We apply these statistics, alongside some conventional tree statistics, to phylogenetic trees from three very different viruses (HIV, dengue fever and measles), from the same virus in different epidemiological scenarios (influenza A and HIV) and from simulation models known to produce trees with different shapes. Using mutual information and supervised learning algorithms, we find that the statistics adapted from network science perform as well as or better than conventional statistics. We describe their distributions and prove some basic results about their extreme values in a tree. We conclude that network science-based tree shape summaries are a promising addition to the toolkit of tree shape features. All our shape summaries, as well as functions to select the most discriminating ones for two sets of trees, are freely available as an R package at http://github.com/Leonardini/treeCentrality.
| S-EPMC8699983 | biostudies-literature
Project description:BackgroundHost population structure is a key determinant of pathogen and infectious disease transmission patterns. Pathogen phylogenetic trees are useful tools to reveal the population structure underlying an epidemic. Determining whether a population is structured or not is useful in informing the type of phylogenetic methods to be used in a given study. We employ tree statistics derived from phylogenetic trees and machine learning classification techniques to reveal an underlying population structure.ResultsIn this paper, we simulate phylogenetic trees from both structured and non-structured host populations. We compute eight statistics for the simulated trees, which are: the number of cherries; Sackin, Colless and total cophenetic indices; ladder length; maximum depth; maximum width, and width-to-depth ratio. Based on the estimated tree statistics, we classify the simulated trees as from either a non-structured or a structured population using the decision tree (DT), K-nearest neighbor (KNN) and support vector machine (SVM). We incorporate the basic reproductive number ([Formula: see text]) in our tree simulation procedure. Sensitivity analysis is done to investigate whether the classifiers are robust to different choice of model parameters and to size of trees. Cross-validated results for area under the curve (AUC) for receiver operating characteristic (ROC) curves yield mean values of over 0.9 for most of the classification models.ConclusionsOur classification procedure distinguishes well between trees from structured and non-structured populations using the classifiers, the two-sample Kolmogorov-Smirnov, Cucconi and Podgor-Gastwirth tests and the box plots. SVM models were more robust to changes in model parameters and tree size compared to KNN and DT classifiers. Our classification procedure was applied to real -world data and the structured population was revealed with high accuracy of [Formula: see text] using SVM-polynomial classifier.
| S-EPMC8579572 | biostudies-literature
Project description:Amplicon-based analysis of 16S rRNA genes and transcripts was used to assess the effect of tree species composition on soil bacterial community structure and function in a temperate deciduous forest. Samples were collected from mono and mixed stands of Fagus sylvatica (beech), Carpinus betulus (hornbeam), Tilia sp. (lime), and Quercus sp. (oak) in spring, summer, and autumn. Soil bacterial community exhibited similar taxonomic composition at total (DNA-based) and potentially active community (RNA-based) level, with fewer taxa present at active community level. Members of Rhizobiales dominated at both total and active bacterial community level, followed by members of Acidobacteriales, Solibacterales, Rhodospirillales, and Xanthomonadales. Bacterial communities at total and active community level showed a significant positive correlation with tree species identity (mono stands) and to a lesser extent with tree species richness (mixed stands). Approximately 58 and 64% of indicator operational taxonomic units (OTUs) showed significant association with only one mono stand at total and active community level, respectively, indicating a strong impact of tree species on soil bacterial community composition. Soil C/N ratio, pH, and P content similarly exhibited a significant positive correlation with soil bacterial communities, which was attributed to direct and indirect effects of forest stands. Seasonality was the strongest driver of predicted metabolic functions related to C fixation and degradation, and N metabolism. Carbon and nitrogen metabolic processes were significantly abundant in spring, while C degradation gene abundances increased from summer to autumn, corresponding to increased litterfall and decomposition. The results revealed that in a spatially homogenous forest soil, tree species diversity and richness are dominant drivers of structure and composition in soil bacterial communities.
| S-EPMC6629791 | biostudies-literature
Project description:This study proposed a novel methodology to classify the shape of gaps using landscape indices and multivariate statistics. Patch-level indices were used to collect the qualified shape and spatial configuration characteristics for canopy gaps in the Lienhuachih Experimental Forest in Taiwan in 1998 and 2002. Non-hierarchical cluster analysis was used to assess the optimal number of gap clusters and canonical discriminant analysis was used to generate the discriminant functions for canopy gap classification. The gaps for the two periods were optimally classified into three categories. In general, gap type 1 had a more complex shape, gap type 2 was more elongated and gap type 3 had the largest gaps that were more regular in shape. The results were evaluated using Wilks' lambda as satisfactory (p?<?0.001). The agreement rate of confusion matrices exceeded 96%. Differences in gap characteristics between the classified gap types that were determined using a one-way ANOVA showed a statistical significance in all patch indices (p?=?0.00), except for the Euclidean nearest neighbor distance (ENN) in 2002. Taken together, these results demonstrated the feasibility and applicability of the proposed methodology to classify the shape of a gap.
| S-EPMC5128799 | biostudies-other
Project description:The estimation of visual motion has long been studied as a paradigmatic neural computation, and multiple models have been advanced to explain behavioral and neural responses to motion signals. A broad class of models, originating with the Reichardt correlator model, proposes that animals estimate motion by computing a temporal cross-correlation of light intensities from two neighboring points in visual space. These models provide a good description of experimental data in specific contexts but cannot explain motion percepts in stimuli lacking pairwise correlations. Here, we develop a theoretical formalism that can accommodate diverse stimuli and behavioral goals. To achieve this, we treat motion estimation as a problem of Bayesian inference. Pairwise models emerge as one component of the generalized strategy for motion estimation. However, correlation functions beyond second order enable more accurate motion estimation. Prior expectations that are asymmetric with respect to bright and dark contrast use correlations of both even and odd orders, and we show that psychophysical experiments using visual stimuli with symmetric probability distributions for contrast cannot reveal whether the subject uses odd-order correlators for motion estimation. This result highlights a gap in previous experiments, which have largely relied on symmetric contrast distributions. Our theoretical treatment provides a natural interpretation of many visual motion percepts, indicates that motion estimation should be revisited using a broader class of stimuli, demonstrates how correlation-based motion estimation is related to stimulus statistics, and provides multiple experimentally testable predictions.
| S-EPMC3150910 | biostudies-literature
Project description:Across the Tree of Life, there are dramatic differences in species numbers among groups. However, the factors that explain the differences among the deepest branches have remained unknown. We tested whether multicellularity and sexual reproduction might explain these patterns, since the most species-rich groups share these traits. We found that groups with multicellularity and sexual reproduction have accelerated rates of species proliferation (diversification), and that multicellularity has a stronger effect than sexual reproduction. Patterns of species richness among clades are then strongly related to these differences in diversification rates. Taken together, these results help explain patterns of biodiversity among groups of organisms at the very broadest scales. They may also help explain the mysterious preponderance of sexual reproduction among species (the 'paradox of sex') by showing that organisms with sexual reproduction proliferate more rapidly.
| S-EPMC8316805 | biostudies-literature
Project description:BackgroundMarkers of left atrial (LA) shape may improve the prediction of postablation outcomes in atrial fibrillation (AF). Correlations to LA volume and AF persistence limit their incremental value over current clinical predictors.ObjectiveTo develop a shape score independent from AF persistence and LA volume using shape-based statistics, and to test its ability to predict postablation outcome.MethodsPreablation computed tomography (CT) images from 141 patients with paroxysmal (57%) or persistent (43%) AF were segmented. Deformation of an average LA shape into each patient encoded patient-specific shape. Local analysis investigates regional differences between patient groups. Linear regression was used to remove shape variations related to LA volume and AF persistence, and to build a shape score to predict postablation outcome. Cross-validation was performed to evaluate its accuracy.ResultsAblation failure rate was 23% over a median 12-month follow-up. Regions associated with ablation failure mostly consisted of a large area on posteroinferior LA, mitral isthmus, and left inferior vein. On univariate analysis, strongest predictors were AF persistence (P = .005), LA indexed volume (P = .02), and the proposed shape score (P = .001). On multivariate analysis, all 3 were independent predictors of ablation failure, with the LA shape score showing the highest predictive value (odds ratio [OR] = 6.2 [2.5-15.8], P < .001), followed by LA indexed volume (OR = 3.1 [1.2-7.9], P = .019) and AF persistence (OR = 2.9 [1.2-7.6], P = .022).ConclusionPosteroinferior LA, mitral isthmus, and left inferior vein are the regions whose shape have the highest impact on outcome. LA shape predicts AF ablation failure independently from, and more accurately than, atrial volume and AF persistence.
| S-EPMC8703187 | biostudies-literature
Project description:Cophylogeny is the congruence of phylogenetic relationships between two different groups of organisms due to their long-term interaction. We investigated the use of tree shape distance measures to quantify the degree of cophylogeny. We implemented a reverse-time simulation model of pathogen phylogenies within a fixed host tree, given cospeciation probability, host switching, and pathogen speciation rates. We used this model to evaluate 18 distance measures between host and pathogen trees including two kernel distances that we developed for labeled and unlabeled trees, which use branch lengths and accommodate different size trees. Finally, we used these measures to revisit published cophylogenetic studies, where authors described the observed associations as representing a high or low degree of cophylogeny. Our simulations demonstrated that some measures are more informative than others with respect to specific coevolution parameters especially when these did not assume extreme values. For real datasets, trees' associations projection revealed clustering of high concordance studies suggesting that investigators are describing it in a consistent way. Our results support the hypothesis that measures can be useful for quantifying cophylogeny. This motivates their usage in the field of coevolution and supports the development of simulation-based methods, i.e., approximate Bayesian computation, to estimate the underlying coevolutionary parameters.
| S-EPMC6618157 | biostudies-literature
Project description:Single nucleotide polymorphisms (SNPs) are useful markers for phylogenetic studies owing in part to their ubiquity throughout the genome and ease of collection. Restriction site associated DNA sequencing (RADseq) methods are becoming increasingly popular for SNP data collection, but an assessment of the best practises for using these data in phylogenetics is lacking. We use computer simulations, and new double digest RADseq (ddRADseq) data for the lizard family Phrynosomatidae, to investigate the accuracy of RAD loci for phylogenetic inference. We compare the two primary ways RAD loci are used during phylogenetic analysis, including the analysis of full sequences (i.e., SNPs together with invariant sites), or the analysis of SNPs on their own after excluding invariant sites. We find that using full sequences rather than just SNPs is preferable from the perspectives of branch length and topological accuracy, but not of computational time. We introduce two new acquisition bias corrections for dealing with alignments composed exclusively of SNPs, a conditional likelihood method and a reconstituted DNA approach. The conditional likelihood method conditions on the presence of variable characters only (the number of invariant sites that are unsampled but known to exist is not considered), while the reconstituted DNA approach requires the user to specify the exact number of unsampled invariant sites prior to the analysis. Under simulation, branch length biases increase with the amount of missing data for both acquisition bias correction methods, but branch length accuracy is much improved in the reconstituted DNA approach compared to the conditional likelihood approach. Phylogenetic analyses of the empirical data using concatenation or a coalescent-based species tree approach provide strong support for many of the accepted relationships among phrynosomatid lizards, suggesting that RAD loci contain useful phylogenetic signal across a range of divergence times despite the presence of missing data. Phylogenetic analysis of RAD loci requires careful attention to model assumptions, especially if downstream analyses depend on branch lengths.
| S-EPMC4604835 | biostudies-literature
Project description:Diversifying planted forests by increasing genetic and species diversity is often promoted as a method to improve forest resilience to climate change and reduce pest and pathogen damage. In this study, we used a young tree diversity experiment replicated at two sites in the UK to study the impacts of tree diversity and tree provenance (geographic origin) on the oak (Quercus robur) insect herbivore community and a specialist biotrophic pathogen, oak powdery mildew. Local UK, French, and Italian provenances were planted in monocultures, provenance mixtures, and species mixes, allowing us to test whether: (a) local and nonlocal provenances differ in their insect herbivore and pathogen communities, and (b) admixing trees leads to associational effects on insect herbivore and pathogen damage. Tree diversity had variable impacts on foliar organisms across sites and years, suggesting that diversity effects can be highly dependent on environmental context. Provenance identity impacted upon both herbivores and powdery mildew, but we did not find consistent support for the local adaptation hypothesis for any group of organisms studied. Independent of provenance, we found tree vigor traits (shoot length, tree height) and tree apparency (the height of focal trees relative to their surroundings) were consistent positive predictors of powdery mildew and insect herbivory. Synthesis. Our results have implications for understanding the complex interplay between tree identity and diversity in determining pest damage, and show that tree traits, partially influenced by tree genotype, can be important drivers of tree pest and pathogen loads.
| S-EPMC6686283 | biostudies-literature