Cargando…
GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes
The [Formula: see text] ratio provides evidence of adaptation or functional constraint in protein-coding genes by quantifying the relative excess or deficit of amino acid-replacing versus silent nucleotide variation. Inexpensive sequencing promises a better understanding of parameters, such as [Form...
Autor principal: | |
---|---|
Formato: | Online Artículo Texto |
Lenguaje: | English |
Publicado: |
Oxford University Press
2020
|
Materias: | |
Acceso en línea: | https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7403622/ https://www.ncbi.nlm.nih.gov/pubmed/32167543 http://dx.doi.org/10.1093/molbev/msaa069 |
_version_ | 1783566978229731328 |
---|---|
author | Wilson, Daniel J |
author_facet | Wilson, Daniel J |
author_sort | Wilson, Daniel J |
collection | PubMed |
description | The [Formula: see text] ratio provides evidence of adaptation or functional constraint in protein-coding genes by quantifying the relative excess or deficit of amino acid-replacing versus silent nucleotide variation. Inexpensive sequencing promises a better understanding of parameters, such as [Formula: see text] , but analyzing very large data sets poses a major statistical challenge. Here, I introduce genomegaMap for estimating within-species genome-wide variation in [Formula: see text] , and I apply it to 3,979 genes across 10,209 tuberculosis genomes to characterize the selection pressures shaping this global pathogen. GenomegaMap is a phylogeny-free method that addresses two major problems with existing approaches: 1) It is fast no matter how large the sample size and 2) it is robust to recombination, which causes phylogenetic methods to report artefactual signals of adaptation. GenomegaMap uses population genetics theory to approximate the distribution of allele frequencies under general, parent-dependent mutation models. Coalescent simulations show that substitution parameters are well estimated even when genomegaMap’s simplifying assumption of independence among sites is violated. I demonstrate the ability of genomegaMap to detect genuine signatures of selection at antimicrobial resistance-conferring substitutions in Mycobacterium tuberculosis and describe a novel signature of selection in the cold-shock DEAD-box protein A gene deaD/csdA. The genomegaMap approach helps accelerate the exploitation of big data for gaining new insights into evolution within species. |
format | Online Article Text |
id | pubmed-7403622 |
institution | National Center for Biotechnology Information |
language | English |
publishDate | 2020 |
publisher | Oxford University Press |
record_format | MEDLINE/PubMed |
spelling | pubmed-74036222020-08-07 GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes Wilson, Daniel J Mol Biol Evol Resources The [Formula: see text] ratio provides evidence of adaptation or functional constraint in protein-coding genes by quantifying the relative excess or deficit of amino acid-replacing versus silent nucleotide variation. Inexpensive sequencing promises a better understanding of parameters, such as [Formula: see text] , but analyzing very large data sets poses a major statistical challenge. Here, I introduce genomegaMap for estimating within-species genome-wide variation in [Formula: see text] , and I apply it to 3,979 genes across 10,209 tuberculosis genomes to characterize the selection pressures shaping this global pathogen. GenomegaMap is a phylogeny-free method that addresses two major problems with existing approaches: 1) It is fast no matter how large the sample size and 2) it is robust to recombination, which causes phylogenetic methods to report artefactual signals of adaptation. GenomegaMap uses population genetics theory to approximate the distribution of allele frequencies under general, parent-dependent mutation models. Coalescent simulations show that substitution parameters are well estimated even when genomegaMap’s simplifying assumption of independence among sites is violated. I demonstrate the ability of genomegaMap to detect genuine signatures of selection at antimicrobial resistance-conferring substitutions in Mycobacterium tuberculosis and describe a novel signature of selection in the cold-shock DEAD-box protein A gene deaD/csdA. The genomegaMap approach helps accelerate the exploitation of big data for gaining new insights into evolution within species. Oxford University Press 2020-08 2020-03-13 /pmc/articles/PMC7403622/ /pubmed/32167543 http://dx.doi.org/10.1093/molbev/msaa069 Text en © The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. http://creativecommons.org/licenses/by/4.0/ This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited. |
spellingShingle | Resources Wilson, Daniel J GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes |
title |
GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes |
title_full |
GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes |
title_fullStr |
GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes |
title_full_unstemmed |
GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes |
title_short |
GenomegaMap: Within-Species Genome-Wide d(N)/d(S) Estimation from over 10,000 Genomes |
title_sort | genomegamap: within-species genome-wide d(n)/d(s) estimation from over 10,000 genomes |
topic | Resources |
url | https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7403622/ https://www.ncbi.nlm.nih.gov/pubmed/32167543 http://dx.doi.org/10.1093/molbev/msaa069 |
work_keys_str_mv | AT wilsondanielj genomegamapwithinspeciesgenomewidedndsestimationfromover10000genomes AT genomegamapwithinspeciesgenomewidedndsestimationfromover10000genomes |