Population genetics

1.1 Definition and scope

Population genetics is the branch of genetics that studies the distribution and change of genetic variation within populations over time. It applies mathematical and statistical models to understand how evolutionary forces—such as mutation, genetic drift, gene flow, natural selection, and non-random mating—alter allele and genotype frequencies. The field bridges Mendelian inheritance and Darwinian evolution, providing quantitative predictions for processes ranging from microevolution to speciation. Its scope includes measuring genetic diversity, quantifying population structure, inferring demographic history, and detecting signatures of selection.

1.2 Historical development

1.2.1 Early contributions: Hardy, Weinberg, Fisher, Haldane, Wright

The foundations of population genetics were laid in the early twentieth century. In 1908, Godfrey Hardy and Wilhelm Weinberg independently derived the principle of allele frequency equilibrium under ideal conditions, now known as the Hardy–Weinberg principle. Ronald A. Fisher, J. B. S. Haldane, and Sewall Wright then developed the mathematical theory of evolution by natural selection, mutation, and drift. Fisher’s *The Genetical Theory of Natural Selection* (1930), Haldane’s work on mutation–selection balance, and Wright’s concept of adaptive landscapes and F‑statistics established the field as a rigorous discipline.

1.2.2 Modern synthesis and neutral theory

During the 1930s–1940s, population genetics became a core component of the modern synthesis, unifying Mendelian genetics with Darwinian evolution and paleontology. Later, in the 1960s, Motoo Kimura proposed the neutral theory of molecular evolution, arguing that most observed genetic variation is selectively neutral and evolves primarily through mutation and genetic drift. This shift spurred the development of molecular population genetics and coalescent theory, and continues to influence contemporary research.

1.3 Basic concepts: alleles, genotypes, and allele frequencies

An allele is one of two or more alternative forms of a gene at a given locus. A genotype is the combination of alleles an individual carries at a locus. The allele frequency is the proportion of a specific allele among all alleles at that locus in a population, while the genotype frequency is the proportion of individuals carrying a particular genotype. These frequencies are the raw data for population genetic analyses, and changes in them over generations are the core subject of study.

2.1 Types of genetic variation

2.1.1 Single nucleotide polymorphisms (SNPs)

SNPs are the most common type of genetic variation, involving a single base-pair change in DNA sequence. They occur throughout the genome and are often used as markers in population genetic studies due to their abundance and ease of genotyping.

2.1.2 Insertions, deletions, and structural variants

Insertions and deletions (indels) add or remove one or more nucleotides. Larger structural variants include inversions, translocations, and complex rearrangements. These can affect gene function and population differentiation, though they are less frequent than SNPs.

2.1.3 Copy number variants and microsatellites

Copy number variants (CNVs) are duplications or deletions of DNA segments ranging from hundreds to millions of base pairs. Microsatellites (short tandem repeats) consist of repeated 1–6 bp motifs; their high mutation rate makes them useful for studying recent population history and kinship.

2.2 Measures of genetic diversity

2.2.1 Heterozygosity (observed and expected)

Observed heterozygosity is the proportion of individuals that are heterozygous at a locus. Expected heterozygosity, calculated under Hardy–Weinberg assumptions, is the probability that two randomly chosen alleles are different. Comparison of observed and expected values can reveal departures from random mating.

2.2.2 Nucleotide diversity (π)

Nucleotide diversity π is the average number of nucleotide differences per site between two randomly chosen DNA sequences in a population. It reflects both mutation rate and effective population size and is a common measure of within-population genetic variation.

2.2.3 Watterson's estimator (θ)

Watterson’s estimator θ is based on the number of segregating sites in a sample, adjusted for sample size. Under the neutral model with constant population size, it provides an estimate of the population mutation rate 4Nₑμ. Differences between π and θ can indicate population size changes or selection.

2.3 Molecular markers and genotyping methods

Molecular markers include SNPs, microsatellites, and restriction fragment length polymorphisms (RFLPs). Genotyping methods range from traditional gel-based assays (e.g., RFLP analysis) to high-throughput technologies such as DNA microarrays, genotyping‑by‑sequencing (GBS), and whole‑genome sequencing. These methods enable large-scale surveys of genetic variation across populations.

3.1 Assumptions of an ideal population

The Hardy–Weinberg principle describes a population in which allele and genotype frequencies remain constant across generations. The assumptions are: large population size (no genetic drift), random mating, no mutation, no gene flow, and no natural selection. Real populations rarely meet all conditions, but the principle serves as a null model against which evolutionary change can be tested.

