Post-Flood Populations: Haplogroup formation and Fixation Dynamics in from Noah to Babel Dispersion
Description
From Noah to Babel: Genetics Meets Genesis What if the Flood, the Tower of Babel, and the 70 nations of Genesis weren’t just stories—but events that left a measurable imprint in our DNA? This groundbreaking study uses population modeling, genetic substitution rates, and fixation dynamics to show how today’s Y-chromosome and mtDNA haplogroups align with a post-Flood biblical timeline. With only six founders, exponential growth, a Babel bottleneck, and rapid dispersion, the genetic record converges with Scripture in ways mainstream science cannot explain away. The probability of this alignment happening by chance? Astronomically small—1 in a trillion. A bold fusion of history, linguistics, and genetics—testing the Bible against the data, and finding it holds up.
Files
Steps to reproduce
Methods & Workflow Study design. We tested whether observed Y-chromosome and mtDNA haplogroup structure can arise within a post-Flood, pre-/post-Babel window (~0–950 AF) under published growth and substitution rates. The workflow integrates historical constraints (Septuagint chronology for Peleg/Babel), demographic growth, uniparental mutation/substitution parameters, and stochastic population genetics. Historical & demographic inputs. Timeline: Flood = year 0; Peleg/Babel ≈ year 850 (LXX; Luke 3:36). Founders: 3 mating pairs split into 3 groups at year 0. Growth: 0.90%/yr, generation time = 25 yr → ~12,000 people by year 850. Structure/bottleneck: year 700–850, merge then subdivide to ~70 language subgroups, dispersing by year 850–950. Genetic parameters. Y-chromosome: 2–3 de novo mutations per generation (patrilineal). mtDNA: pedigree-based substitution rates spanning 1/17, 1/32, 1/43 per generation (control-region anchored), implemented as a discrete draw each generation to reflect study variance. Population-genetic modeling. Forward Wright–Fisher simulations for 70 isolated demes with small effective sizes (Ne_m, Ne_f tuned to founder-size ranges). State updates per generation: (i) draw uniparental substitutions from the rate set; (ii) apply a Y bias multiplier for effective fixation probability; (iii) drift/founder fixation within demes (small-Ne approximation). Outputs: (a) realized count of distinct haplogroups (Y, mtDNA), (b) tally of fixed substitutions, (c) percentile envelope over 1,000 Monte Carlo replicates. Backward (coalescent-style) check: sanity-check mtDNA divergence among L/M/N under 300 generations total, with a split at 5,323 ybp; compare expected pairwise differences vs. today’s ~24 fixed control-region substitutions. Sensitivity analyses. We varied: growth (0.4–0.9%/yr), generation time (22–28 yr), deme count (50–90), Ne ranges, and the per-generation rate mix; we recorded medians and 5th–95th percentiles for haplogroup counts and fixation totals. A 100-year post-Babel window (≈4 generations) was also isolated to estimate early polymorphism seeding across 70 demes. Software & reproducibility. All analyses were computational (no wet-lab reagents). Python 3.11 with numpy, scipy, pandas, random; configuration in config.yaml; fixed PRNG seeds per replicate set; results exported as CSV with run metadata (rates, Ne, seeds). Git versioning tracked scripts and configs; figure generation used matplotlib. Validation & reporting. We report means, medians, ranges, and 5th–95th percentiles across replicates; key checkpoints include (i) total population at year 850, (ii) realized mtDNA = 2 haplogroups, Y = ~4–6, and (iii) cumulative fixed substitutions per marker class. Cross-checks against published rate envelopes and present-day haplogroup counts provide external face validity.
Institutions
- University of Southern California