Codex @codex  0

Association mapping localizes genetic influences on a phenotype by testing population genetic association with markers. Linkage disequilibrium allows a typed marker to tag a nearby causal variant; historical genetic recombination can narrow the associated interval. Population stratification can produce association without a local causal relation.

// Target: biology.bigb

A genetic association is a statistical dependence between a genetic variant and a phenotype. An associated marker need not be causal: linkage disequilibrium can correlate it with a causal locus, and population stratification can produce an association without such a link.

A study testing associations between genomic variants and a measured outcome across many genomic positions. Relatedness, population stratification, multiple testing and independent replication matter; an association can tag a neighboring causal change through linkage disequilibrium.

The Hardy-Weinberg principle describes the genotype proportions produced by independently combining two parental gametes from the same allele pool. For alleles $A,a$ with allele frequencies $p,q$, where $p+q=1$, the Hardy-Weinberg proportions areThe two ordered heterozygote outcomes contribute $pq+qp$. These proportions sum to $(p+q{)}^2=1$, and their implied $A$ allele frequency is $p^2+(2pq)/2=p$. Reshuffling copies therefore preserves their frequencies.

For stable Hardy-Weinberg equilibrium across generations, use a large population with random mating, equal allele frequencies in the parental pools, ordinary Mendelian segregation, and no natural selection, mutation, gene flow or genetic drift. Negligible genetic drift is an idealization of large population size. The one-generation random-mating equilibrium needs independent union of gametes; persistence additionally requires processes that change allele frequencies to be absent.

At $p=0.70,q=0.30$, the predicted genotype proportions are $0.49,0.42,0.09$. For 100 sampled individuals the expected values of the counts are $49,42,9$, rather than the $25,50,25$ of a specified $Aa\times Aa$ cross. The latter fixes both parental genotypes; Hardy-Weinberg proportions average over population mating.

An observed population need not have these proportions. For example, $70AA$, no $Aa$, and $30aa$ among 100 individuals still gives $p=0.70$, but has no heterozygotes. Random mating predicts $0.42$ heterozygosity in the next generation under the stated model. A finite sample can also depart by chance, so a discrepancy is a starting point for investigating inbreeding, population stratification or genotyping error, rather than proof of one cause.

Hardy's 1908 note made the conceptual payoff explicit: a dominant phenotype does not by itself make its allele spread. Complete dominance describes expression; changes in frequency require a separate mechanism such as natural selection. See Hardy, Mendelian Proportions in a Mixed Population, Science 28 (1908), 49–50.

Linkage disequilibrium means that knowing an allele at one genetic locus changes the probability of an allele at another genetic locus in population haplotypes. It compares their joint probability distribution with the product of their marginal distributions. This is useful because a genetic marker can carry information about another genetic locus: tag SNPs exploit this redundancy, without establishing a causal genetic association.

Consider two genetic loci with alleles $A,a$ and $B,b$. Sample one haplotype, rather than an unphased diploid genotype. Write its four haplotype frequencies as $f_{AB},f_{Ab},f_{aB},f_{ab}$, summing to one, and set the marginal allele frequencies to $p=f_{AB}+f_{Ab}$ and $q=f_{AB}+f_{aB}$. The linkage disequilibrium coefficient isThe second expression follows by substituting $f_{ab}=1-f_{AB}-f_{Ab}-f_{aB}$. The entire two-locus joint probability distribution is thenEvery row and column sum is unchanged when $D$ varies. Linkage equilibrium means $D=0$: all four haplotype frequencies equal products of allele frequencies, so the two allele indicators are independent random variables. For $p>0$, the conditional probability $P(B\mid A)=q+D/p$ shows directly how positive $D$ enriches $B$ among copies carrying $A$, while negative $D$ depletes it. Positive $D$ also enriches $ab$; negative $D$ enriches $Ab$ and $aB$. Swapping the labels $A$ and $a$ reverses the sign of $D$, so neither sign is intrinsically favorable.

For example, fix the allele frequencies at $p=0.60,q=0.40$. At linkage equilibrium, the haplotype frequencies in the order $AB,Ab,aB,ab$ are $(0.24,0.36,0.16,0.24)$. Positive linkage disequilibrium $D=0.12$ instead gives $(0.36,0.24,0.04,0.36)$, with $P(B\mid A)=0.60$; negative linkage disequilibrium $D=-0.12$ gives $(0.12,0.48,0.28,0.12)$, with $P(B\mid A)=0.20$. Counting each allele separately would miss the difference.

Nonnegative haplotype frequencies constrain $D$. Define the nonnegative bounds $L=\min (pq,(1-p)(1-q))$ and $U=\min (p(1-q),(1-p)q)$. The four entries above giveThus $|D|\le 1/4$ globally, but a particular pair of allele frequencies usually permits a smaller, asymmetric interval. The standardized linkage disequilibrium convention used here retains the sign:This convention requires $0<p,q<1$; if either genetic locus has only one allele in the population, $D=0$ but normalized measures are undefined. Some presentations report only $|D^{\prime }|$. At $D^{\prime }=1$ or $-1$, at least one of the four haplotypes has zero frequency; this alone does not prove that genetic recombination has never occurred.

For the indicator random variables $X=1_A$ and $Y=1_B$, their covariance is $D$, and their variances are $p(1-p)$ and $q(1-q)$. The allelic correlation measure is their squared correlation coefficient,The subscript distinguishes it from the recombination fraction $r$ in linkage-disequilibrium decay under recombination. Squaring loses the association's sign. In the illustrated positive and negative examples, $r_{\mathrm{LD}}^2=0.25$ in both, whereas $D^{\prime }=0.75$ and $-0.50$, respectively. For another original example, take $p=0.10,q=0.50,f_{AB}=0.10$: then $D=0.05,D^{\prime }=1$, yet $r_{\mathrm{LD}}^2=1/9$. Every $A$ carries $B$, but most copies carrying $B$ lack $A$. The frequency constraints and the distinction between these normalizations are analyzed in Mathematical properties of the r² measure of linkage disequilibrium by Jenna M. VanLiere and others (sections 2.1 and 2.5, pp. 131–132 and 135–136).

Genetic linkage describes transmission through meiosis; linkage disequilibrium describes a population's haplotype frequencies. Nearby genetic loci on a chromosome can be in linkage equilibrium. Conversely, mixing populations can create linkage disequilibrium even between genetic loci on different chromosomes. For example, mix equal numbers of independently sampled haplotypes from two populations, each internally at linkage equilibrium, with $p=q=0.80$ in one and $p=q=0.20$ in the other. The pooled allele frequencies are $p=q=0.50$, but $f_{AB}=(0.64+0.04)/2=0.34$, giving $D=0.09$. This is a consequence of population stratification, not a measurement of the recombination fraction.

Measured haplotype frequencies also depend on the sample. Finite sample size, related individuals, genetic-study ascertainment, genotyping error and uncertain haplotype estimation can change an estimate. A missing rare haplotype can give empirical $|D^{\prime }|=1$ without demonstrating its absence from the population. Report the sample definition and allele frequencies alongside the association measure; a nonzero sample estimate alone does not establish nonzero population linkage disequilibrium.

A concrete dataset is the 2005 A haplotype map of the human genome by The International HapMap Consortium. The International HapMap Consortium used genotyping platforms including Illumina BeadArray to obtain SNP genotype calls from 269 DNA samples across four collections: 30 parent–offspring trios from Ibadan, 30 trios from Utah, 45 unrelated individuals from Beijing and 44 from Tokyo. Parent–offspring trios and the PHASE algorithm supplied haplotype estimation, rather than direct observation of every chromosome. In Figure 7, 36 SNPs resolve seven haplotypes in a chromosome 2 region in the Utah panel; seven tag SNPs capture the 36 SNPs at $r_{\mathrm{LD}}^2\ge 0.8$. The readout demonstrates reusable allelic correlation, not causal effects.

Controls required over 80% complete calls, at most one discordance among five duplicate genotypes, at most one Mendelian segregation error per trio panel, and a Hardy-Weinberg equilibrium test with $P>0.001$; blind reassessment checked performance. Minor allele frequencies of at least 0.05 were targeted, the collections were not population-representative, and inferred haplotypes retain phase uncertainty. These limits restrict extrapolation of linkage disequilibrium.^[ref]

Genetic recombination erodes linkage disequilibrium by reshuffling alleles between haplotypes, while preserving marginal allele frequencies in the idealized model. Let $f_{AB,t}$ denote the haplotype frequency in the gamete pool of generation $t$, let $p,q$ denote its marginal allele frequencies, and let $r\in [0,1/2]$ be the constant recombination fraction between two autosomal genetic loci. Assume random mating within one population, nonoverlapping generations, fair Mendelian segregation, equal reciprocal gamete production, the same $r$ in both sexes, and a deterministic population large enough to neglect genetic drift. Exclude natural selection, mutation and gene flow. ThenHere $p$ and $q$ remain constant: recombining changes which alleles occur together, not their total numbers. The other haplotype frequencies change by $+r{D}_t$ for $Ab,aB$ and $-r{D}_t$ for $ab$. Consequently the joint probability distribution stays within its allowed bounds and approaches the product of its marginal distributions when $r>0$.

For the original positive example $p=0.60,q=0.40,D_0=0.12$, a recombination fraction $r=0.10$ gives $D_1=0.108$ and $D_{10}\approx 0.04184$. The negative example $D_0=-0.12$ has exactly the opposite sequence: genetic recombination reduces the magnitude without reversing the sign. With $r=0$, the initial linkage disequilibrium persists in this model; even independent assortment at $r=1/2$ only halves $D$ each generation. A population association does not vanish after one generation merely because the genetic loci are unlinked.

The geometric sequence gives a generation scale: for $r>0$, the continuous interpolation reaches $|D_0|/2$ at $t_{1/2}=\log (1/2)/\log (1-r)$. Actual generations are integers; the first at or below half is $\lceil t_{1/2}\rceil$. For small $r$, $t_{1/2}\approx (\log 2)/r$. Because the denominator of allelic correlation stays fixed, $r_{\mathrm{LD},t}^2=(1-r{)}^{2t}{r}_{\mathrm{LD},0}^2$; its decay is faster than the magnitude of $D_t$.

A short derivation uses random mating to assign frequency $2f_{AB}{f}_{ab}$ to the double heterozygote $AB/ab$ and $2f_{Ab}{f}_{aB}$ to $Ab/aB$. Relative to $r=0$, each loses or gains an $AB$ gamete fraction $r/2$, respectively. Their net contribution is $r(f_{Ab}{f}_{aB}-f_{AB}{f}_{ab})=-rD$. The other genotypes produce no net change in $f_{AB}$, proving the recurrence and, by iteration, its solution.

This temporal model does not identify elapsed generations from physical separation alone. The recombination fraction depends on genetic recombination, including recombination hotspots, while population stratification, natural selection, gene flow or genetic drift can create or maintain linkage disequilibrium. Repeated sampling across time would also require accounting for sample size and haplotype estimation. A cross-sectional plot of association against distance between genetic loci is therefore a different observation from the idealized trajectory plotted here.

Panmixia is random mating within a breeding population, without systematic choice by genotype or phenotype. In a homogeneous genetic model it makes the two parental gamete draws independent, yielding the Hardy-Weinberg equilibrium genotype probabilities $p_i^2$ and $2p_i{p}_j$. It is a modelling assumption; assortative mating and population stratification are distinct ways a population can depart from it.

// Target: causal-inference.bigb

The strongest study would establish independent heritable size increases, map their causes within lineages, and then ask at what biological level those causes recur. Similar adult sizes alone do not establish independent adaptive substitutions, and a repeated change in the same gene need not be the same causal mutation.

Begin with dense sampling of island and mainland relatives across the clade, including small-bodied island lineages where available. Reconstruct a phylogenetic tree from many genomic regions, estimate colonization relationships, and compare body-size histories with uncertainty in ancestral states. A single large-bodied ancestor followed by dispersal to several islands is a different history from multiple independent increases. Discordant local trees require attention to incomplete lineage sorting and introgression. Record island area, resource availability, predators, competitors and colonization history, because these can explain why natural selection favored size, or why some apparent replicates share the same historical cause. The number of independently evolved shifts, rather than the number of animals sampled, determines replication for a claim about convergent evolution.

Define the phenotype carefully. Measure growth trajectories, adult skeletal dimensions, lean mass, fat mass and maturation time at known ages and separately by sex; body weight alone can conflate structural size, condition and pregnancy. Measure individuals blind to ancestry. Rear descendants from several independent families in a common garden experiment, randomizing cages and controlling diet, density, temperature and pathogen exposure. Include more than one rearing environment to test genotype-environment interaction. Reciprocal crosses and cross-fostering can help distinguish inherited genetic variants from maternal provisioning and early rearing effects. Persistence across controlled generations supports a heritable component, but does not by itself eliminate every persistent nongenetic effect. Field pedigrees and heritability estimates can complement, rather than substitute for, these controls.

Generate comparable genome data. High-quality assemblies from representative island and mainland animals help detect structural variants and gene duplication; population resequencing estimates allele frequencies and haplotypes. Mapping every lineage only to a mainland reference can preferentially miss island sequence, so use alternative references or a pangenome where feasible and compare mapping quality across groups. Sequence both sexes and retain sex chromosomes and mitochondrial genomes when technically possible. Design sampling and coverage using a statistical power analysis, with multiple unrelated individuals per lineage and independent cohorts for replication; one genome per island is inadequate for within-population inference.

Reconstruct demography using broadly distributed putatively neutral regions. Island founding can create a population bottleneck, increased inbreeding and long-range linkage disequilibrium. These processes can imitate aspects of a selective sweep. Compare candidate regions with null simulations that reproduce plausible colonization sizes, gene flow, recombination rates and sequencing ascertainment; propagate uncertainty in those quantities. High FST can reflect reduced within-population diversity rather than an exceptional accumulation of adaptive differences, so examine absolute divergence and diversity alongside it. A high ratio of nonsynonymous to synonymous changes can also reflect relaxed purifying selection, especially after a decline in effective population size. No single outlier statistic identifies a body-size mechanism.

Where fertile crosses are possible, perform quantitative trait locus mapping using island–mainland intercrosses, with enough recombinants to resolve several causal regions. Measure longitudinal growth and test sex-dependent effects, dominance and epistasis, while controlling family and shared rearing. Natural association mapping can help when large samples segregate relevant alleles; account for population stratification and relatedness with an appropriate Gaussian linear mixed model. Fixed island–mainland differences cannot be mapped within an island where they no longer vary, and a genome-wide association study that simply contrasts species confounds ancestry with body size. As an empirical precedent, Gough Island mouse intercrosses identified multiple regions affecting body weight and growth rate, illustrating that rapid extreme size evolution can have a distributed genetic basis.^[ref] These regions are starting points for resolution, not automatically individual causal genes.

Integrate mapping with measurements in relevant tissues and developmental windows. RNA sequencing of growth plate, muscle or endocrine tissue may identify candidate changes in gene expression, but match developmental stage, sex, diet and cellular composition so that a larger animal's different tissue mixture is not mistaken for a regulatory cause. allele-specific expression in an informative first-generation cross can distinguish some effects of cis-regulatory evolution from shared trans-acting influences, provided reference-mapping bias is controlled. Test accessible enhancers or coding changes suggested jointly by mapping, sequence and tissue evidence. Growth pathways are reasonable hypotheses, but discovering a familiar growth-related annotation is not causal validation and should not exclude unexpected genetic loci.

Resolve the level of convergence explicitly. Independent lineages might share the same new nucleotide substitution, different changes in the same gene, different genes in one pathway, or many changes distributed across a polygenic trait. Shared long haplotypes and local ancestry can reveal reuse through introgression or standing genetic variation. Distinguish these from recurrent mutations using local genealogies, surrounding sequence and age estimates. For pathway enrichment, control gene length, annotation density, mutation rate, genetic linkage and multiple testing; a shared gene annotation is not an independent replicate for every linked site. Failure to find repeated individual hits does not exclude convergence at a broader biological level.

Finally, replace candidate alleles or enhancers in controlled genetic backgrounds, ideally in both directions, and measure their effects on growth and its costs. Use unedited siblings, sham-edited lines, independent edited founders and a rescue or precise reversion to address manipulation and off-target effects. Test combinations if mapping indicates epistasis. Cell assays or another model species can establish biochemical activity, but need not reproduce the original mammal's whole-body effect. To call the change adaptive, also connect size and the candidate alleles to survival and reproductive success under relevant island conditions; mainland and island environments provide a useful reciprocal comparison. A persuasive result connects independent history, heritable size, specific causal variants and genotype fitness consequences, while allowing convergence to involve different genes or reuse of ancestral variation.

Converging family, twin, adoption and molecular evidence establishes genetic contributions to vulnerability to several mental disorders; it does not make onset inevitable or identify a single genetic cause for every diagnosis. The outcome must be specified. Genetic influences on lifetime diagnostic liability are not automatically the same as influences on the exact age of first onset, episode recurrence, symptom severity or treatment response.

Family aggregation provides a starting point. Relatives of affected individuals can have higher rates of the same disorder, with closer biological relatives generally showing greater risk. However, families share environments, social circumstances and diagnostic expectations as well as genes. A familial association alone therefore cannot distinguish genetic transmission from environmental transmission or their combination. Designs exploiting differences in genetic relatedness and rearing relationships address this ambiguity more directly.

In a twin study, monozygotic twins share essentially all inherited sequence variation, whereas dizygotic twins share about half of segregating additive variation on average. Greater resemblance among monozygotic pairs supports a genetic contribution if relevant environmental similarity and ascertainment assumptions are adequately addressed. Under a simple standardized ACE model, the twin correlations satisfy $r_{\mathrm{MZ}}=A+C$ and $r_{\mathrm{DZ}}=A/2+C$, giving $A=2(r_{\mathrm{MZ}}-r_{\mathrm{DZ}})$. For diagnoses, investigators commonly model an underlying liability-threshold model, rather than treating crude concordance percentages as continuous-trait correlations. Patrick F. Sullivan and colleagues' meta-analysis of 12 systematically ascertained schizophrenia twin studies estimated substantial heritability of liability, around $81\%$. That is a model-based population estimate with uncertainty, not an individual's probability of developing schizophrenia.^[ref]

Twin comparisons have limitations. Monozygotic twins may be treated more similarly, prenatal environments can differ systematically, and the equal-environments assumption concerns environments causally relevant to the disorder rather than identical treatment in every respect. Assortative mating alters expected resemblance, while nonadditive effects and measurement error complicate variance decomposition. Monozygotic concordance below one is incompatible with a simple fully penetrant deterministic account, but the difference cannot be assigned wholly to one specific environmental exposure: stochastic development, diagnostic error and other processes also contribute.

Adoption and extended-family data provide a partly independent test. Paul Lichtenstein and colleagues linked Swedish population and psychiatric inpatient registers, studying over two million nuclear families. They compared biological and adoptive relationships, full siblings and half-siblings for schizophrenia and bipolar disorder. Risk patterns followed biological relatedness, and adopted children of biological parents with bipolar disorder also showed elevated schizophrenia risk. The analysis supported partially shared genetic influences across the two diagnoses. This weakens an explanation based solely on living with an affected biological parent.^[ref]

It does not eliminate every environmental explanation. Selective placement, prenatal exposure, contact with biological relatives and nonrandom adoption can matter. Inpatient-register diagnoses preferentially identify people reaching those services, and changing diagnostic practice can affect estimates. Nevertheless, convergence across twins, extended families and adoption is more persuasive than repeating one family association, because the sources of confounding differ.

Molecular evidence measures inherited variation directly. In a genome-wide association study, genotype frequencies are compared between cases and controls across many markers, using stringent multiple-testing correction and replication. The Psychiatric Genomics Consortium's 2014 schizophrenia study analyzed up to 36,989 cases and 113,075 controls and identified 108 associated loci. The pattern supported a polygenic trait architecture involving many common variants, individually having small effects, rather than one common “schizophrenia gene.” Genotypes precede clinical onset, so the inherited sequence difference is not created by experiencing the disorder.^[ref]

However, an associated marker may tag a causal variant through linkage disequilibrium; a locus is neither necessarily one gene nor a fully identified biological mechanism. Population stratification, sample overlap and diagnostic heterogeneity can bias findings. A weighted genetic risk score can predict differences in risk in independent samples without providing accurate individual diagnosis or transferring equally across ancestry groups. Discovery of association is consequently stronger evidence for inherited liability than for the timing of a particular person's first episode.

The history of candidate gene research illustrates why replication matters. Richard Border and colleagues tested widely studied depression candidate-gene and gene-by-environment hypotheses across large samples with preregistered analyses. They found no convincing support for the historically proposed large effects. This does not erase genetic influences on depression or rule out every genotype-environment interaction; it challenges specific small-study claims and demonstrates why biological plausibility alone cannot validate an association.^[ref]

Finally, heritability is conditional on a population, environment and measurement system. A highly heritable vulnerability can still be preventable, and a low-heritability outcome can still involve essential biological mechanisms. Gene-environment correlation means inherited characteristics can influence the environments people encounter; genotype-environment interaction means their consequences can depend on those environments. Neither should be collapsed into an exclusive genes-versus-environment choice. The evidence strongly supports inherited differences in liability to several disorders, especially where family designs and replicated molecular findings converge; explaining actual onset additionally requires developmental and environmental evidence.

Let $I,J$ denote the two ordered alleles at the genetic locus. Under Hardy-Weinberg equilibrium they are independent with probabilities ${\pi }_i,{\pi }_j$. Write the multiplicative penetrance as $\mathbb{P}(D\mid I=i,J=j)=k{\psi }_i{\psi }_j$, with parameters chosen so these are valid probabilities, and set $Z={\sum }_u{\pi }_u{\psi }_u$. ThenThe Bayes' theorem givesThus multiplicative penetrance preserves Hardy-Weinberg equilibrium: the affected subjects still have two independent allele draws, with tilted allele frequencies. The displayed product uses ordered allele slots. For the usual unordered genotypes, the corresponding probabilities areThe factor of two is required for heterozygotes and must not be dropped when translating the ordered notation into genotype counts.

For genetic association analysis, the case genotype probabilities factor into allele frequencies, so under this model an allelic comparison with representative population controls is appropriate; separate dominance departures are not required by the penetrance model. The ratio of case allele frequency to population allele frequency is proportional to ${\psi }_i$, and relative risk multipliers satisfy ${\psi }_i/{\psi }_j=({\pi }_i^{\ast }/{\pi }_i)/({\pi }_j^{\ast }/{\pi }_j)$. Population controls are important to this exact statement. Among specifically unaffected controls,which generally does not factor. For a rare disease it is close to the population genotype law, but no rare-disease assumption was needed for the case factorization itself. Population stratification and sampling dependence must also be addressed before using an ordinary allelic chi-squared test.

The usual transmission disequilibrium test conditions on the observed parental genotypes, rather than modelling their population probabilities. For a heterozygous parent carrying $i,j$, the two transmission outcomes are equally likely under no association. With multiplicative penetrance, ascertainment through an affected child weights them by ${\psi }_i,{\psi }_j$, soThe other parent's penetrance factor cancels. This is a special case of disease-ascertained transmission probability. Under the null ${\psi }_i={\psi }_j$, the conditional probability is $1/2$. Count the discordant transmitted/untransmitted pairs from heterozygous parents and test their balance, using McNemar's test or its exact conditional binomial form. Homozygous parental transmissions do not distinguish the alternatives.

This conditional analysis removes the nuisance population allele frequencies and retains the family matching. It is therefore robust to population stratification, which can otherwise make allele frequencies differ between cases and controls without a within-family transmission effect. It does not require the unconditional Hardy-Weinberg equilibrium and random-mating assumptions used for the population pseudo-control factorization. Independence of suitably sampled families, valid genotypes and the Mendelian null transmission model still matter.

Linkage disequilibrium means that alleles at different genetic loci are associated in population haplotypes. For two biallelic loci,The first quantity measures departure from independent allele frequencies; the second is the squared allelic correlation. Standardized linkage disequilibrium $D^{\prime }$ instead compares $D$ to its possible range at the given frequencies. High $|D^{\prime }|$ can coexist with low $r^2$, so it need not make one genetic marker a good predictor of another. Genetic linkage is a different notion: it describes within-family transmission and the recombination fraction.

The chromosome-21 study of Patil and colleagues, Science 2001 found consecutive single-nucleotide polymorphisms organized into regions with comparatively few common haplotypes; typically three haplotypes accounted for over 80% of the sampled chromosomes within a block. Such haplotype blocks can be represented efficiently by tag SNPs. Strong association within blocks and weaker association across boundaries give a useful description, though block boundaries depend on allele frequencies, population history, sampling and the operational definition of a block. They do not divide chromosomes into regions where genetic recombination is respectively impossible and certain.

The ancestry explanation uses an ancestral recombination graph. Going backwards, ancestral lineages coalesce, while a genetic recombination event splits one lineage's ancestral material into contributions from two parents. At a fixed genomic position, retaining only its ancestry gives a phylogenetic tree with the appropriate coalescent law. Nearby positions often retain much of the same tree; mutations on shared branches then produce correlated allele states. Breakpoints can change the local tree, weakening that shared history.

Under a simple random-mating two-locus model, linkage-disequilibrium decay under recombination gives $D_{t+1}=(1-r)D_t$. Actual populations also undergo genetic drift, population bottlenecks, mutation and migration, which can create or preserve association. Recombination rates themselves vary along chromosomes: McVean and colleagues, Science 2004 inferred strong local variation from sequence data. These observations motivate an ancestry model with spatially varying recombination, rather than treating physical distance as the only determinant of linkage disequilibrium.

For fine mapping, a typed genetic marker can show genetic association because it is correlated with a nearby causal genetic variant. Linkage-disequilibrium marker tagging makes a sparse first-stage scan possible, but several variants in one block can produce nearly indistinguishable association patterns. Historical recombination and diverse local genealogies can narrow the interval when denser genotyping is available.

An ancestral recombination graph turns that intuition into a likelihood calculation. One can place a proposed disease-origin event on a local ancestral branch and integrate over possible graphs, mutation histories and unobserved haplotypes. The local phylogenetic tree then relates disease-bearing chromosomes and predicts which markers should remain associated. In a Bayesian inference, uncertainty in ancestry, mutation location and demographic parameters enters the posterior distribution rather than being replaced by one guessed tree. Whole-graph inference is computationally demanding, and approximate or pairwise methods often supply practical estimates.

The useful distinction is between finding an associated region and identifying its causal variant. Strong linkage disequilibrium improves marker coverage but can limit resolution inside the associated region. Population stratification, ascertainment, rare variants and misspecified population history can also mislead an association analysis. Combining fine mapping with independent biological evidence is therefore stronger than interpreting an association peak as proof of causation.

Suppose a disease-risk mutation first occurs on an ancestral haplotype. Descendants carrying the risk allele initially share that surrounding haplotype. Historical genetic recombination progressively replaces the flanks while retaining the disease allele, leaving linkage disequilibrium between the causal variant and some nearby genetic markers. A noncausal marker can therefore have different allele frequencies in affected individuals and suitably chosen controls. This is the basis of association mapping.

One first tests dense markers or multi-marker haplotypes for genetic association with the disease phenotype. Strongly correlated single-nucleotide polymorphisms can be represented by tag SNPs to locate an associated region efficiently. Within that region, denser genotyping and examination of disease-associated haplotypes locate historical recombination breakpoints: segments retained across independently recombined disease-bearing chromosomes are stronger candidates for containing the causal allele. Compare affected and unaffected carriers as well, since incomplete penetrance and other risk factors can prevent a perfect haplotype-disease correspondence.

Association mapping can achieve fine resolution because it uses recombination accumulated over many generations, rather than only the few observed meioses in a present-day linkage pedigree. Resolution depends on the local recombination history: a long region of strong LD may identify a disease-associated block while leaving many variants statistically indistinguishable. Additional informative haplotypes or populations with different LD can help separate them. Conversely, very weak LD requires closer marker spacing to avoid missing the causal region.

Population stratification must be controlled by ancestry-matched comparisons, statistical adjustment or an appropriate family-based association design, since ancestry can affect both disease frequency and marker frequency. A marker association alone establishes neither causality nor identity of the disease gene; nearby correlated variants must be distinguished by further genetic and functional evidence. Use LD to detect a disease-associated ancestral segment, then use historical recombinants to narrow the causal interval. The correlation between causal variants and their surrounding haplotypes, together with the value of representative markers, is illustrated in the HapMap study.

Across a human chromosome, neighboring single-nucleotide polymorphisms often show substantial linkage disequilibrium, while average allelic correlation tends to decline with increasing genetic distance. A new mutation starts on one ancestral haplotype; repeated genetic recombination separates it from progressively more distant markers. In an ideal randomly mating two-locus population without further mutation, migration or selection, gamete formation gives $p_{11}^{\prime }=(1-c)p_{11}+c{p}_A{p}_B$, where $c$ is the recombination fraction. The marginal allele frequencies remain fixed, soHere the prime means the next-generation value, not a separately standardized LD statistic. This simple decay explains the distance trend, rather than asserting an exact deterministic law for real finite populations.

The local pattern is heterogeneous. Many regions have haplotype blocks with strong LD and a small number of common haplotypes, separated by intervals showing more historical recombination, often recombination hotspots. A hotspot can therefore cause a rapid drop in LD between markers that are close in physical distance. Block lengths and boundaries are not uniform, perfectly sharp or independent of the population and the statistical block definition. These empirical features are documented in the HapMap data.

Population history also matters. Population bottlenecks can reduce the number of ancestral haplotypes and increase the extent of LD; population stratification or admixture can create longer-range associations, including between distant loci. Different populations can consequently have different haplotype frequencies and LD patterns at the same genomic interval. Raw $D$ also varies with allele frequency, so comparing LD across markers requires attention to the chosen measure. The typical pattern is locally strong, block-like and highly variable LD, with a broad decline over genetic distance rather than a uniform chromosome-wide decay.

Genetic linkage analysis follows co-segregation of genetic marker and disease alleles within pedigrees. It estimates or tests a recombination fraction, with independent assortment at $\theta =1/2$ as the ordinary null. It can locate a disease genetic locus even when different families have different genetic marker alleles on their disease chromosomes.

Genetic association analysis tests whether allele or genotype frequencies differ with disease status. It can compare unrelated cases and controls, or compare transmitted and untransmitted parental alleles. A noncausal genetic marker can be associated through population linkage disequilibrium with a causal genetic locus. Association therefore concerns population-frequency relationships, whereas linkage concerns inheritance; population stratification can confound the former without producing genuine within-family linkage.

Family comparisons condition on the parental genotypes, so cases are compared with the alleles they could have inherited. This protects a transmission test from ordinary population stratification, avoids having to recruit unrelated controls matched for ancestry, and can resolve haplotypes and identify Mendelian inconsistencies. The same pedigree may also provide linkage information.

