On estimating evolutionary probabilities of population variants.
Ontology highlight
ABSTRACT:
SUBMITTER: Patel R
PROVIDER: S-EPMC6593550 | biostudies-literature | 2019 Jun
REPOSITORIES: biostudies-literature
Publications
On estimating evolutionary probabilities of population variants.
BMC evolutionary biology 20190625 1
<h4>Background</h4>The evolutionary probability (EP) of an allele in a DNA or protein sequence predicts evolutionarily permissible (ePerm; EP ≥ 0.05) and forbidden (eForb; EP < 0.05) variants. EP of an allele represents an independent evolutionary expectation of observing an allele in a population based solely on the long-term substitution patterns captured in a multiple sequence alignment. In the neutral theory, EP and population frequencies can be compared to identify neutral and non-neutral a ...[more]
Similar Datasets
Project description:The global rise in nature-based recreation increases the need for research on visitor activity use and interaction especially for multi-use trail systems. Conflict often arises during negatively perceived physical encounters (i.e., direct observation) of different user groups. Our study addresses these encounters on a winter multi-use refuge in Fairbanks, Alaska. Our goal was to develop a method that generates spatially and temporally explicit estimates of trail occupancy and encounter probabilities among different user groups. We used trail cameras with optic alteration to protect individual identity. We monitored winter recreational activity from November 2019 to April 2020 (n = 133 days) and sorted users into three user groups: 1) motor-powered, 2) dog-powered, and 3) human-powered. We calculated the total number of occurrences and proportion of activity across all user groups at each camera location. We identified hotspots of activity overlap (e.g., near trail access points), and peak times (14:01-15:00), days (Saturdays and Sundays), and months (December, February, and March) that may have had higher potential for physical encounters and conflict. We used multiplication and addition probability rules to estimate two probabilities: 1) the probability of user groups occupying individual trail segments, and 2) the probability of encounter between different user groups. We scaled up these probability estimates both temporally (hourly and daily) and spatially (refuge quadrant and refuge-wide). Researchers can adapt our novel method to any recreational trail system to identify locations with potential for congestion and conflict. This method can help inform management that improves visitor experience and overall trail user satisfaction.Management implicationsWe provide managers of recreational trail systems with a quantitative, objective, and noninvasive method to monitor activity among trail user groups. This method can be altered both spatially and temporally to fit any recreational trail system's research questions. These questions may involve congestion, trail carrying capacity, or user group and wildlife encounters. Our method advances current knowledge of trail use dynamics by quantifying the extent of activity overlap between different user groups that may be prone to conflict. Managers can use this information to incorporate relevant management strategies to mitigate congestion and conflict for their own recreational trail system.
| S-EPMC10215815 | biostudies-literature
Project description:Motivation: Computational inference of genome organization based on Hi-C sequencing has greatly aided the understanding of chromatin and nuclear organization in three dimensions (3D). However, existing computational methods fail to address the cell population heterogeneity. Here we describe a probabilistic modeling-based method called CscoreTool-M that infers multiple 3D genome sub-compartments from Hi-C data. Results: The compartment scores inferred using CscoreTool-M represents the probability of a genomic region locating in a specific sub-compartment. Compared to published methods, CscoreTool-M is more accurate in inferring sub-compartments corresponding to both active and repressed chromatin. The compartment scores calculated by CscoreTool-M also help to quantify the levels of heterogeneity in sub-compartment localization within cell populations. By comparing proliferating cells and terminally differentiated non-proliferating cells, we show that the proliferating cells have higher genome organization heterogeneity, which is likely caused by cells at different cell-cycle stages. By analyzing 10 sub-compartments, we found a sub-compartment containing chromatin potentially related to the early-G1 chromatin regions proximal to the nuclear lamina in HCT116 cells, suggesting the method can deconvolve cell cycle stage-specific genome organization among asynchronously dividing cells. Finally, we show that CscoreTool-M can identify sub-compartments that contain genes enriched in housekeeping or cell-type-specific functions. Availability: https://github.com/scoutzxb/CscoreTool-M
2023-05-31 | GSE224282 | GEO
Project description:Ecological and evolutionary change is generated by variation in individual performance. Biologists have consequently long been interested in decomposing change measured at the population level into contributions from individuals, the traits they express and the alleles they carry. We present a novel method of estimating individual contributions to population growth and changes in distributions of quantitative traits and alleles. An individual's contribution to population growth is an individual's realized annual fitness. We demonstrate how the quantities we develop can be used to address a range of empirical questions, and provide an application to a detailed dataset of Soay sheep. The approach provides results that are consistent with those obtained using lifetime estimates of individual performance, yet is substantially more powerful as it allows lifetime performance to be decomposed into annual survival and fecundity contributions.
| S-EPMC1560056 | biostudies-literature
Project description:Genome-wide Association Studies (GWAS) result in millions of summary statistics ("z-scores") for single nucleotide polymorphism (SNP) associations with phenotypes. These rich datasets afford deep insights into the nature and extent of genetic contributions to complex phenotypes such as psychiatric disorders, which are understood to have substantial genetic components that arise from very large numbers of SNPs. The complexity of the datasets, however, poses a significant challenge to maximizing their utility. This is reflected in a need for better understanding the landscape of z-scores, as such knowledge would enhance causal SNP and gene discovery, help elucidate mechanistic pathways, and inform future study design. Here we present a parsimonious methodology for modeling effect sizes and replication probabilities, relying only on summary statistics from GWAS substudies, and a scheme allowing for direct empirical validation. We show that modeling z-scores as a mixture of Gaussians is conceptually appropriate, in particular taking into account ubiquitous non-null effects that are likely in the datasets due to weak linkage disequilibrium with causal SNPs. The four-parameter model allows for estimating the degree of polygenicity of the phenotype and predicting the proportion of chip heritability explainable by genome-wide significant SNPs in future studies with larger sample sizes. We apply the model to recent GWAS of schizophrenia (N = 82,315) and putamen volume (N = 12,596), with approximately 9.3 million SNP z-scores in both cases. We show that, over a broad range of z-scores and sample sizes, the model accurately predicts expectation estimates of true effect sizes and replication probabilities in multistage GWAS designs. We assess the degree to which effect sizes are over-estimated when based on linear-regression association coefficients. We estimate the polygenicity of schizophrenia to be 0.037 and the putamen to be 0.001, while the respective sample sizes required to approach fully explaining the chip heritability are 10(6) and 10(5). The model can be extended to incorporate prior knowledge such as pleiotropy and SNP annotation. The current findings suggest that the model is applicable to a broad array of complex phenotypes and will enhance understanding of their genetic architectures.
| S-EPMC4754432 | biostudies-literature
Project description:Understanding the biology and conducting effective conservation of migratory species requires an understanding of migratory connectivity - the geographic linkages of populations between stages of the annual cycle. Unfortunately, for most species, we are lacking such information. The North American Bird Banding Laboratory (BBL) houses an extensive database of marking, recaptures and recoveries, and such data could provide migratory connectivity information for many species. To date, however, few species have been analyzed for migratory connectivity largely because heterogeneous re-encounter probabilities make interpretation problematic. We accounted for regional variation in re-encounter probabilities by borrowing information across species and by using effort covariates on recapture and recovery probabilities in a multistate capture-recapture and recovery model. The effort covariates were derived from recaptures and recoveries of species within the same regions. We estimated the migratory connectivity for three tern species breeding in North America and over-wintering in the tropics, common (Sterna hirundo), roseate (Sterna dougallii), and Caspian terns (Hydroprogne caspia). For western breeding terns, model-derived estimates of migratory connectivity differed considerably from those derived directly from the proportions of re-encounters. Conversely, for eastern breeding terns, estimates were merely refined by the inclusion of re-encounter probabilities. In general, eastern breeding terns were strongly connected to eastern South America, and western breeding terns were strongly linked to the more western parts of the nonbreeding range under both models. Through simulation, we found this approach is likely useful for many species in the BBL database, although precision improved with higher re-encounter probabilities and stronger migratory connectivity. We describe an approach to deal with the inherent biases in BBL banding and re-encounter data to demonstrate that this large dataset is a valuable source of information about the migratory connectivity of the birds of North America.
| S-EPMC4063466 | biostudies-literature
Project description:BackgroundWe review three common methods to estimate predicted probabilities following confounder-adjusted logistic regression: marginal standardization (predicted probabilities summed to a weighted average reflecting the confounder distribution in the target population); prediction at the modes (conditional predicted probabilities calculated by setting each confounder to its modal value); and prediction at the means (predicted probabilities calculated by setting each confounder to its mean value). That each method corresponds to a different target population is underappreciated in practice. Specifically, prediction at the means is often incorrectly interpreted as estimating average probabilities for the overall study population, and furthermore yields nonsensical estimates in the presence of dichotomous confounders. Default commands in popular statistical software packages often lead to inadvertent misapplication of prediction at the means.MethodsUsing an applied example, we demonstrate discrepancies in predicted probabilities across these methods, discuss implications for interpretation and provide syntax for SAS and Stata.ResultsMarginal standardization allows inference to the total population from which data are drawn. Prediction at the modes or means allows inference only to the relevant stratum of observations. With dichotomous confounders, prediction at the means corresponds to a stratum that does not include any real-life observations.ConclusionsMarginal standardization is the appropriate method when making inference to the overall population. Other methods should be used with caution, and prediction at the means should not be used with binary confounders. Stata, but not SAS, incorporates simple methods for marginal standardization.
| S-EPMC4052139 | biostudies-literature
Project description:BackgroundWe recently reported about the derivation of a diagnostic probability function for acute coronary syndrome (ACS). The present study aims to validate the probability function as a rule-out criterion in a new sample of patients.Methods186 patients presenting with chest pain and/or dyspnea at one of the three participating hospitals' emergency rooms in Switzerland were included in the study. In these patients, information on a set of pre-specified variables was collected and a predicted probability of ACS was calculated for each patient. Approximately two weeks after the initial visit in the emergency room, patients were contacted by phone to assess whether a diagnosis of ACS was established.ResultsOf the 186 patients included in the study, 31 (17%) had an acute coronary syndrome. A risk probability for ACS below 2% was considered a rule-out criterion for ACS, leading to a sensitivity of 87% and a specificity of 17% of the rule. The characteristics of the study patients were compared to the cases from which the probability function was derived, and considerable deviations were found in some of the variables.ConclusionsThe proposed probability function, with a 2% cut-off for ruling out ACS works quite well if the patient data lie within the ranges of values of the original vignettes. If the observations deviate too much from these ranges, the predicted probabilities for ACS should be seen with caution.
| S-EPMC4289321 | biostudies-literature
Project description:BackgroundGenomic DNA has been shaped by mutational processes through evolution. The cellular machinery for error correction and repair has left its marks in the nucleotide composition along with structural and functional constraints. Therefore, the probability of observing a base in a certain position in the human genome is highly context-dependent.ResultsHere we develop context-dependent nucleotide models. We first investigate models of nucleotides conditioned on sequence context. We develop a bidirectional Markov model that use an average of the probability from a Markov model applied to both strands of the sequence and thus depends on up to 14 bases to each side of the nucleotide. We show how the genome predictability varies across different types of genomic regions. Surprisingly, this model can predict a base from its context with an average of more than 50% accuracy. For somatic variants we show a tendency towards higher probability for the variant base than for the reference base. Inspired by DNA substitution models, we develop a model of mutability that estimates a mutation matrix (called the alpha matrix) on top of the nucleotide distribution. The alpha matrix can be estimated from a much smaller context than the nucleotide model, but the final model will still depend on the full context of the nucleotide model. With the bidirectional Markov model of order 14 and an alpha matrix dependent on just one base to each side, we obtain a model that compares well with a model of mutability that estimates mutation probabilities directly conditioned on three nucleotides to each side. For somatic variants in particular, our model fits better than the simpler model. Interestingly, the model is not very sensitive to the size of the context for the alpha matrix.ConclusionsOur study found strong context dependencies of nucleotides in the human genome. The best model uses a context of 14 nucleotides to each side. Based on these models, a substitution model was constructed that separates into the context model and a matrix dependent on a small context. The model fit somatic variants particularly well.
| S-EPMC8802520 | biostudies-literature
Project description:Some individuals die before a trait is measured or expressed (the invisible fraction), and some relevant traits are not measured in any individual (missing traits). This paper discusses how these concepts can be cast in terms of missing data problems from statistics. Using missing data theory, I show formally the conditions under which a valid evolutionary inference is possible when the invisible fraction and/or missing traits are ignored. These conditions are restrictive and unlikely to be met in even the most comprehensive long-term studies. When these conditions are not met, many selection and quantitative genetic parameters cannot be estimated accurately unless the missing data process is explicitly modelled. Surprisingly, this does not seem to have been attempted in evolutionary biology. In the case of the invisible fraction, viability selection and the missing data process are often intimately linked. In such cases, models used in survival analysis can be extended to provide a flexible and justified model of the missing data mechanism. Although missing traits pose a more difficult problem, important biological parameters can still be estimated without bias when appropriate techniques are used. This is in contrast to current methods which have large biases and poor precision. Generally, the quantitative genetic approach is shown to be superior to phenotypic studies of selection when invisible fractions or missing traits exist because part of the missing information can be recovered from relatives.
| S-EPMC2596846 | biostudies-other
Project description:Maximum running speed is an important locomotor parameter for many animals-predators as well as prey-and is thus of interest to palaeobiologists wishing to reconstruct the behavioural ecology of extinct species. A variety of approaches have been tried in the past including anatomical comparisons, bone scaling and strength, safety factors and ground reaction force analyses. However, these approaches are all indirect and an alternative approach is to create a musculoskeletal model of the animal and see how fast it can run. The major advantage of this approach is that all assumptions about the animal's morphology and physiology are directly addressed, whereas the exact same assumptions are hidden in the indirect approaches. In this paper, we present simple musculoskeletal models of three extant and five extinct bipedal species. The models predict top speed in the extant species with reasonably good agreement with accepted values, so we conclude that the values presented for the five extinct species are reasonable predictions given the modelling assumptions made. Improved musculoskeletal models and better estimates of soft tissue parameters will produce more accurate values. Limited sensitivity analysis is performed on key muscle parameters but there is considerable scope for extending this in the future.
| S-EPMC2279215 | biostudies-other