3.2 Hardy–Weinberg equilibrium equation

3.2.1 Derivation for two alleles

For a locus with two alleles, A and a, with frequencies p and q (p + q = 1), random mating yields genotype frequencies AA = p², Aa = 2pq, and aa = q². These frequencies remain constant if the assumptions hold. The equation p² + 2pq + q² = 1 is the Hardy–Weinberg equilibrium.

3.2.2 Multiple alleles and polyploidy

For multiple alleles, the equilibrium genotype frequency of any homozygous genotype is the square of its allele frequency, and heterozygotes have frequency 2 times the product of the two allele frequencies. In polyploid organisms (e.g., tetraploids), higher-order binomial expansions apply, with genotype frequencies following the multinomial distribution.

3.3 Testing for deviations from equilibrium

Deviations are assessed by comparing observed genotype frequencies to those expected under Hardy–Weinberg using a chi‑square test or Fisher’s exact test. Significant deviations may indicate inbreeding, population substructure, selection, or genotyping errors. Many software packages include such tests.

3.4 Applications of Hardy–Weinberg

3.4.1 Estimation of allele frequencies

In cases where only phenotype data are available (e.g., recessive traits), allele frequencies can be estimated assuming Hardy–Weinberg equilibrium. For a recessive homozygote frequency of q², the recessive allele frequency q = √(q²) and the dominant allele frequency p = 1 – q.

3.4.2 Detection of selection or non-random mating

A consistent excess of homozygotes in a large, randomly mating population suggests inbreeding or population subdivision, while an excess of heterozygotes may indicate heterozygote advantage (overdominance) or assortative mating. Hardy–Weinberg tests thus provide a first filter for detecting evolutionary forces.

4.1 Mutation

4.1.1 Mutation rates and types

Mutation rates vary across loci and species; typical per‑base rates are ~10⁻⁸–10⁻⁹ per generation. Types include point mutations (transitions, transversions), insertions, deletions, and larger structural changes. Most mutations are neutral or deleterious, but a small fraction are beneficial.

4.1.2 Effect on allele frequencies

Mutation introduces new alleles at low frequencies. Over many generations, recurrent mutation can slowly change allele frequencies, but its effect per generation is very small compared to other forces unless population size is extremely large or selection is weak.

4.1.3 Mutation–selection balance

When a deleterious allele is recurrently introduced by mutation and removed by purifying selection, an equilibrium frequency is reached. The equilibrium allele frequency of a fully recessive deleterious allele (selection coefficient s) is approximately √(μ/s), where μ is the mutation rate.

4.2 Genetic drift

4.2.1 Random fluctuations in finite populations

Genetic drift is the change in allele frequencies due to random sampling of gametes each generation. In small populations, frequencies can fluctuate dramatically; eventually, an allele either becomes fixed (frequency = 1) or lost (frequency = 0). Drift reduces genetic diversity over time.

4.2.2 Effective population size (Ne)

Effective population size is the size of an ideal population that experiences the same rate of drift as the actual population. It is often smaller than the census size due to factors such as unequal sex ratios, variation in reproductive success, and population bottlenecks. Ne determines the strength of drift.

4.2.3 Founder effect and population bottlenecks

A founder effect occurs when a new population is established by a small number of individuals, leading to reduced genetic diversity and altered allele frequencies. A bottleneck is a temporary drastic reduction in population size, also causing loss of rare alleles and a deficit of heterozygosity relative to equilibrium.

4.2.4 Wright–Fisher model

The Wright–Fisher model is a discrete‑time model of drift: each generation, a finite population of N diploid individuals produces a large pool of gametes, from which the next generation is randomly sampled. Under this model, the variance in allele frequency change is p(1‑p)/(2N) and the probability of ultimate fixation for a neutral allele equals its initial frequency.

4.3 Gene flow (migration)

4.3.1 Island model and continent–island model

The island model (Wright) assumes a set of subpopulations exchanging migrants at a constant rate among all pairs. The continent–island model simplifies to one-way gene flow from a large continent into a smaller island. In both cases, migration reduces genetic differentiation and can counteract drift or selection.

4.3.2 Effects of gene flow on population differentiation

Gene flow homogenizes allele frequencies across populations, decreasing FST values. When migration is high, populations tend to remain similar; when migration is low or absent, differentiation increases by drift and local adaptation. Gene flow can also introduce advantageous alleles.

