How Often Do Protein Genes Navigate Valleys of Low Fitness?

To escape from local fitness peaks, a population must navigate across valleys of low fitness. How these transitions occur, and what role they play in adaptation, have been subjects of active interest in evolutionary genetics for almost a century. However, to our knowledge, this problem has never bee...

Full description

Bibliographic Details
Main Authors: Erik D. Nelson, Nick V. Grishin
Format: Article
Language:English
Published: MDPI AG 2019-04-01
Series:Genes
Subjects:
Online Access:https://www.mdpi.com/2073-4425/10/4/283
Description
Summary:To escape from local fitness peaks, a population must navigate across valleys of low fitness. How these transitions occur, and what role they play in adaptation, have been subjects of active interest in evolutionary genetics for almost a century. However, to our knowledge, this problem has never been addressed directly by considering the evolution of a gene, or group of genes, as a whole, including the complex effects of fitness interactions among multiple loci. Here, we use a precise model of protein fitness to compute the probability <inline-formula> <math display="inline"> <semantics> <mrow> <mi>P</mi> <mo>(</mo> <mi>s</mi> <mo>,</mo> <mo>&#916;</mo> <mi>t</mi> <mo>)</mo> </mrow> </semantics> </math> </inline-formula> that an allele, randomly sampled from a population at time <i>t</i>, has crossed a fitness valley of depth <i>s</i> during an interval <inline-formula> <math display="inline"> <semantics> <mfenced separators="" open="[" close="]"> <mi>t</mi> <mo>&#8722;</mo> <mo>&#916;</mo> <mi>t</mi> <mo>,</mo> <mi>t</mi> </mfenced> </semantics> </math> </inline-formula> in the immediate past. We study populations of model genes evolving under equilibrium conditions consistent with those in mammalian mitochondria. From this data, we estimate that genes encoding small protein motifs navigate fitness valleys of depth <inline-formula> <math display="inline"> <semantics> <mrow> <mn>2</mn> <mi>N</mi> <mi>s</mi> <mo>≳</mo> <mn>30</mn> </mrow> </semantics> </math> </inline-formula> with probability <inline-formula> <math display="inline"> <semantics> <mrow> <mi>P</mi> <mo>≳</mo> <mn>0.1</mn> </mrow> </semantics> </math> </inline-formula> on a time scale of human evolution, where <i>N</i> is the (mitochondrial) effective population size. The results are consistent with recent findings for Watson&#8211;Crick switching in mammalian mitochondrial tRNA molecules.
ISSN:2073-4425