Abstract:The development of analytical methods for Genome-wide Association Studies (GWAS) has outpaced the evolution of simulation techniques and pipelines. This disparity underscores the importance of innovative simulation methods that can keep pace with the rapidly increasing scale of GWAS. The median sample size of GWAS over the past ten years has exceeded 50,000 individuals, a trend that emphasizes the need for simulation tools capable of generating data on a similar or larger scale. This paper introduces a novel method, the small-group originating (SGO) model, utilizing the SLiM software for simulating individual-level GWAS data. Our standardized protocol facilitates the generation of tens of thousands of pseudo-individuals with millions of variants from small (30-90) open-access datasets. SGO stands out, especially when compared to the widely-used resampling method in HapGen, showcasing superior simulation efficiency for large sample sizes (> 13,000) of unrelated individuals. This capability is particularly relevant given the current trajectory towards larger GWAS, necessitating tools that can simulate datasets reflective of this growth. Additionally, SGO provides customization options and can model dynamic life cycles and mating across generations, positioning it as a highly promising alternative for GWAS simulations. In a case study, sensitivity analyses of chromosome-level principal component analysis and kinship coefficient estimation were conducted. The results highlighted the poor robustness of chromosome-level quality control (QC) indexes and the uneven distribution of population structure across chromosomes and ancestries, advocating for the caution against relying solely on chromosome-level QC statistics. With its flexible and efficient approach to generating pseudo GWAS data, our standardized SGO protocol emerges as a crucial asset for method development, power analysis, and benchmarking in GWAS research. It is especially vital in the context of accommodating the demands for large-scale simulations, aligning with the current and future scale of GWAS.

GPOPSIM: a simulation tool for whole-genome genetic data

The theory on and software simulating large-scale genomic data for genotype-by-environment interactions

G2P: a Genome-Wide-Association-Study simulation tool for genotype simulation, phenotype simulation and power evaluation.

Gips: A Software Guide to Sequencing-Based Direct Gene Cloning in Forward Genetics Studies

An Overview of Population Genetic Data Simulation

PGsim: A Comprehensive and Highly Customizable Personal Genome Simulator.

Gsdcreator: An Efficient And Comprehensive Simulator For Genarating Ngs Data With Population Genetic Information

SimHOEPI: A resampling simulator for generating single nucleotide polymorphism data with a high‐order epistasis model

EpiSIM: Simulation of Multiple Epistasis, Linkage Disequilibrium Patterns and Haplotype Blocks for Genome-Wide Interaction Analysis

GAPIT Version 2: an Enhanced Integrated Tool for Genomic Association and Prediction

SimBit: A high performance, flexible and easy-to-use population genetic simulator

Small-group originating model: Optimized individual-level GWAS simulation featured by SLiM and using open-access data

HAPSIMU: a Genetic Simulation Platform for Population-Based Association Studies

py_ped_sim - A flexible forward genetic simulator for complex family pedigree analysis

SPS: A Simulation Tool for Calculating Power of Set‐Based Genetic Association Tests

pgainsim: an R-package to assess the mode of inheritance for quantitative trait loci in GWAS

A review: simulation tools for genome-wide interaction studies

HPOSim: an R package for phenotypic similarity measure and enrichment analysis based on the human phenotype ontology

Simulating Realistic Genomic Data with Rare Variants

SSPOS: a Software for Simulating, Predicting and Optimizing Social Genetic Effects

Pysim-sv: a Package for Simulating Structural Variation Data with GC-biases