4.3.3 Admixture and clines

Admixture occurs when previously isolated populations interbreed, producing individuals with mixed ancestry. Clines are gradual changes in allele frequencies across geographic space, often resulting from a balance between gene flow and selection or from secondary contact after isolation.

4.4 Natural selection

4.4.1 Modes of selection: directional, stabilizing, disruptive

Directional selection favors one extreme phenotype, shifting allele frequencies toward that extreme. Stabilizing selection favors intermediate phenotypes, reducing variance and maintaining a stable mean. Disruptive selection favors both extremes, potentially leading to polymorphism or speciation.

4.4.2 Selection coefficients and fitness

Fitness is a measure of reproductive success. The selection coefficient s quantifies the relative disadvantage of a genotype: for a deleterious allele, fitness of the homozygote is 1‑s, and for a dominant beneficial allele, the advantage is often expressed as 1+s. Larger	s	leads to faster allele frequency change.

4.4.3 Heterozygote advantage (overdominance)

When heterozygotes have higher fitness than both homozygotes, both alleles can be maintained at stable equilibrium frequencies. Classic examples include the sickle‑cell allele in malaria‑endemic regions and the MHC locus in vertebrates.

4.4.4 Balancing selection and frequency-dependent selection

Balancing selection maintains multiple alleles in a population. Frequency‑dependent selection, a form of balancing selection, occurs when the fitness of a genotype depends on its frequency (e.g., rare advantage in host–parasite systems). This can sustain polymorphism indefinitely.

4.4.5 Selective sweeps and hitchhiking

A selective sweep occurs when a beneficial mutation rapidly increases in frequency, reducing variation at linked neutral sites by hitchhiking. The result is a local loss of diversity and a skewed allele frequency spectrum. Hard sweeps involve a single de novo mutation; soft sweeps may involve standing variation or multiple origins.

4.5 Non-random mating

4.5.1 Inbreeding and inbreeding depression

Inbreeding (mating between relatives) increases homozygosity across the genome. It can expose deleterious recessive alleles, reducing fitness—a phenomenon called inbreeding depression. The inbreeding coefficient F measures the probability that two alleles at a locus are identical by descent.

4.5.2 Assortative and disassortative mating

Assortative mating occurs when individuals mate with phenotypically similar partners (e.g., same height), increasing homozygosity at loci affecting that trait. Disassortative mating (mating with dissimilar partners) increases heterozygosity. These patterns affect genotype frequencies without altering allele frequencies.

4.5.3 Selfing and its genetic consequences

Self‑fertilization (selfing) is common in many plants. It rapidly reduces heterozygosity: after one generation of complete selfing, heterozygosity is halved. Over many generations, selfing leads to nearly complete homozygosity, shifting the population away from Hardy–Weinberg proportions and increasing the expression of recessive alleles.

5.1 F-statistics (FST, FIS, FIT)

5.1.1 Definition and interpretation

F‑statistics quantify the partitioning of genetic variation. FIS measures inbreeding within subpopulations (deficit of heterozygotes relative to Hardy–Weinberg), FST measures differentiation among subpopulations, and FIT measures overall inbreeding relative to the total population. FST ranges from 0 (no differentiation) to 1 (fixed differences).

5.1.2 Estimation from genetic data

FST is commonly estimated from allele frequencies using the formula Σ(πT‑πS)/πT, where πT is total diversity and πS is average within‑subpopulation diversity. Other estimators (e.g., Weir and Cockerham’s θ) correct for sample size and are implemented in many software packages.

5.2 Models of population structure

5.2.1 Island model

In Wright’s infinite island model, a large number of subpopulations exchange migrants at a constant rate m each generation. Under equilibrium, FST ≈ 1/(1+4Nₑm). This model provides a simple prediction of the relationship between gene flow and differentiation.

5.2.2 Stepping-stone model

The stepping‑stone model arranges populations in a linear or grid pattern, with gene flow occurring only between neighboring demes. This yields a pattern of isolation by distance, where differentiation increases with geographic distance.

5.2.3 Isolation-by-distance

Isolation‑by‑distance describes the positive correlation between genetic differentiation (e.g., FST) and geographic distance. Under a stepping‑stone or continuous habitat model, populations exchange fewer genes with distant neighbors, leading to gradually increasing differentiation.

5.3 Clustering and principal component analysis

Clustering algorithms (e.g., STRUCTURE, ADMIXTURE) assign individuals to inferred genetic clusters based on multilocus genotype data, often assuming Hardy–Weinberg within clusters. Principal component analysis (PCA) reduces the genetic covariance matrix to a few axes that summarize population structure, revealing gradients or discrete groups.

5.4 Gene flow and barriers to migration

Gene flow can be limited by physical barriers (mountains, oceans), ecological differences (habitat preferences), or behavioral factors (mating preferences). Identifying barriers is important for understanding speciation and conservation. Methods based on coalescent theory can estimate historical migration rates between populations.

6.1 Basic coalescent process

6.1.1 Kingman's coalescent

Kingman’s coalescent models the ancestry of a sample of genes backward in time. In a constant‑size population, the waiting time until two lineages coalesce (merge into a common ancestor) is exponentially distributed with rate 1/(2Nₑ). The genealogy of a sample is a random tree with coalescent intervals that increase as fewer lineages remain.

6.1.2 Time to most recent common ancestor (TMRCA)

The expected time to the most recent common ancestor (TMRCA) of a sample of n genes is 2Nₑ(1‑1/n) generations. The actual TMRCA varies stochastically. Coalescent theory provides the distribution of TMRCA, which can be used to infer effective population size and demographic events.

6.2 Coalescent with recombination

Recombination breaks the genome into segments with different genealogies. The ancestral recombination graph (ARG) extends the coalescent by allowing lineages to recombine, creating multiple correlated trees along the chromosome. This is essential for understanding linkage disequilibrium and mapping.

6.3 Coalescent with population size changes

When population size changes (e.g., bottlenecks, expansions), coalescent rates vary accordingly. During a bottleneck, coalescent times are shortened; during expansion, coalescent intervals lengthen. Demographic inference uses the distribution of coalescent times from genetic data to estimate such changes.

6.4 Applications of coalescent theory

6.4.1 Demographic inference

Coalescent methods can estimate past population sizes, growth rates, divergence times, and migration rates. Programs like ∂a∂i (diffusion approximation for demographic inference) and MSMC (multiple sequentially Markovian coalescent) apply coalescent models to site‑frequency spectra or whole‑genome sequences.

6.4.2 Detection of natural selection

Coalescent simulations under a neutral model provide null distributions for test statistics. Deviations in the site frequency spectrum or in linkage disequilibrium patterns can indicate selection, as selection alters the shape of the genealogy (e.g., leading to an excess of rare variants under a recent sweep).

7.1 Neutral theory and nearly neutral theory

7.1.1 Kimura's neutral theory

Motoo Kimura’s neutral theory posits that most molecular evolution and within‑species variation are caused by mutation and genetic drift, not natural selection. The rate of neutral substitution equals the neutral mutation rate, independent of population size. This theory explains the observed patterns of molecular clocks and the high level of polymorphism.

7.1.2 Molecular clock and rate of evolution

The molecular clock is the roughly constant rate at which neutral substitutions accumulate over time. For proteins, the rate is often measured as substitutions per site per year. The nearly neutral theory (Ohta) extends this by considering slightly deleterious mutations that behave as neutral in small populations but are selected against in large ones.

7.2 Tests of neutrality

7.2.1 Tajima's D

Tajima’s D compares two estimators of θ: π (nucleotide diversity) and Watterson’s estimator (based on segregating sites). Under neutrality and constant population size, D ≈ 0. Negative D indicates an excess of rare variants (e.g., population expansion or purifying selection); positive D suggests balanced polymorphism or population bottleneck.

7.2.2 Fu and Li's tests

Fu and Li’s tests (D*, F*) use the number of singletons (mutations appearing only once in a sample) to distinguish selection from demography. These tests are sensitive to recent sweeps (excess of singletons) and to population growth (singleton excess) but are often used together with Tajima’s D.

7.2.3 McDonald–Kreitman test

This test compares the ratio of non‑synonymous to synonymous fixed differences between species to the same ratio within species. A significant excess of non‑synonymous fixed differences indicates positive selection acting across species; an excess of non‑synonymous polymorphism suggests balancing selection or slightly deleterious mutations.

7.3 Patterns of polymorphism and divergence

7.3.1 Site frequency spectrum

The site frequency spectrum (SFS) is the distribution of allele frequencies of polymorphisms in a sample. The SFS reflects demographic history and selection: population expansion yields an excess of low‑frequency variants; a selective sweep creates a local departure from the neutral SFS.

7.3.2 Linkage disequilibrium and haplotype structure

