Dataset Information

Accounting for spatial sampling patterns in Bayesian phylogeography.

ABSTRACT: Statistical phylogeography provides useful tools to characterize and quantify the spread of organisms during the course of evolution. Analyzing georeferenced genetic data often relies on the assumption that samples are preferentially collected in densely populated areas of the habitat. Deviation from this assumption negatively impacts the inference of the spatial and demographic dynamics. This issue is pervasive in phylogeography. It affects analyses that approximate the habitat as a set of discrete demes as well as those that treat it as a continuum. The present study introduces a Bayesian modeling approach that explicitly accommodates for spatial sampling strategies. An original inference technique, based on recent advances in statistical computing, is then described that is most suited to modeling data where sequences are preferentially collected at certain locations, independently of the outcome of the evolutionary process. The analysis of georeferenced genetic sequences from the West Nile virus in North America along with simulated data shows how assumptions about spatial sampling may impact our understanding of the forces shaping biodiversity across time and space.

SUBMITTER: Guindon S

PROVIDER: S-EPMC8719894 | biostudies-literature | 2021 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Accounting for spatial sampling patterns in Bayesian phylogeography.

Guindon Stéphane S De Maio Nicola N

Proceedings of the National Academy of Sciences of the United States of America 20211201 52

Statistical phylogeography provides useful tools to characterize and quantify the spread of organisms during the course of evolution. Analyzing georeferenced genetic data often relies on the assumption that samples are preferentially collected in densely populated areas of the habitat. Deviation from this assumption negatively impacts the inference of the spatial and demographic dynamics. This issue is pervasive in phylogeography. It affects analyses that approximate the habitat as a set of disc ...[more]

PMID: 34930835

Similar Datasets

Project description:BackgroundData collected to inform time variations in natural population size are tainted by sampling error. Ignoring sampling error in population dynamics models induces bias in parameter estimators, e.g., density-dependence. In particular, when sampling errors are independent among populations, the classical estimator of the synchrony strength (zero-lag correlation) is biased downward. However, this bias is rarely taken into account in synchrony studies although it may lead to overemphasizing the role of intrinsic factors (e.g., dispersal) with respect to extrinsic factors (the Moran effect) in generating population synchrony as well as to underestimating the extinction risk of a metapopulation.Methodology/principal findingsThe aim of this paper was first to illustrate the extent of the bias that can be encountered in empirical studies when sampling error is neglected. Second, we presented a space-state modelling approach that explicitly accounts for sampling error when quantifying population synchrony. Third, we exemplify our approach with datasets for which sampling variance (i) has been previously estimated, and (ii) has to be jointly estimated with population synchrony. Finally, we compared our results to those of a standard approach neglecting sampling variance. We showed that ignoring sampling variance can mask a synchrony pattern whatever its true value and that the common practice of averaging few replicates of population size estimates poorly performed at decreasing the bias of the classical estimator of the synchrony strength.Conclusion/significanceThe state-space model used in this study provides a flexible way of accurately quantifying the strength of synchrony patterns from most population size data encountered in field studies, including over-dispersed count data. We provided a user-friendly R-program and a tutorial example to encourage further studies aiming at quantifying the strength of population synchrony to account for uncertainty in population size estimates.

Project description:A challenge to understanding biological diversification is accounting for community-scale processes that cause multiple, co-distributed lineages to co-speciate. Such processes predict non-independent, temporally clustered divergences across taxa. Approximate-likelihood Bayesian computation (ABC) approaches to inferring such patterns from comparative genetic data are very sensitive to prior assumptions and often biased toward estimating shared divergences. We introduce a full-likelihood Bayesian approach, ecoevolity, which takes full advantage of information in genomic data. By analytically integrating over gene trees, we are able to directly calculate the likelihood of the population history from genomic data, and efficiently sample the model-averaged posterior via Markov chain Monte Carlo algorithms. Using simulations, we find that the new method is much more accurate and precise at estimating the number and timing of divergence events across pairs of populations than existing approximate-likelihood methods. Our full Bayesian approach also requires several orders of magnitude less computational time than existing ABC approaches. We find that despite assuming unlinked characters (e.g., unlinked single-nucleotide polymorphisms), the new method performs better if this assumption is violated in order to retain the constant characters of whole linked loci. In fact, retaining constant characters allows the new method to robustly estimate the correct number of divergence events with high posterior probability in the face of character-acquisition biases, which commonly plague loci assembled from reduced-representation genomic libraries. We apply our method to genomic data from four pairs of insular populations of Gekko lizards from the Philippines that are not expected to have co-diverged. Despite all four pairs diverging very recently, our method strongly supports that they diverged independently, and these results are robust to very disparate prior assumptions.

Project description:The temporal variability of the thalamus in functional networks may provide valuable insights into the pathophysiology of schizophrenia. To address the complexity of the role of the thalamic nuclei in psychosis, we introduced micro-co-activation patterns (μCAPs) and employed this method on the human genetic model of schizophrenia 22q11.2 deletion syndrome (22q11.2DS). Participants underwent resting-state functional MRI and a data-driven iterative process resulting in the identification of six whole-brain μCAPs with specific activity patterns within the thalamus. Unlike conventional methods, μCAPs extract dynamic spatial patterns that reveal partially overlapping and non-mutually exclusive functional subparts. Thus, the μCAPs method detects finer foci of activity within the initial seed region, retaining valuable and clinically relevant temporal and spatial information. We found that a μCAP showing co-activation of the mediodorsal thalamus with brain-wide cortical regions was expressed significantly less frequently in patients with 22q11.2DS, and its occurrence negatively correlated with the severity of positive psychotic symptoms. Additionally, activity within the auditory-visual cortex and their respective geniculate nuclei was expressed in two different μCAPs. One of these auditory-visual μCAPs co-activated with salience areas, while the other co-activated with the default mode network (DMN). A significant shift of occurrence from the salience+visuo-auditory-thalamus to the DMN + visuo-auditory-thalamus μCAP was observed in patients with 22q11.2DS. Thus, our findings support existing research on the gatekeeping role of the thalamus for sensory information in the pathophysiology of psychosis and revisit the evidence of geniculate nuclei hyperconnectivity with the audio-visual cortex in 22q11.2DS in the context of dynamic functional connectivity, seen here as the specific hyper-occurrence of these circuits with the task-negative brain networks.

Dataset Information

Accounting for spatial sampling patterns in Bayesian phylogeography.

Publications

Accounting for spatial sampling patterns in Bayesian phylogeography.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets