Aoife McLysaght Dept. of Genetics Trinity College Dublin

Transcription

1 Aoife McLysaght Dept. of Genetics Trinity College Dublin

2 Evolution of genome arrangement Evolution of genome content.

3 Evolution of genome arrangement Gene order changes Inversions, translocations Evolution of genome content.

4 Evolution of genome arrangement Gene order changes Inversions, translocations Evolution of genome content Gene gain (sequence divergence, duplication, recombination, horizontal transfer) Gene loss (deletion).

5 Evolution of genome arrangement Gene order changes Inversions, translocations Evolution of genome content Gene gain (sequence divergence, duplication, recombination, horizontal transfer) Gene loss (deletion) One or more genes per event

6 Translate knowledge from sequenced or model genomes to organism of interest Positional cloning of genes Use probes designed in one genome to detect a target in another genome Improve model parameters for phylogenetic inference from genome arrangement

7 Not just a bag of genes Genome organisation contains information Order of Hox genes corresponds to spatial pattern of gene expression Clustering of housekeeping genes By observation of allowed changes gain understanding of genomic constraints and plasticity

8 Greater power to detect change Precision Can infer lineage in which change occurred Detect direction and rate of change More genomes also increase computational burden

9 20 completely sequenced genomes kb containing ~200 genes

10 Double-stranded DNA viruses, no RNA stage Replicate in the host cytoplasm Entomopox insect infecting Chordopox vertebrate infecting Orthopox subset of chordopox which includes smallpox (variola) and vaccinia

11 How are these genomes arranged? How has genome content changed? Is the rate of change constant?

12 How are these genomes arranged? How has genome content changed? Is the rate of change constant? Can we detect adaptive genome evolution?

13 Significant sequence similarity How significant? over a long stretch of the protein How long?

14 e-20 1e-10 1e-5 1 e-value threshold Minimum aligned proportion

15 Complete linkage Single-link clustering Our method

16 C A E D B

17 C A E G F D B

18 C A E F J D C B E G I D H B

19 4042 total proteins 3384 proteins classified into 875 groups 813 complete linkage 521 groups of 1 member 150 groups of 2 members members

20

21 34 orthologues present in all genomes

22 34 orthologues present in all genomes

23 92 orthologues present in all orthopox genomes

24 Examine phylogenetic spread of a group of orthologues Assign gene gain and loss events to branches in the phylogeny

25

26

27

28

29

30 Tested for uniform rate of gene acquisition Assume a molecular clock

31 Tested for uniform rate of gene acquisition Assume a molecular clock Are gene acquisition events distributed randomly throughout the tree?

32 Tested for uniform rate of gene acquisition Assume a molecular clock Are gene acquisition events distributed randomly throughout the tree? Simulations

33 Significant deficit Significant excess

34 Slower rate of amino acid substitution within this clade (leading to abberantly short branch lengths) Takezaki relative rate test Branch lengths from synonymous distances

35 Slower rate of amino acid substitution within this clade (leading to abberantly short branch lengths) Takezaki relative rate test Branch lengths from synonymous distances Increased rate of gene gain Increased selection for the retention of gained genes

36 Extensive sequence divergence Recombination Horizontal transfer

37 AMV-EPB_034 inhibitor of apoptosis from Amsacta moorei entomopoxvirus (AMV-EPB) GenBank sequence inhibitor of apoptosis from Bombyx mori (silkworm) BLAST e-value 9e-81 Amsacta moorei entomopoxvirus infects Amsacta moorei (Red Hairy Caterpillar) Bombyx and Amsacta both Order Lepidoptera

38 AMV-EPB_034 inhibitor of apoptosis from Amsacta moorei entomopoxvirus (AMV-EPB) GenBank sequence inhibitor of apoptosis from Bombyx mori (silkworm) BLAST e-value 9e-81 Amsacta moorei entomopoxvirus infects Amsacta moorei (Red Hairy Caterpillar) Bombyx and Amsacta both Order Lepidoptera 62% of best non-viral GenBank hits are from same taxonomic Class as viral host

39 Events are not independent Depend on previous (in time) gain and loss events of the gene family Requires a probabilistic model?

40

41 Selection for diversification Positive selection Characteristic of host-parasite co-evolution

42 GGG GAG GCG GUG GGA GAA Glu GCA GUA GGC GAC GCC GUC GGU Gly GAU Asp GCU Ala GUU Val AGG AAG ACG AUG Met AGA Arg AAA Lys ACA AUA AGC AAC ACC AUC AGU Ser AAU Asn ACU Thr AUU Ile CGG CAG CCG CUG CGA CAA Gln CCA CUA CGC CAC CCC CUC CGU Arg CAU His CCU Pro CUU Leu UGG Trp UAG ter UCG UUG UGA ter UAA ter UCA UUA Leu UGC UAC UCC UUC UGU Cys UAU Tyr UCU Ser UUU Phe

43 Two classes of DNA substitutions Synonymous (DNA change without amino acid change) Nonsynonymous (DNA change causing amino acid change) Neutral equal frequencies Conservative selection fewer nonsynonymous substitutions Positive selection more nonsynonymous substitutions

44 Two classes of DNA substitutions Synonymous (DNA change without amino acid change) Nonsynonymous (DNA change causing amino acid change) Neutral equal frequencies Conservative selection fewer nonsynonymous substitutions Positive selection more nonsynonymous substitutions

45 Two classes of DNA substitutions Synonymous (DNA change without amino acid change) Nonsynonymous (DNA change causing amino acid change) Neutral equal frequencies Conservative selection fewer nonsynonymous substitutions Positive selection more nonsynonymous substitutions

46 Two classes of DNA substitutions Synonymous (DNA change without amino acid change) Nonsynonymous (DNA change causing amino acid change) Neutral equal frequencies Conservative selection fewer nonsynonymous substitutions Positive selection more nonsynonymous substitutions

47 204 groups of orthologues Maximum liklihood test for positive selection (PAML) Significantly higher frequency of nonsynonymous substitutions

48 Detected positive selection on 26 genes Examples: Membrane glycoprotein Haemagluttinin Immunoglobulin domain protein

49 13 genes are unique to orthopox clade Significantly more than expected (P < 0.05) Disproportionate frequency of positive selection on genes gained within the orthopox lineage

50 Association of positive selection on protein sequences and increased rate of gene acquisition

51 Association of positive selection on protein sequences and increased rate of gene acquisition Adaptive significance of gene acquisition? Mimic host defences Avoid host recognition Block cell death

52 The rate of genome evolution is not constant The rate of gene acquisition has increased in the orthopox lineage Orthopox lineage is also has an increased frequency of positive selection Possible adaptive significance of genome evolution

53 University of California, Irvine Brandon Gaut Pierre Baldi