Linkage disequilibrium (LD) is the non‑random association of alleles at different loci. LD decays with recombination and is influenced by selection, admixture, and population bottlenecks. The extent of LD (e.g., r², D′) is used for mapping traits and inferring historical recombination rates.

8.1 Detecting selection from genomic data

8.1.1 FST outlier tests

Outlier analysis identifies loci with FST values significantly higher than the neutral background. Such loci are candidates for local adaptation (divergent selection). Methods may use a genome‑wide empirical distribution or a model based on island model expectations.

8.1.2 Extended haplotype homozygosity (EHH, iHS)

EHH measures the decay of LD around a core SNP. A long‑range, high‑frequency haplotype suggests a recent selective sweep (i.e., the beneficial allele rose to high frequency faster than recombination could break haplotypes). The integrated haplotype score (iHS) compares these signals across the genome.

8.1.3 dN/dS ratios

Within coding sequences, the ratio of non‑synonymous to synonymous substitution rates (dN/dS) indicates selection: dN/dS < 1 purifying selection, = 1 neutrality, > 1 positive selection. This is typically applied between species but can also be used to detect intra‑species selective constraint.

8.2 Adaptive evolution and selective sweeps

8.2.1 Hard and soft sweeps

A hard sweep results from a single beneficial mutation that arises and becomes fixed, eliminating linked variation. A soft sweep occurs when multiple copies of a beneficial allele are present (either from standing variation or recurrent mutation). Soft sweeps leave weaker signatures of reduced diversity.

8.2.2 Polygenic adaptation

Many traits are influenced by many loci. Polygenic adaptation refers to small allele‑frequency shifts across many loci, collectively producing a phenotypic change. This can be detected by methods that combine association statistics across the genome (e.g., polygenic risk score comparisons).

8.3 Local adaptation and clinal variation

Local adaptation occurs when populations evolve different phenotypes in response to different environments. Clinal variation—smooth changes in allele frequency along an environmental gradient—often reflects a balance between selection and gene flow. Methods like BayEnv and LFMM identify allele–environment correlations.

9.1 Linkage disequilibrium and association mapping

9.1.1 Genome-wide association studies (GWAS)

GWAS test millions of SNPs across the genome for statistical associations with a trait of interest. Significant SNPs often lie in or near genes influencing the trait. Population structure must be accounted for to avoid spurious associations. The effect sizes of associated variants are typically small.

9.1.2 Fine-mapping causal variants

Once a genomic region is associated, fine‑mapping uses dense genotyping, imputation, and statistical methods (e.g., Bayesian credible sets) to narrow the set of candidate causal variants. Information from functional genomics (e.g., regulatory elements) helps pinpoint the true causal site.

9.2 Heritability and its estimation

9.2.1 Broad-sense and narrow-sense heritability

Broad‑sense heritability (H²) is the proportion of phenotypic variance due to all genetic effects (additive, dominance, epistasis). Narrow‑sense heritability (h²) is the proportion due to additive genetic effects alone. h² determines the response to selection and is estimated from twin studies, pedigrees, or genomic relationships (GREML).

9.2.2 Partitioning heritability

Genomic partitioning attributes heritability to functional categories (e.g., coding regions, enhancers) or to chromosomes. This reveals whether variants in certain classes contribute disproportionately to trait variation. Methods use mixed linear models with multiple variance components.

9.3 Evolutionary quantitative genetics

9.3.1 Breeder's equation and response to selection

The breeder’s equation, R = h²S, predicts the response to selection R from the selection differential S (difference between selected parents and population mean) and narrow‑sense heritability. It is foundational for artificial selection in agriculture and for understanding natural selection on quantitative traits.

9.3.2 G matrix and multivariate evolution

The G matrix summarizes additive genetic variances and covariances among multiple traits. Multivariate selection changes the mean of a vector of traits, and the G matrix determines the evolutionary trajectory. Constraints imposed by G (e.g., genetic correlations) can limit or channel adaptation.

10.1 Conservation genetics

10.1.1 Identifying evolutionarily significant units (ESUs)

ESUs are populations or groups that merit separate conservation priority. Population genetic data (e.g., FST, phylogenies) help delineate ESUs based on historical isolation or adaptive distinctiveness.

10.1.2 Managing genetic diversity in threatened species

Low genetic diversity increases extinction risk due to inbreeding depression and reduced adaptability. Conservation programs monitor heterozygosity, allelic richness, and effective population size to guide captive breeding and habitat management.

10.1.3 Inbreeding avoidance in captive breeding

