Uppsala University and GMI

What are pleiotropic mechanisms?

Broad sense pleiotropy

  • Alleles affect multiple traits

Genotype-by-environment interactions

  • Falconer (1952): G-by-E is a special case of pleiotropy

"Three types" of pleiotropy

Negative pleiotropy

Allele increases one trait, but decreases the other

"Three types" of pleiotropy

Positive pleiotropy

Allele 2 increases both traits

"Three types" of pleiotropy

Non-pleiotropy

  • Effect in context A, no effect in context B
  • Also called conditional neutrality

"Three types" of pleiotropy

Which of these mechanisms contibute to adaptation?

Mechanisms are a continuum

  • Pleiotropy "types" are special cases
  • Quantitative questions:
    • How costly is the wrong allele?
    • How difficult would it be to compensate?

How do we usually study pleiotropy?

Null-hypothesis tests cause bias

What we usually do

Example: mapping fitness QTL in A. thaliana (Ellis et al, in prep.)

Effect on fitness of Swedish allele in Italy (red) and Sweden (blue)

Plot allelic effects on two axes

(fictional data)

Get p-values in each context

Which loci have p<0.05 in one or both contexts?

Bias towards strongest effects

Hill & Zhang (2012), Genetics 190:1131–1137

Ignores quantitative information

Underestimates pleiotropy

Harder to get two significant results than one

Aims

Null-hypothesis tests:

  • exclude weak but important effects
  • throw away quantitative information
  • underestimate pleiotropy

Aims

Null-hypothesis tests:

  • excludes weak, but important effects
  • throw away quantitative information
  • underestimate pleiotropy

It would be better to:

  • Use all the data
  • Quantify rather than classify
  • Estimate effect sizes and their uncertainty

Quantifying pleiotropy

Pleiotropic effects as vectors

Vector length \(z\) \(\rightarrow\) magnitude of effect

Pleiotropic effects as vectors

Angle reflects mechanism

Sine wave of the angle

\(\sin 2 \theta\)

Positive pleiotropy

\(\sin 2 \theta\)

Negative pleiotropy

\(\sin 2 \theta\)

Non-pleiotropy

\(\sin 2 \theta\)

Linearise the function

\[ \tau = \frac{\arcsin (\sin(2 \theta))}{\arcsin(1)} \]

Uncertainty in \(\tau\) and \(z\)

  • Non-parametric bootstrap of effects in each context
  • Distribution of \(\tau\) and vector lengths

Uncertainty in \(\tau\) and \(z\)

  • Average over bootstraps in downstream analyses.

Greater confidence in longer vectors

  • Does not depend on the direction of the effect!

\(\tau\) quantifies pleiotropic mechanism

  • 1 = positive pleiotropy
  • 0 = non-pleiotropy (conditional neutrality)
  • -1 = negative pleiotropy

\(\tau\) quantifies pleiotropic mechanism

  • 1 = positive pleiotropy
  • 0 = non-pleiotropy (conditional neutrality)
  • -1 = negative pleiotropy
  • (matches direction of genetic correlation)

\(\tau\) quantifies pleiotropic mechanism

  • 1 = positive pleiotropy
  • 0 = non-pleiotropy (conditional neutrality)
  • -1 = negative pleiotropy
  • (matches direction of genetic correlation)
  • And everything in between!

Application to data

GxE for fitness and local adaptation

Local adaptation in plants

  • Data from 4 reciprocal transplant experiments:
    • Arabidopsis thaliana (Ågren et al. 2013)
    • Avena barbata (Latta et al. 2009)
    • Boechera stricta (Anderson et al. 2014)
    • Hordeum spontaneum (Verhoeven et al. 2004)
  • F2-derived linkage-mapping populations
  • Genotyped at 100-350 markers
  • Log relative fitness of each marker

Widespread trade-offs in A. thaliana

  • Histrogram of \(\tau\) across all markers

Most of the genome shows pleiotropy

  • Trade-offs across most of chromosomes 3, 4, 5
  • More extensive than were detected through mapping.

Substantial global superiority

  • Histograms shifted right

Substantial global superiority

  • Dispersal limits spread?
  • Genetic load?

Global effects in Boechera

  • Anderson et al. (2013), "Genetic trade‐offs and conditional neutrality contribute to local adaptation", Mol. Ecol.

Summary and outlook

Summary

  • Pleiotropy is a continuum
  • Statistical tests underestimate pleiotropy

Summary

  • Pleiotropy is a continuum
  • Statistical tests underestimate pleiotropy

  • Using \(\tau\) allows you to quantify pleiotropic mechanisms
  • Estimate effect and quantify uncertainty
  • Does not depend on the direction of the effect

Applications

  • Applicable whenever you can quantify 'effects':
    • Phenotypes of individuals/strains/lines
    • Allelic effects
    • Wherever you can measure a 'reaction norm'

Applications

  • Applicable whenever you can quantify pairs of 'effects':
    • Phenotypes of individuals/strains/sexes
    • Allelic effects
    • Wherever you can measure a 'reaction norm'
  • True for pleiotropy and GxE

Applications

  • Applicable whenever you can estimate a reaction norm:
    • Phenotypes of individuals/strains/lines
    • Allelic effects
    • Wherever you can measure a 'reaction norm'
  • True for pleiotropy and GxE
  • R package sintillate available:

Extensions

  • Correlations in large phenotype databases
  • Extensible to multivariate pleiotropy, with caution
  • Quantify epistasis and dominance

Acknowledgements

  • A. thaliana project:
    • Jon Ågren (Uppsala)
    • Doug Schemske (MSU)
    • Chris Oakley (MSU, Purdue University)
    • Many helpers, especially Mattias Vass, Linus Vikström, Jenny Glans
  • And Magnus Nordborg and the GMI

Survival and fecundity

  • Fitness components in A. thaliana (Ågren et al. 2013)
    • Reciprocal transplants in Italy and Sweden
    • Recombinant inbred lines + local parents
    • Line mean survival and fecundity
    • 3 years

RIL means are positively correlated

  • RIL = recombinant inbred line, selfed from an F2 cross

Correlation reflects survival

  • RIL = recombinant inbred line, selfed from an F2 cross

Parental phenotypes follow the RILs

  • Link parents, RIL phenotypes, and marker effects.

Genetic basis of trait correlations

  • Correlations due to pleiotropy
    • Same genes affect both traits
    • \(\tau > 0\)
  • Or:
    • Different genes affect each trait
    • \(\tau \sim 0\)

Effect sizes

  • Log relative fitness of the Swedish allele at each marker:
    • Univariate
    • Includes effects of linked loci

Distributions of \(\tau\)

  • Observed (bars) with bootstrap means (points)

Non-pleiotropy in 2009

  • \(\tau \sim 0 \rightarrow\) non-pleiotropy

Positive pleiotropy in 2010

  • Distribution shifted right \(\rightarrow\) positive pleiotropy

A bit less pleiotropy in 2011

  • \(\tau\) between 0 and 1.

Loci predict parental phenotypes

  • Observed \(\tau\) for parental phenotypes ± 95% confidence intervals

Pleiotropy across the genome

  • Grey lines = 95% confidence intervals from 1000 bootstraps ## Rescaling vector norm

  • Scale: relative fitness

Italy vs Sweden

Italy vs Sweden

Relative fitness