SARS-CoV-2 serology across scales: a framework for unbiased seroprevalence estimation incorporating antibody kinetics and epidemic recency.
Ontology highlight
ABSTRACT:
SUBMITTER: Takahashi S
PROVIDER: S-EPMC8452112 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Serosurveys are a key resource for measuring severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) population exposure. A growing body of evidence suggests that asymptomatic and mild infections (together making up over 95% of all infections) are associated with lower antibody titers than severe infections. Antibody levels also peak a few weeks after infection and decay gradually. We developed a statistical approach to produce estimates of cumulative incidence from raw seroprevalence survey results that account for these sources of spectrum bias. We incorporate data on antibody responses on multiple assays from a postinfection longitudinal cohort, along with epidemic time series to account for the timing of a serosurvey relative to how recently individuals may have been infected. We applied this method to produce estimates of cumulative incidence from 5 large-scale SARS-CoV-2 serosurveys across different settings and study designs. We identified substantial differences between raw seroprevalence and cumulative incidence of over 2-fold in the results of some surveys, and we provide a tool for practitioners to generate cumulative incidence estimates with preset or custom parameter values. While unprecedented efforts have been launched to generate SARS-CoV-2 seroprevalence estimates over this past year, interpretation of results from these studies requires properly accounting for both population-level epidemiologic context and individual-level immune dynamics.
| S-EPMC10472487 | biostudies-literature
Project description:Timescales characterize the pace of change for many dynamic processes in nature. Timescales are usually estimated by fitting the exponential decay of data autocorrelation in the time or frequency domain. Here we show that this standard procedure often fails to recover the correct timescales due to a statistical bias arising from the finite sample size. We develop an alternative approach which estimates timescales by fitting the sample autocorrelation or power spectrum with a generative model based on a mixture of Ornstein-Uhlenbeck (OU) processes using adaptive approximate Bayesian computations (aABC). Our method accounts for finite sample size and noise in data and returns a posterior distribution of timescales that quantifies the estimation uncertainty and can be used for model selection. We demonstrate the accuracy of our method on synthetic data and illustrate its application to recordings from primate cortex. We provide a customizable Python package implementing our framework with different generative models suitable for diverse applications.
| S-EPMC9835171 | biostudies-literature
Project description:Knowledge of the genetic relatedness between individuals is important in many research areas in quantitative genetics, conservation genetics, forensics, evolution, and ecology. In the absence of pedigree records, relatedness can be estimated from genetic marker data using a number of estimators. These estimators, however, make the critical assumption of a large random mating population without genetic structures. The assumption is frequently violated in the real world where geographic/social structures or nonrandom mating usually lead to genetic structures. In this study, I investigated two approaches to the estimation of relatedness between a pair of individuals from a subpopulation due to recent common ancestors (i.e., relatedness is defined and measured with the current focal subpopulation as reference). The indirect approach uses the allele frequencies of the entire population with and without accounting for the population structure, and the direct approach uses the allele frequencies of the current focal subpopulation. I found by simulations that currently widely applied relatedness estimators are upwardly biased under the indirect approach, but can be modified to become unbiased and more accurate by using Wright's F(st) to account for population structures. However, the modified unbiased estimators under the indirect approach are clearly inferior to the unmodified original estimators under the direct approach, even when small samples are used in estimating both allele frequencies and relatedness.
| S-EPMC3063680 | biostudies-other
Project description:Within-host genetic diversity and large transmission bottlenecks confound phylodynamic inference of epidemiological dynamics. Conventional phylodynamic approaches assume that nodes in a time-scaled pathogen phylogeny correspond closely to the time of transmission between hosts that are ancestral to the sample. However, when hosts harbor diverse pathogen populations, node times can substantially pre-date infection times. Imperfect bottlenecks can cause lineages sampled in different individuals to coalesce in unexpected patterns. To address realistic violations of standard phylodynamic assumptions we developed a new inference approach based on a multi-scale coalescent model, accounting for nonlinear epidemiological dynamics, heterogeneous sampling through time, non-negligible genetic diversity of pathogens within hosts, and imperfect transmission bottlenecks. We apply this method to HIV-1 and Ebola virus (EBOV) outbreak sequence data, illustrating how and when conventional phylodynamic inference may give misleading results. Within-host diversity of HIV-1 causes substantial upwards bias in the number of infected hosts using conventional coalescent models, but estimates using the multi-scale model have greater consistency with reported number of diagnoses through time. In contrast, we find that within-host diversity of EBOV has little influence on estimated numbers of infected hosts or reproduction numbers, and estimates are highly consistent with the reported number of diagnoses through time. The multi-scale coalescent also enables estimation of within-host effective population size using single sequences from a random sample of patients. We find within-host population genetic diversity of HIV-1 p17 to be 2Nμ=0.012 (95% CI 0.0066-0.023), which is lower than estimates based on HIV envelope serial sequencing of individual patients.
| S-EPMC5400386 | biostudies-literature
Project description:Population sizing from still aerial pictures is of wide applicability in ecological and social sciences. The problem is long standing because current automatic detection and counting algorithms are known to fail in most cases, and exhaustive manual counting is tedious, slow, difficult to verify and unfeasible for large populations. An alternative is to multiply population density with some reference area but, unfortunately, sampling details, handling of edge effects, etc., are seldom described. For the first time we address the problem using principles of geometric sampling. These principles are old and solid, but largely unknown outside the areas of three dimensional microscopy and stereology. Here we adapt them to estimate the size of any population of individuals lying on an essentially planar area, e.g. people, animals, trees on a savanna, etc. The proposed design is unbiased irrespective of population size, pattern, perspective artifacts, etc. The implementation is very simple-it is based on the random superimposition of coarse quadrat grids. Also, an objective error assessment is often lacking. For the latter purpose the quadrat counts are often assumed to be independent. We demonstrate that this approach can perform very poorly, and we propose (and check via Monte Carlo resampling) a new theoretical error prediction formula. As far as efficiency, counting about 50 (100) individuals in 20 quadrats, can yield relative standard errors of about 8% (5%) in typical cases. This fact effectively breaks the barrier hitherto imposed by the current lack of automatic face detection algorithms, because semiautomatic sampling and manual counting becomes an attractive option.
| S-EPMC4633052 | biostudies-literature
Project description:Accurate estimates of seasonal infection risk can be used by animal health officials to predict future disease risk and improve understanding of the mechanisms driving disease dynamics. It can be difficult to estimate seasonal infection risk in wildlife disease systems because surveillance assays typically target antibodies (serosurveillance), which are not necessarily indicative of current infection, and serosurveillance sampling is often opportunistic. Recently developed methods estimate past time of infection from serosurveillance data using quantitative serological assays that indicate the amount of antibodies in a serology sample. However, current methods do not account for common opportunistic and uneven sampling associated with serosurveillance data. We extended the framework of survival analysis to improve estimates of seasonal infection risk from serosurveillance data across population and regional scales. We found that accounting for the right-censored nature of quantitative serology samples greatly improved estimates of seasonal infection risk, even when sampling was uneven in time. Survival analysis can also be used to account for common challenges when estimating infection risk from serology data, such as biases induced by host demography and continually elevated antibodies following infection. The framework developed herein is widely applicable for estimating seasonal infection risk from serosurveillance data in humans, wildlife, and livestock.
| S-EPMC6940506 | biostudies-literature
Project description:Linkage disequilibrium (LD) is used to infer evolutionary history, to identify genomic regions under selection, and to dissect the relationship between genotype and phenotype. In each case, we require accurate estimates of LD statistics from sequencing data. Unphased data present a challenge because multilocus haplotypes cannot be inferred exactly. Widely used estimators for the common statistics r2 and D2 exhibit large and variable upward biases that complicate interpretation and comparison across cohorts. Here, we show how to find unbiased estimators for a wide range of two-locus statistics, including D2, for both single and multiple randomly mating populations. These unbiased statistics are particularly well suited to estimate effective population sizes from unlinked loci in small populations. We develop a simple inference pipeline and use it to refine estimates of recent effective population sizes of the threatened Channel Island Fox populations.
| S-EPMC7038669 | biostudies-literature
Project description:ObjectivesThe objective of this study is to extend the UNAIDS incidence estimation model, the UNAIDS Estimation and Projection Package (EPP), so that it can incorporate data from incidence assays.MethodsWe propose combining the likelihood of the incidence assay data with the likelihood of other data, in a manner that is consistent with the biomarker-based incidence estimator using incidence assay data. Two calibrating parameters specify the performance of the incidence assays: the false recent rate and the mean duration of recent infection. We then use synthetic data, based on prevalence data obtained from antenatal clinic surveillances, and in some cases household surveys, from 24 countries, to examine the impact of including incidence assay data, under circumstances wherein the incidence assay data imply the same or a different incidence rate as that inferred from the prevalence data alone, and wherein incorrect calibrating parameters for the incidence assay data are used.ResultsUsing incidence assay data, in addition to prevalence data, can improve estimate by narrowing uncertainty intervals in derived HIV incidence estimates, and by providing information on levels or trends in incidence that were not apparent in the prevalence data alone. However, the effect is relatively modest if the sample size of the incidence assay survey is small and results can be biased if the calibrating parameters for the incidence assay data are not known accurately.ConclusionIncorporating information from incidence assays in the manner proposed has the potential to improve estimates. Further work will examine in more detail the circumstances under which the contribution of incidence assay data would be most valuable.
| S-EPMC4247275 | biostudies-literature
Project description:In this article, we address the problem of measuring and analyzing sensation, the subjective magnitude of one's experience. We do this in the context of the method of triads: The sensation of the stimulus is evaluated via relative judgments of the following form: "Is stimulus \(S_i\) more similar to stimulus \(S_j\) or to stimulus \(S_k\)?" We propose to use ordinal embedding methods from machine learning to estimate the scaling function from the relative judgments. We review two relevant and well-known methods in psychophysics that are partially applicable in our setting: nonmetric multidimensional scaling (NMDS) and the method of maximum likelihood difference scaling (MLDS). Considering various scaling functions, we perform an extensive set of simulations to demonstrate the performance of the ordinal embedding methods. We show that in contrast to existing approaches, our ordinal embedding approach allows, first, to obtain reasonable scaling functions from comparatively few relative judgments and, second, to estimate multidimensional perceptual scales. In addition to the simulations, we analyze data from two real psychophysics experiments using ordinal embedding methods. Our results show that in the one-dimensional perceptual scale, our ordinal embedding approach works as well as MLDS, while in higher dimensions, only our ordinal embedding methods can produce a desirable scaling function. To make our methods widely accessible, we provide an R-implementation and general rules of thumb on how to use ordinal embedding in the context of psychophysics.
| S-EPMC7533746 | biostudies-literature
Project description:Accurate estimate of relatedness is important for genetic data analyses, such as heritability estimation and association mapping based on data collected from genome-wide association studies. Inaccurate relatedness estimates may lead to biased heritability estimations and spurious associations. Individual-level genotype data are often used to estimate kinship coefficient between individuals. The commonly used sample correlation-based genomic relationship matrix (scGRM) method estimates kinship coefficient by calculating the average sample correlation coefficient among all single nucleotide polymorphisms (SNPs), where the observed allele frequencies are used to calculate both the expectations and variances of genotypes. Although this method is widely used, a substantial proportion of estimated kinship coefficients are negative, which are difficult to interpret. In this paper, through mathematical derivation, we show that there indeed exists bias in the estimated kinship coefficient using the scGRM method when the observed allele frequencies are regarded as true frequencies. This leads to negative bias for the average estimate of kinship among all individuals, which explains the estimated negative kinship coefficients. Based on this observation, we propose an unbiased estimation method, UKin, which can reduce kinship estimation bias. We justify our improved method with rigorous mathematical proof. We have conducted simulations as well as two real data analyses to compare UKin with scGRM and three other kinship estimating methods: rGRM, tsGRM, and KING. Our results demonstrate that both bias and root mean square error in kinship coefficient estimation could be reduced by using UKin. We further investigated the performance of UKin, KING, and three GRM-based methods in calculating the SNP-based heritability, and show that UKin can improve estimation accuracy for heritability regardless of the scale of SNP panel.
| S-EPMC9727941 | biostudies-literature