Pedigree analyses and molecular markers identify related individuals to minimize inbreeding in captive populations. Genetic management strategies include equalizing family sizes and swapping individuals among breeding facilities.

10.2 Human population genetics

10.2.1 Out-of-Africa migration and global expansion

Genetic evidence supports a single origin of modern humans in Africa, followed by migrations into Eurasia and the rest of the world. Patterns of diversity (e.g., declining heterozygosity with distance from Africa) and high‑coverage ancient genomes refine the timings and routes of these expansions.

10.2.2 Admixture and population structure in humans

Human populations show varying degrees of admixture (e.g., in the Americas, South Asia). Population structure (e.g., in Europe) reflects geography, language, and historical events. These patterns are important for correcting stratification in medical studies.

10.2.3 Genetic epidemiology and disease mapping

Population genetics informs the design of case‑control studies for complex diseases. For example, mapping by admixture linkage disequilibrium (MALD) uses differences in allele frequencies between ancestral populations to localize disease genes.

10.3 Agricultural and domestication genetics

10.3.1 Genetic improvement of crops and livestock

Marker‑assisted selection and genomic selection use population genetic principles to accelerate breeding. SNP arrays and genotyping-by-sequencing enable prediction of breeding values for traits like yield and disease resistance.

10.3.2 Domestication signatures in genomes

Comparison of domestic vs. wild genomes reveals loci under selection during domestication (e.g., the *sh1* gene in maize, the *MC1R* gene in dog coat color). Such signatures often exhibit reduced diversity and elevated FST relative to background.

10.4 Forensic genetics

10.4.1 DNA fingerprinting and ancestry inference

Short tandem repeats (STRs) are used for individual identification. Ancestry‑informative markers (AIMs) from population genetics can estimate an individual’s biogeographical origin, assisting law enforcement.

10.4.2 Population databases and match probabilities

Forensic databases (e.g., CODIS) contain allele frequencies from reference populations. Match probabilities are calculated using the product rule or corrections for population substructure (e.g., θ‑correction) to avoid overestimating the rarity of a match.

11.1 Simulation tools (e.g., ms, SLiM)

Forward‑time simulators like SLiM and coalescent simulators like ms allow researchers to generate genetic data under specified demographic and selection models. These simulations are used to test statistical methods, estimate power, and compare observed data to null distributions.

11.2 Statistical inference software (e.g., PLINK, ADMIXTURE, STRUCTURE)

PLINK provides tools for GWAS, LD calculations, and basic population genetics. ADMIXTURE and STRUCTURE estimate individual ancestry proportions under a model of K ancestral populations. Other software like BEAST and MCMCcoal estimate phylogenies and demographic parameters using Bayesian inference.

11.3 Likelihood and Bayesian approaches

Maximum‑likelihood methods (e.g., in ∂a∂i, fastsimcoal2) fit demographic models to the SFS or other summary statistics. Bayesian approaches (e.g., approximate Bayesian computation, ABC) compare simulated data to observed data via distance metrics, useful for complex models where likelihoods are intractable.

11.4 Machine learning in population genetics

Supervised learning (e.g., random forests, neural networks) can classify populations, detect selection, and infer demographic parameters from genomic data. Methods like mix‐models and autoencoders are increasingly applied to large genetic datasets to capture nonlinear relationships.

12.1 Ancient DNA and temporal population genetics

Advances in sequencing ancient DNA allow direct observation of allele frequency changes over thousands of years. Temporal data provide opportunities to estimate selection coefficients and migration rates with higher precision and to uncover population turnover events (e.g., the Neolithic transition in Europe).

12.2 Population genomics in non-model organisms

Cheaper sequencing technologies have expanded population genetics beyond traditional model organisms. Studies of natural populations (e.g., wild fish, trees, insects) now routinely generate genome‑wide data, enabling tests of local adaptation, speciation genomics, and landscape genetics.

12.3 Integrative models of demography and selection

New methods jointly infer demography, migration, and selection in a unified framework. For example, forward simulations with SLiM can incorporate realistic genomic features and environments, while approximate Bayesian computation can estimate complex scenarios from summary statistics.

12.4 Ethical considerations in population genetic research

Population genetics raises ethical issues regarding consent, privacy, and the potential misuse of ancestry or disease‑risk information. Researchers must ensure transparent communication of results, avoid reinforcing racial stereotypes, and respect the rights of indigenous communities whose genetic samples are studied. Standards for data sharing and governance continue to evolve.