Operations Research Methods for Estimating the Population Size of Neuron Types.
Ontology highlight
ABSTRACT:
SUBMITTER: Attili SM
PROVIDER: S-EPMC7748248 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:In many species a fundamental feature of genetic diversity is that genetic similarity decays with geographic distance; however, this relationship is often complex, and may vary across space and time. Methods to uncover and visualize such relationships have widespread use for analyses in molecular ecology, conservation genetics, evolutionary genetics, and human genetics. While several frameworks exist, a promising approach is to infer maps of how migration rates vary across geographic space. Such maps could, in principle, be estimated across time to reveal the full complexity of population histories. Here, we take a step in this direction: we present a method to infer maps of population sizes and migration rates associated with different time periods from a matrix of genetic similarity between every pair of individuals. Specifically, genetic similarity is measured by counting the number of long segments of haplotype sharing (also known as identity-by-descent tracts). By varying the length of these segments we obtain parameter estimates associated with different time periods. Using simulations, we show that the method can reveal time-varying migration rates and population sizes, including changes that are not detectable when using a similar method that ignores haplotypic structure. We apply the method to a dataset of contemporary European individuals (POPRES), and provide an integrated analysis of recent population structure and growth over the last ∼3,000 years in Europe.
| S-EPMC6347299 | biostudies-literature
Project description:It is known that the effective population size (Ne) and the mutation rate (u) vary across the genome. Here, we show that ignoring this heterogeneity may lead to biased estimates of past demography. To solve the problem, we develop new methods for jointly inferring past changes in population size and detecting variation in Ne and u between loci. These methods rely on either polymorphism data alone or both polymorphism and divergence data. In addition to inferring demography, we can use the methods to study a variety of questions: 1) comparing sex chromosomes with autosomes (for finding evidence for male-driven evolution, an unequal sex ratio, or sex-biased demographic changes) and 2) analyzing multilocus data from within autosomes or sex chromosomes (for studying determinants of variability in Ne and u). Simulations suggest that the methods can provide accurate parameter estimates and have substantial statistical power for detecting difference in Ne and u. As an example, we use the methods to analyze a polymorphism data set from Drosophila simulans. We find clear evidence for rapid population expansion. The results also indicate that the autosomes have a higher mutation rate than the X chromosome and that the sex ratio is probably female-biased. The new methods have been implemented in a user-friendly package.
| S-EPMC6409433 | biostudies-literature
Project description:Two new methods for quantifying the rapidity of action potential onset have lower relative standard deviations and better distinguish neuron cell types than current methods. Action potentials (APs) in most central mammalian neurons exhibit sharp onset dynamics. The main views explaining such an abrupt onset differ. Some studies suggest sharp onsets reflect cooperative sodium channels activation, while others suggest they reflect AP backpropagation from the axon initial segment. However, AP onset rapidity is defined subjectively in these studies, often using the slope at an arbitrary value on the phase plot. Thus, we proposed more systematic methods using the membrane potential's second-time derivative ([Formula: see text]) peak width. Here, the AP rapidity was measured for four different cortical and hippocampal neuron types using four quantification methods: the inverse of full-width at the half maximum of the [Formula: see text] peak (IFWd2), the inverse of half-width at the half maximum of the [Formula: see text] peak (IHWd2), the phase plot slope, and the error ratio method. The IFWd2 and IHWd2 methods show the smallest variation among neurons of the same type. Furthermore, the AP rapidity, using the [Formula: see text] peak width methods, significantly differentiates between different types of neurons, indicating that AP rapidity can be used to classify neuron types. The AP rapidity measured using the IFWd2 method was able to differentiate between all four neuron types analyzed. Therefore, the [Formula: see text] peak width methods provide another sensitive tool to investigate the mechanisms impacting the AP onset dynamics.
| S-EPMC8032118 | biostudies-literature
Project description:Understanding HIV incidence, the rate at which new infections occur in populations, is critical for tracking and surveillance of the epidemic. In this article, we derive methods for determining sample sizes for cross-sectional surveys to estimate incidence with sufficient precision. We further show how to specify sample sizes for two successive cross-sectional surveys to detect changes in incidence with adequate power. In these surveys biomarkers such as CD4 cell count, viral load, and recently developed serological assays are used to determine which individuals are in an early disease stage of infection. The total number of individuals in this stage, divided by the number of people who are uninfected, is used to approximate the incidence rate. Our methods account for uncertainty in the durations of time spent in the biomarker defined early disease stage. We find that failure to account for this uncertainty when designing surveys can lead to imprecise estimates of incidence and underpowered studies. We evaluated our sample size methods in simulations and found that they performed well in a variety of underlying epidemics. Code for implementing our methods in R is available with this article at the Biometrics website on Wiley Online Library.
| S-EPMC4715554 | biostudies-literature
Project description:Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.
| S-EPMC7042749 | biostudies-literature
Project description:By combining high-throughput sequencing (HTS) with experimental evolution, we can observe the within-host dynamics of pathogen variants of biomedical or ecological interest. We studied the evolutionary dynamics of five variants of Potato virus Y (PVY) in 15 doubled-haploid lines of pepper. All plants were inoculated with the same mixture of virus variants and variant frequencies were determined by HTS in eight plants of each pepper line at each of six sampling dates. We developed a method for estimating the intensities of selection and genetic drift in a multi-allelic Wright-Fisher model, applicable whether these forces are strong or weak, and in the absence of neutral markers. This method requires variant frequency determination at several time points, in independent hosts. The parameters are the selection coefficients for each PVY variant and four effective population sizes Ne at different time-points of the experiment. Numerical simulations of asexual haploid Wright-Fisher populations subjected to contrasting genetic drift (Ne ? [10, 2000]) and selection (|s| ? [0, 0.15]) regimes were used to validate the method proposed. The experiment in closely related pepper host genotypes revealed that viruses experienced a considerable diversity of selection and genetic drift regimes. The resulting variant dynamics were accurately described by Wright-Fisher models. The fitness ranks of the variants were almost identical between host genotypes. By contrast, the dynamics of Ne were highly variable, although a bottleneck was often identified during the systemic movement of the virus. We demonstrated that, for a fixed initial PVY population, virus effective population size is a heritable trait in plants. These findings pave the way for the breeding of plant varieties exposing viruses to stronger genetic drift, thereby slowing virus adaptation.
| S-EPMC5720836 | biostudies-literature
Project description:Coalescent theory combined with statistical modeling allows us to estimate effective population size fluctuations from molecular sequences of individuals sampled from a population of interest. When sequences are sampled serially through time and the distribution of the sampling times depends on the effective population size, explicit statistical modeling of sampling times improves population size estimation. Previous work assumed that the genealogy relating sampled sequences is known and modeled sampling times as an inhomogeneous Poisson process with log-intensity equal to a linear function of the log-transformed effective population size. We improve this approach in two ways. First, we extend the method to allow for joint Bayesian estimation of the genealogy, effective population size trajectory, and other model parameters. Next, we improve the sampling time model by incorporating additional sources of information in the form of time-varying covariates. We validate our new modeling framework using a simulation study and apply our new methodology to analyses of population dynamics of seasonal influenza and to the recent Ebola virus outbreak in West Africa.
| S-EPMC7580988 | biostudies-literature
Project description:Technological and methodological advances have facilitated the use of genetic data to infer census population size (Nc) in natural populations, particularly where traditional mark-and-recapture is challenging. The effective number of breeders (Nb) describes how many adults effectively contribute to a cohort and is often correlated with Nc. Predicting Nc from Nb or vice versa in species with overlapping generations has important implications for conservation by permitting (i) estimation of the more difficult to quantify variable and (ii) inferences of Nb/Nc relationships in related species lacking data. We quantitatively synthesized Nb/Nc relationships in three salmonid fishes where sufficient data have recently accumulated. Mixed-effects models were analysed in which each variable was included as a dependent variable or predictor term (Nb from Nc and vice versa). Species-dependent Nb/Nc slope estimates were significantly positive in two of three species. Variation in species slopes was likely due to varying life histories and reinforce caution when inferring Nb/Nc from taxonomically related species. Models provided maximum probable estimates for Nb and Nc for two species. However, study, population and year effects explained substantial amounts of variation (39%-57%). Consequently, prediction intervals were wide and included or were close to zero for all population sizes and species; model predictive utility was limited. Cost-benefit trade-offs when estimating Nb and/or Nc were also discussed using a real-world system example. Our findings based on salmonids suggest that no short cuts currently exist when estimating population size and researchers should focus on quantifying the variable of interest or be aware of caveats when inferring the desired variable because of cost or logistics. We caution that the salmonid species examined share life-history traits that may obscure relationships between Nb and Nc. Sufficient data on other taxa were unavailable; additional research examining Nb/Nc relationships in species with potentially relevant life-history trait differences (e.g., differing survival curves) is needed.
| S-EPMC5680432 | biostudies-literature
Project description:Estimates of wildlife population size are critical for conservation and management, but accurate estimates are difficult to obtain for many species. Several methods have recently been developed that estimate abundance using kinship relationships observed in genetic samples, particularly parent-offspring pairs. While these methods are similar to traditional Capture-Mark-Recapture, they do not need physical recapture, as individuals are considered recaptured if a sample contains one or more close relatives. This makes methods based on genetically-identified parent-offspring pairs particularly interesting for species for which releasing marked animals back into the population is not desirable or not possible (e.g., harvested fish or game species). However, while these methods have successfully been applied in commercially important fish species, in the absence of life-history data, they are making several assumptions unlikely to be met for harvested terrestrial species. They assume that a sample contains only one generation of parents and one generation of juveniles of the year, while more than two generations can coexist in the hunting bags of long-lived species, or that the sampling probability is the same for each individual, an assumption that is violated when fecundity and/or survival depend on sex or other individual traits. In order to assess the usefulness of kin-based methods to estimate population sizes of terrestrial game species, we simulated population pedigrees of two different species with contrasting demographic strategies (wild boar and red deer), applied four different methods and compared the accuracy and precision of their estimates. We also performed a sensitivity analysis, simulating population pedigrees with varying fecundity characteristics and various levels of harvesting to identify optimal conditions of applicability of each method. We showed that all these methods reached the required levels of accuracy and precision to be effective in wildlife management under simulated circumstances (i.e., for species within a given range of fecundity and for a given range of sampling intensity), while being robust to fecundity variation. Despite the potential usefulness of the methods for terrestrial game species, care is needed as several biases linked to hunting practices still need to be investigated (e.g., when hunting bags are biased toward a particular group of individuals).
| S-EPMC10105560 | biostudies-literature
Project description:BackgroundThe intracellular accommodation of arbuscular mycorrhizal (AM) fungi involves a profound molecular reprogramming of the host cell architecture and metabolism, based on the activation of a symbiotic signaling pathway. In analogy with other plant biotrophs, AM fungi are reported to trigger cell cycle reactivation in their host tissues, possibly in support of the enhanced metabolic demand required for the symbiosis.ResultsWe here compare the efficiency of three Fiji/ImageJ image analysis plugins in localizing and quantifying the increase in nuclear size - a hallmark of recursive events of endoreduplication - in M. truncatula roots colonized by the AM fungus Gigaspora margarita. All three approaches proved to be versatile and upgradeable, allowing the investigation of nuclear changes in a complex tissue; 3D Object Counter provided more detailed information than both TrackMate and Round Surface Detector plugins. On this base we challenged 3D Object Counter with two case studies: verifying the lack of endoreduplication-triggering responses in Medicago truncatula mutants with a known non-symbiotic phenotype; and analysing the correlation in space and time between the induction of cortical cell division and endoreduplication upon AM colonization. Both case studies revealed important biological aspects. Mutant phenotype analyses have demonstrated that the knock-out mutation of different key genes in the symbiotic signaling pathway block AM-associated endoreduplication. Furthermore, our data show that cell divisions occur during initial stages of root colonization and are followed by recursive activation of the endocycle in preparation for arbuscule accommodation.ConclusionsIn conclusion, our results indicate 3D Object Counter as the best performing Fiji/ImageJ image analysis script in plant root thick sections and its application highlighted endoreduplication as a major feature of the AM pre-penetration response in root cortical cells.
| S-EPMC6500585 | biostudies-literature