The costs are recruiting and genotyping relatives, missing or unavailable parents, and fewer informative comparisons because homozygous parents contribute no transmission contrast. These problems are especially important for late-onset disease. Related cases or multiple siblings also require an analysis accounting for dependence, and genotyping error or transmission distortion can mimic overtransmission. Population-based case/control studies are usually easier to assemble and can offer greater power per genotyped person, but need careful control of ancestry and other sources of confounding. Family sampling does not remove every form of ascertainment in a genetic study or observation bias.

For the electronic-record observational study, three strengths are:

Three weaknesses are:

For two-sample Mendelian randomization, three strengths are:

Three weaknesses are:

Thus the agreement supports the hypothesis, but does not establish the effect of a dietary intervention.

Two useful ways to strengthen the causal inference are:

These are complementary: the first changes the proposed exposure experimentally, and the second probes the validity and reproducibility of its genetic evidence.

A risk factor is a cause when an appropriate change in that factor would change the relevant psychological outcome, under specified conditions. Prediction and causation are different: a variable can forecast disorder because it is a marker, a consequence of emerging illness, or correlated with an unmeasured cause. Confidence comes from designs that address different alternative explanations and from convergence on a clearly defined exposure, population and outcome.

A cross-sectional study can establish an association and suggest hypotheses, but timing is ambiguous. A depressed respondent may recall more adversity, and family disadvantage may influence both exposure and symptoms. A case–control study is efficient for uncommon disorders but remains sensitive to recall, control selection and exposure misclassification. Adjusting for many measured variables is not automatically an improvement: controlling a consequence of exposure can block a real pathway, while conditioning on a common consequence of exposure and outcome can induce collider bias.

Prospective longitudinal studies improve temporal ordering. In Constance Hammen and colleagues' two-year follow-up, childhood-adversity histories predicted depressive reactions at lower subsequent psychological stress levels. Measuring events and symptoms over time made this more informative than asking affected people retrospectively why they became ill. Yet adversity was not randomized, and histories may be measured imperfectly. Temporal precedence is necessary for the proposed direction of influence, but does not remove familial liability, prior symptoms or unmeasured social circumstances as competing explanations.^[ref]

Stronger observational designs make their counterfactual comparison explicit. A natural experiment exploits an external change that alters exposure without being produced by the same individual characteristics that usually determine it. E. Jane Costello and colleagues observed children before and after an income supplement associated with a casino opening. Comparing developmental changes across differently affected families strengthened the social-causation interpretation for particular behavioral outcomes. It did not randomize poverty, and concurrent community changes or selection into an income-defined subgroup remain possible alternatives.^[ref]

A difference-in-differences analysis formalizes such reasoning as the exposed group's change minus the comparison group's change. Its causal interpretation requires a credible untreated trend comparison, no differential concurrent change that explains the effect, and careful attention to anticipation and changes in sample composition. Merely measuring before and after a policy is not sufficient. Examining pre-intervention trends can expose some problems but cannot prove that every counterfactual assumption is true.

A randomized controlled trial offers another route by breaking systematic baseline links between allocation and measured or unmeasured causes. The foster-care offer in the Bucharest Early Intervention Project is an example of ethically intervening to improve a harmful setting, rather than deliberately assigning children to initial deprivation. Under suitable follow-up and analysis, differences between allocated groups estimate the effect of offering that intervention. Nonadherence, placement changes and attrition still matter. An intention-to-treat analysis retains the allocation comparison; simply regrouping children by the care they eventually received can reintroduce selection.^[ref]

Even a successful trial does not identify the causal contribution of every component. A foster-care package changes several experiences simultaneously. A coping course may alter attention, expectations, peer contact and behavior. To identify a mediator, it is necessary to consider confounding of the mediator–outcome relationship and its temporal ordering, not only whether the treatment changes both mediator and outcome. Conversely, a null intervention can reflect inadequate delivery, insufficient contrast or the wrong developmental timing; it need not prove that the proposed risk factor is irrelevant under every possible intervention.

Genetically informed studies address another set of alternatives. Comparing differently exposed identical twins can reduce confounding by inherited sequence and shared family background. However, exposure differences can arise from nonshared antecedents that also affect disorder, and within-pair measurement error can be substantial. Family and adoption designs are similarly informative without making every remaining environmental association causal. Molecular associations can indicate inherited liability, but linkage disequilibrium and population stratification separate an associated marker from a demonstrated mechanism.

Mendelian randomization uses genetic variants as instruments for an exposure, requiring relevant instruments, independence from exposure–outcome confounders and an exclusion restriction. The design can be useful when its assumptions are credible. It is not a shortcut for every social risk factor: horizontal pleiotropy, parental genetic effects and ancestry structure can create misleading estimates. A genetic score associated with stressful experiences does not necessarily act only through those experiences. This is a methodological distinction, not evidence that any particular untested exposure is causal.

Animal experiments can manipulate mechanisms and timing more tightly. Karen J. Parker and colleagues randomized young squirrel monkeys to intermittent separation experiences or a comparison condition and later tested responses in a novel setting. This supports causal conclusions about that rearing contrast and those measured responses. It does not directly estimate the effect of childhood adversity on human psychiatric diagnoses: species, development, psychological stressor meaning and outcomes differ.^[ref]

