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.

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

GPOPSIM: a simulation tool for whole-genome genetic data

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

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

Mixed Linear Model Approaches of Association Mapping for Complex Traits Based on Omics Variants

Methods for predicting superior genotypes under multiple environments based on QTL effects

Mapping QTLs with digenic epistasis under multiple environments and predicting heterosis based on QTL effects

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

A review: simulation tools for genome-wide interaction studies

Mapping the genetic architecture of complex traits in experimental populations.

GbyE: an integrated tool for genome widely association study and genome selection based on genetic by environmental interaction

An Overview of Population Genetic Data Simulation

Enhancing Genomic Selection by Fitting Large-Effect SNPs As Fixed Effects and a Genotype-by-environment Effect Using a Maize BC1F3:4 Population.

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

[Analytic Methods for Seed Models with Genotype X Environment Interactions].

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

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

G × EBLUP: A Novel Method for Exploring Genotype by Environment Interactions and Genomic Prediction.

GEM: scalable and flexible gene–environment interaction analysis in millions of samples

GAPIT: genome association and prediction integrated tool.

G2P Provides an Integrative Environment for Multi-model Genomic Selection Analysis to Improve Genotype-to-phenotype Prediction.