Genome-wide inference of ancestral recombination graphs.
The complex correlation structure of a collection of orthologous DNA sequences is uniquely captured by the "ancestral recombination graph" (ARG), a complete record of coalescence and recombination events in the history of the sample. However, existing methods for ARG inference are computat...
| Main Authors: | , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Public Library of Science (PLoS)
2014-01-01
|
| Series: | PLoS Genetics |
| Online Access: | http://europepmc.org/articles/PMC4022496?pdf=render |
| _version_ | 1818138610734465024 |
|---|---|
| author | Matthew D Rasmussen Melissa J Hubisz Ilan Gronau Adam Siepel |
| author_facet | Matthew D Rasmussen Melissa J Hubisz Ilan Gronau Adam Siepel |
| author_sort | Matthew D Rasmussen |
| collection | DOAJ |
| description | The complex correlation structure of a collection of orthologous DNA sequences is uniquely captured by the "ancestral recombination graph" (ARG), a complete record of coalescence and recombination events in the history of the sample. However, existing methods for ARG inference are computationally intensive, highly approximate, or limited to small numbers of sequences, and, as a consequence, explicit ARG inference is rarely used in applied population genomics. Here, we introduce a new algorithm for ARG inference that is efficient enough to apply to dozens of complete mammalian genomes. The key idea of our approach is to sample an ARG of [Formula: see text] chromosomes conditional on an ARG of [Formula: see text] chromosomes, an operation we call "threading." Using techniques based on hidden Markov models, we can perform this threading operation exactly, up to the assumptions of the sequentially Markov coalescent and a discretization of time. An extension allows for threading of subtrees instead of individual sequences. Repeated application of these threading operations results in highly efficient Markov chain Monte Carlo samplers for ARGs. We have implemented these methods in a computer program called ARGweaver. Experiments with simulated data indicate that ARGweaver converges rapidly to the posterior distribution over ARGs and is effective in recovering various features of the ARG for dozens of sequences generated under realistic parameters for human populations. In applications of ARGweaver to 54 human genome sequences from Complete Genomics, we find clear signatures of natural selection, including regions of unusually ancient ancestry associated with balancing selection and reductions in allele age in sites under directional selection. The patterns we observe near protein-coding genes are consistent with a primary influence from background selection rather than hitchhiking, although we cannot rule out a contribution from recurrent selective sweeps. |
| first_indexed | 2024-12-11T10:14:56Z |
| format | Article |
| id | doaj.art-bbc106f8dc4d452d993072ec627fe12a |
| institution | Directory Open Access Journal |
| issn | 1553-7390 1553-7404 |
| language | English |
| last_indexed | 2024-12-11T10:14:56Z |
| publishDate | 2014-01-01 |
| publisher | Public Library of Science (PLoS) |
| record_format | Article |
| series | PLoS Genetics |
| spelling | doaj.art-bbc106f8dc4d452d993072ec627fe12a2022-12-22T01:11:39ZengPublic Library of Science (PLoS)PLoS Genetics1553-73901553-74042014-01-01105e100434210.1371/journal.pgen.1004342Genome-wide inference of ancestral recombination graphs.Matthew D RasmussenMelissa J HubiszIlan GronauAdam SiepelThe complex correlation structure of a collection of orthologous DNA sequences is uniquely captured by the "ancestral recombination graph" (ARG), a complete record of coalescence and recombination events in the history of the sample. However, existing methods for ARG inference are computationally intensive, highly approximate, or limited to small numbers of sequences, and, as a consequence, explicit ARG inference is rarely used in applied population genomics. Here, we introduce a new algorithm for ARG inference that is efficient enough to apply to dozens of complete mammalian genomes. The key idea of our approach is to sample an ARG of [Formula: see text] chromosomes conditional on an ARG of [Formula: see text] chromosomes, an operation we call "threading." Using techniques based on hidden Markov models, we can perform this threading operation exactly, up to the assumptions of the sequentially Markov coalescent and a discretization of time. An extension allows for threading of subtrees instead of individual sequences. Repeated application of these threading operations results in highly efficient Markov chain Monte Carlo samplers for ARGs. We have implemented these methods in a computer program called ARGweaver. Experiments with simulated data indicate that ARGweaver converges rapidly to the posterior distribution over ARGs and is effective in recovering various features of the ARG for dozens of sequences generated under realistic parameters for human populations. In applications of ARGweaver to 54 human genome sequences from Complete Genomics, we find clear signatures of natural selection, including regions of unusually ancient ancestry associated with balancing selection and reductions in allele age in sites under directional selection. The patterns we observe near protein-coding genes are consistent with a primary influence from background selection rather than hitchhiking, although we cannot rule out a contribution from recurrent selective sweeps.http://europepmc.org/articles/PMC4022496?pdf=render |
| spellingShingle | Matthew D Rasmussen Melissa J Hubisz Ilan Gronau Adam Siepel Genome-wide inference of ancestral recombination graphs. PLoS Genetics |
| title | Genome-wide inference of ancestral recombination graphs. |
| title_full | Genome-wide inference of ancestral recombination graphs. |
| title_fullStr | Genome-wide inference of ancestral recombination graphs. |
| title_full_unstemmed | Genome-wide inference of ancestral recombination graphs. |
| title_short | Genome-wide inference of ancestral recombination graphs. |
| title_sort | genome wide inference of ancestral recombination graphs |
| url | http://europepmc.org/articles/PMC4022496?pdf=render |
| work_keys_str_mv | AT matthewdrasmussen genomewideinferenceofancestralrecombinationgraphs AT melissajhubisz genomewideinferenceofancestralrecombinationgraphs AT ilangronau genomewideinferenceofancestralrecombinationgraphs AT adamsiepel genomewideinferenceofancestralrecombinationgraphs |