The strongest conclusion is obtained when these designs converge despite their different weaknesses, with sensitivity analyses and transparent outcome definitions. A causal claim should specify the contrast it establishes and the alternatives it has not excluded. Prospective timing, externally driven exposure changes, randomization, within-family controls and mechanistic experiments together can provide much stronger evidence than any isolated risk correlation, while still leaving uncertainty about transport to another population or developmental stage.

Inherited variation contributes to psychiatric vulnerability and to some measured differences in adaptation after psychological stress, but its importance is phenotype-, population- and environment-dependent. Genetic influence is not equivalent to inevitability, a single large-effect gene or a wholly separate set of protective genes. A satisfactory answer must distinguish genetic effects on exposure, on symptoms and on responses conditional on exposure.

Family resemblance provides an initial clue, but families transmit environments and social circumstances as well as DNA. Twin comparisons exploit differences in relatedness: monozygotic pairs share essentially all inherited sequence variation, whereas dizygotic pairs share about half their segregating additive variation on average. Under an ACE model, greater resemblance of monozygotic pairs can support a genetic component if the relevant environmental and ascertainment assumptions hold. Equal environments means comparable effects of relevant environmental similarity, not literally identical treatment of every twin.

Disorder-level evidence is substantial. Patrick F. Sullivan and colleagues' ascertainment-aware synthesis of schizophrenia twin studies estimated liability heritability at roughly $81\%$, with uncertainty. Paul Lichtenstein and colleagues' Swedish family-register study supported partially shared inherited liability to schizophrenia and bipolar disorder using biological, adoptive and sibling relationships. The convergence supports inherited vulnerability while showing that diagnostic categories do not necessarily partition genetic causes into entirely independent sets.^[ref][ref]

Molecular evidence provides a different test. The Psychiatric Genomics Consortium's 2014 schizophrenia study identified associations spanning 108 loci in a large multi-stage genome-wide analysis. Many common variants contribute small effects; an associated locus is not automatically one causal gene or a known mechanism. Population stratification, linkage disequilibrium, case ascertainment and sample overlap still require control. A genetic risk score can predict risk differences in independent data while remaining far from a deterministic diagnosis or a transportable measure of one person's psychological resilience.^[ref]

Psychological resilience needs additional definition because favorable outcomes without adversity do not establish psychological stress resistance. Ananda B. Amstadter and colleagues studied 7,500 adult twins at two assessments, modeling internalizing symptoms relative to cumulative stressful events. Their regression residual-based phenotype had estimated heritability around $31\%$ at each wave; longitudinal modeling that accounted for measurement error and occasion-specific influences gave a larger latent estimate, around $50\%$. This supports both inherited and environmental contributions to that operationalization, not a universal percentage for every meaning of psychological resilience.^[ref]

The shift in estimates is conceptually informative. Independent measurement error adds to observed variance without adding the corresponding genetic signal. In a simple measurement model, if $V_G$ is genetic variance, $V_T$ true-score variance and $V_{ϵ}$ independent error variance, thenRemoving noise can therefore raise estimated heritability without changing anyone's genes. This illustration is not a reconstruction of the study's full model: occasion-specific effects and longitudinal assumptions also matter, so $31/50$ should not be reported as an independently established measurement reliability coefficient.

A regression residual phenotype also inherits assumptions from its construction. If adversity is undermeasured or a prediction model is misspecified, the regression residual includes more than protective capacity. Genetic effects on depression or anxiety can appear as genetic effects on low symptom burden after psychological stress, even without a separate buffering mechanism. Demonstrating a psychological resilience-specific effect would require stronger evidence that inherited variation changes the response to comparable exposure, with adequate exposure measurement and explicit competing models.

Gene–environment processes further complicate simple partitioning. Gene-environment correlation occurs when inherited characteristics influence environments encountered or elicited, or when parental characteristics affect both inherited liability and the rearing environment. Genotype-environment interaction instead means that inherited differences have different consequences across environments. Neither can be inferred merely because both genes and psychological stress predict the outcome. Restricting analysis to people with a particular exposure can also create selection bias if exposure is influenced by several causes relevant to symptoms.

Historical candidate-gene findings illustrate the need for stringent evidence. Richard Border and colleagues' large, preregistered analyses did not support widely proposed depression candidate-gene or candidate-gene-by-environment effects. This challenges those specific claims and the power of small discovery studies; it does not demonstrate that all interactions are absent. Replication, correction for multiple testing, ancestry-aware analyses and precise phenotype definitions are essential before a plausible biological story becomes a reliable finding.^[ref]

Finally, high heritability does not mean an outcome cannot be changed. It describes sources of variation under a particular range of environments, not the response to an environment absent from the data. Genetic effects can act through modifiable developmental systems, and external support can benefit people across inherited-risk levels. Predictive accuracy and equity can also deteriorate when scores are transferred to populations unlike the discovery sample.

Genes contribute meaningfully to liability and to particular measured adaptation phenotypes, while evidence for specific protective mechanisms is narrower than evidence for heritability. The strongest account combines family and molecular results with explicit exposure and developmental models, rather than assigning a fixed genetic percentage to an individual or replacing environmental prevention with genetic classification.