Essential Site Maintenance: Authorea-powered sites will be updated circa 15:00-17:00 Eastern on Tuesday 5 November.
There should be no interruption to normal services, but please contact us at [email protected] in case you face any issues.

loading page

gscramble: Simulation of admixed individuals without reuse of genetic material
  • +1
  • Eric Anderson,
  • Rachael Giglio,
  • Matthew DeSaix,
  • Timothy Smyser
Eric Anderson
Southwest Fisheries Science Center

Corresponding Author:[email protected]

Author Profile
Rachael Giglio
USDA Animal and Plant Health Inspection Service
Author Profile
Matthew DeSaix
Colorado State University
Author Profile
Timothy Smyser
USDA
Author Profile

Abstract

While a best practice for evaluating the behavior of genetic clustering algorithms on empirical data is to conduct parallel analyses on simulated data, these types of simulation techniques often involve sampling genetic data with replacement. In this paper we demonstrate that sampling with replacement, especially with large marker sets, inflates the perceived statistical power to correctly assign individuals (or the alleles that they carry) back to source populations—a phenomenon we refer to as resampling-induced, spurious power inflation (RISPI). To address this issue, we present gscramble a simulation approach in R for creating biologically informed individual genotypes from empirical data that: 1) samples alleles from populations without replacement, 2) segregates alleles based on species-specific recombination rates. This framework makes it possible to simulate admixed individuals in a way that respects the physical linkage between markers on the same chromosome and which does not suffer from RISPI. This is achieved in gscramble by allowing users to specify pedigrees of varying complexity in order to simulate admixed genotypes, segregating and tracking haplotype blocks from different source populations through those pedigrees, and then sampling—using a variety of permutation schemes—alleles from empirical data into those haplotype blocks. We demonstrate the functionality of gscramble with both simulated and empirical data sets and highlight additional uses of the package that users may find valuable.
15 Jul 2024Submitted to Molecular Ecology Resources
19 Jul 2024Submission Checks Completed
19 Jul 2024Assigned to Editor
19 Jul 2024Review(s) Completed, Editorial Evaluation Pending
14 Aug 2024Reviewer(s) Assigned