Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data
The pairwise sequentially Markovian coalescent (PSMC) method uses the genome sequence of a single individual to estimate demographic history covering a time span of thousands of generations. Although originally designed for whole-genome data, we here use simulations to investigate its applicability...
Main Authors: | , |
---|---|
Format: | Dataset |
Language: | unknown |
Published: |
2016
|
Subjects: | |
Online Access: | https://zenodo.org/record/5022512 https://doi.org/10.5061/dryad.0618v |
id |
ftzenodo:oai:zenodo.org:5022512 |
---|---|
record_format |
openpolar |
spelling |
ftzenodo:oai:zenodo.org:5022512 2023-05-15T16:30:26+02:00 Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data Liu, Shenglin Hansen, Michael M. 2016-10-05 https://zenodo.org/record/5022512 https://doi.org/10.5061/dryad.0618v unknown doi:10.1111/1755-0998.12606 https://zenodo.org/communities/dryad https://zenodo.org/record/5022512 https://doi.org/10.5061/dryad.0618v oai:zenodo.org:5022512 info:eu-repo/semantics/openAccess https://creativecommons.org/publicdomain/zero/1.0/legalcode Pairwise Sequential Markovian Coalescent (PSMC) whole genome sequencing threespine stickleback info:eu-repo/semantics/other dataset 2016 ftzenodo https://doi.org/10.5061/dryad.0618v10.1111/1755-0998.12606 2023-03-10T19:18:22Z The pairwise sequentially Markovian coalescent (PSMC) method uses the genome sequence of a single individual to estimate demographic history covering a time span of thousands of generations. Although originally designed for whole-genome data, we here use simulations to investigate its applicability to reference genome-aligned restriction site associated DNA (RAD) data. We find that RAD data can potentially be used for PSMC analysis, but at present with limitations. The key factor is the proportion (p) of the genome that the RAD data covers. In our simulations, a proportion of 10% can still retain a substantial amount of coalescent information, whereas for 1% estimation becomes unreliable. The performance depends strongly on mutation rate (μ) and recombination rate (r) and is proportional to μ*p/r. When the value of this term is low, increasing the amount of data and number of iterations helps restoring the power of the estimation. We subsequently analyse one whole-genome-sequenced and 17 RAD-sequenced three-spined sticklebacks (Gasterosteus aculeatus) from a lake in Greenland. The whole-genome sequence suggests a relatively recent expansion and decline within ca. 4000–40 000 generations ago, possibly reflecting postglacial expansion and founding of the lake population. RAD data, where chromosomes from 10 individuals are combined, identify a similar pattern. Our study provides guidance about the use of PSMC analysis and suggests measures that can improve its utility for RAD data. Finally, the study shows that RAD loci in general contain coalescent information that can be used for developing more targeted methods. empirical data(psmcfa files)The zip file contains the input files for PSMC from the empirical data of threespine sticklebacks in the study. Seventeen of them are from RAD data, and one of them is from whole genome sequencing data.Generating PSMC input files from simulated dataThe file contains a function written in R. This function randomly samples RAD loci from chromosomes simulated by ms (Hudson ... Dataset Greenland Zenodo Greenland Hudson |
institution |
Open Polar |
collection |
Zenodo |
op_collection_id |
ftzenodo |
language |
unknown |
topic |
Pairwise Sequential Markovian Coalescent (PSMC) whole genome sequencing threespine stickleback |
spellingShingle |
Pairwise Sequential Markovian Coalescent (PSMC) whole genome sequencing threespine stickleback Liu, Shenglin Hansen, Michael M. Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data |
topic_facet |
Pairwise Sequential Markovian Coalescent (PSMC) whole genome sequencing threespine stickleback |
description |
The pairwise sequentially Markovian coalescent (PSMC) method uses the genome sequence of a single individual to estimate demographic history covering a time span of thousands of generations. Although originally designed for whole-genome data, we here use simulations to investigate its applicability to reference genome-aligned restriction site associated DNA (RAD) data. We find that RAD data can potentially be used for PSMC analysis, but at present with limitations. The key factor is the proportion (p) of the genome that the RAD data covers. In our simulations, a proportion of 10% can still retain a substantial amount of coalescent information, whereas for 1% estimation becomes unreliable. The performance depends strongly on mutation rate (μ) and recombination rate (r) and is proportional to μ*p/r. When the value of this term is low, increasing the amount of data and number of iterations helps restoring the power of the estimation. We subsequently analyse one whole-genome-sequenced and 17 RAD-sequenced three-spined sticklebacks (Gasterosteus aculeatus) from a lake in Greenland. The whole-genome sequence suggests a relatively recent expansion and decline within ca. 4000–40 000 generations ago, possibly reflecting postglacial expansion and founding of the lake population. RAD data, where chromosomes from 10 individuals are combined, identify a similar pattern. Our study provides guidance about the use of PSMC analysis and suggests measures that can improve its utility for RAD data. Finally, the study shows that RAD loci in general contain coalescent information that can be used for developing more targeted methods. empirical data(psmcfa files)The zip file contains the input files for PSMC from the empirical data of threespine sticklebacks in the study. Seventeen of them are from RAD data, and one of them is from whole genome sequencing data.Generating PSMC input files from simulated dataThe file contains a function written in R. This function randomly samples RAD loci from chromosomes simulated by ms (Hudson ... |
format |
Dataset |
author |
Liu, Shenglin Hansen, Michael M. |
author_facet |
Liu, Shenglin Hansen, Michael M. |
author_sort |
Liu, Shenglin |
title |
Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data |
title_short |
Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data |
title_full |
Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data |
title_fullStr |
Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data |
title_full_unstemmed |
Data from: PSMC (pairwise sequentially Markovian coalescent) analysis of RAD (restriction site associated DNA) sequencing data |
title_sort |
data from: psmc (pairwise sequentially markovian coalescent) analysis of rad (restriction site associated dna) sequencing data |
publishDate |
2016 |
url |
https://zenodo.org/record/5022512 https://doi.org/10.5061/dryad.0618v |
geographic |
Greenland Hudson |
geographic_facet |
Greenland Hudson |
genre |
Greenland |
genre_facet |
Greenland |
op_relation |
doi:10.1111/1755-0998.12606 https://zenodo.org/communities/dryad https://zenodo.org/record/5022512 https://doi.org/10.5061/dryad.0618v oai:zenodo.org:5022512 |
op_rights |
info:eu-repo/semantics/openAccess https://creativecommons.org/publicdomain/zero/1.0/legalcode |
op_doi |
https://doi.org/10.5061/dryad.0618v10.1111/1755-0998.12606 |
_version_ |
1766020156909158400 |