The Role of the Effective Population Size in Compensatory Evolution

The impact of the effective population size (Ne) on the efficacy of selection has been the focus of many theoretical and empirical studies over the recent years. Yet, the effect of Ne on evolution under epistatic fitness interactions is not well understood. In this study, we compare selective constraints at independently evolving (unpaired) and coevolving (paired) sites in orthologous transfer RNAs (tRNA molecules for vertebrate and drosophilid species pairs of different Ne. We show that patterns of nucleotide variation for the two classes of sites are explained well by Kimura's one- and two-locus models of sequence evolution under mutational pressure. We find that constraints in orthologous tRNAs increase with increasing Ne of the investigated species pair. Thereby, the effect of Ne on the efficacy of selection is stronger at unpaired sites than at paired sites. Furthermore, we identify a “core” set of tRNAs with high structural similarity to tRNAs from all major kingdoms of life and a “peripheral” set with lower similarity. We observe that tRNAs in the former set are subject to higher constraints and less prone to the effect of Ne, whereas constraints in tRNAs of the latter set show a large influence of Ne. Finally, we are able to demonstrate that constraints are relaxed in X-linked drosophilid tRNAs compared with autosomal tRNAs and suggest that Ne is responsible for this difference. The observed effects of Ne are consistent with the hypothesis that evolution of most tRNAs is governed by slightly to moderately deleterious mutations (i.e., |Nes| ≤ 5).

[1]  J. Jurka,et al.  Repbase Update, a database of eukaryotic repetitive elements , 2005, Cytogenetic and Genome Research.

[2]  P. Stadler,et al.  Secondary structure prediction for aligned RNA sequences. , 2002, Journal of molecular biology.

[3]  T. Jukes,et al.  The neutral theory of molecular evolution. , 2000, Genetics.

[4]  D. Petrov,et al.  Similar Levels of X-linked and Autosomal Nucleotide Variation in African and non-African populations of Drosophila melanogaster , 2007, BMC Evolutionary Biology.

[5]  Adam Eyre-Walker,et al.  Patterns of evolutionary constraints in intronic and intergenic DNA of Drosophila. , 2004, Genome research.

[6]  B. Charlesworth.  Effective population size and patterns of molecular evolution and variation , 2009, Nature Reviews Genetics.

[7]  J H Gillespie,et al.  The role of population size in molecular evolution. , 1999, Theoretical population biology.

[8]  Rolf Backofen,et al.  Inferring Noncoding RNA Families and Classes by Means of Genome-Scale Structure-Based Clustering , 2007, PLoS Comput. Biol..

[9]  S. Eddy,et al.  tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. , 1997, Nucleic acids research.

[10]  P. Andolfatto.  Contrasting patterns of X-linked and autosomal nucleotide variation in Drosophila melanogaster and Drosophila simulans. , 2001, Molecular biology and evolution.

[11]  B. Charlesworth,et al.  Effective Population Size and the Faster-X Effect: An Extended Model , 2009, Evolution; international journal of organic evolution.

[12]  W Stephan,et al.  Selection intensity against deleterious mutations in RNA secondary structures and rate of compensatory nucleotide substitutions. , 2001, Genetics.

[13]  T. Ohta.  Slightly Deleterious Mutant Substitutions in Evolution , 1973, Nature.

[14]  L. Duret,et al.  GC Content Evolution of the Human and Mouse Genomes: Insights from the Study of Processed Pseudogenes in Regions of Different Recombination Rates , 2006, Journal of Molecular Evolution.

[15]  P. Andolfatto,et al.  Effective Population Size and the Efficacy of Selection on the X Chromosomes of Two Closely Related Drosophila Species , 2010, Genome biology and evolution.

[16]  Eric Westhof,et al.  Base pairing constraints drive structural epistasis in ribosomal RNA sequences. , 2010, Molecular biology and evolution.

[17]  W. Stephan.  The rate of compensatory evolution. , 1996, Genetics.

[18]  M. Nowak,et al.  Stochastic Tunnels in Evolutionary Dynamics , 2004, Genetics.

[19]  W. Stephan,et al.  Distinctly Different Sex Ratios in African and European Populations of Drosophila melanogaster Inferred From Chromosomewide Single Nucleotide Polymorphism Data , 2007, Genetics.

[20]  D. Weinreich,et al.  Contrasting patterns of nonneutral evolution in proteins encoded in nuclear and mitochondrial genomes. , 2000, Genetics.

[21]  Masatoshi Nei,et al.  Selectionism and neutralism in molecular evolution. , 2005, Molecular biology and evolution.

[22]  M. Lynch.  Scaling expectations for the time to establishment of complex adaptations , 2010, Proceedings of the National Academy of Sciences.

[23]  W. Swanson,et al.  A test for faster X evolution in Drosophila. , 2002, Molecular biology and evolution.

[24]  D. Haussler,et al.  Phylogenetic estimation of context-dependent substitution rates by maximum likelihood. , 2003, Molecular biology and evolution.

[25]  J. Parsch,et al.  On the utility of short intron sequences as a reference for the detection of positive and negative selection in Drosophila. , 2010, Molecular biology and evolution.

[26]  Wolfgang Stephan,et al.  Compensatory evolution of a precursor messenger RNA secondary structure in the Drosophila melanogaster Adh gene , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Daniel J. Gaffney,et al.  Quantifying the slightly deleterious mutation model of molecular evolution. , 2002, Molecular biology and evolution.

[28]  Sean R. Eddy,et al.  Infernal 1.0: inference of RNA alignments , 2009, Bioinform..

[29]  L. Hurst,et al.  Hearing silence: non-neutral evolution at synonymous sites in mammals , 2006, Nature Reviews Genetics.

[30]  T. Ohta,et al.  The Average Number of Generations until Fixation of a Mutant Gene in a Finite Population. , 1969, Genetics.

[31]  J. Parsch,et al.  Comparative sequence analysis and patterns of covariation in RNA secondary structures. , 2000, Genetics.

[32]  E. Westhof,et al.  Geometric nomenclature and classification of RNA base pairs. , 2001, RNA.

[33]  Fyodor A. Kondrashov,et al.  Compensatory evolution in mitochondrial tRNAs navigates valleys of low fitness , 2010, Nature.

[34]  Colin N. Dewey,et al.  Population Genomics: Whole-Genome Analysis of Polymorphism and Divergence in Drosophila simulans , 2007, PLoS biology.

[35]  Sean R. Eddy,et al.  Infernal 1.0: inference of RNA alignments , 2009, Bioinform..

[36]  D. Halligan,et al.  Distributions of selectively constrained sites and deleterious mutation rates in the hominid and murid genomes. , 2010, Molecular biology and evolution.

[37]  Daniel B. Weissman,et al.  The Rate of Fitness-Valley Crossing in Sexual Populations , 2010, Genetics.

[38]  J. Mank,et al.  EFFECTIVE POPULATION SIZE AND THE FASTER‐X EFFECT: EMPIRICAL RESULTS AND THEIR INTERPRETATION , 2010, Evolution; international journal of organic evolution.

[39]  Wolfgang Stephan,et al.  In vivo introduction of unpreferred synonymous codons into the Drosophila Adh gene results in reduced levels of ADH protein. , 2003, Genetics.

[40]  Daniel J. Gaffney,et al.  Effect of the assignment of ancestral CpG state on the estimation of nucleotide substitution rates in mammals , 2008, BMC Evolutionary Biology.

[41]  M. Kimura.  The Neutral Theory of Molecular Evolution: Introduction , 1983 .

[42]  Tom H. Pringle,et al.  The human genome browser at UCSC. , 2002, Genome research.

[43]  Robert D. Finn,et al.  Rfam: updates to the RNA families database , 2008, Nucleic Acids Res..

[44]  T. Connallon.  Adaptive protein evolution of X-linked and autosomal genes in Drosophila: implications for faster-X hypotheses. , 2007, Molecular biology and evolution.

[45]  J. Parsch,et al.  RNA secondary structure and compensatory evolution. , 1999, Genes & genetic systems.

[46]  Scott V Edwards,et al.  SPECIATIONAL HISTORY OF AUSTRALIAN GRASS FINCHES (POEPHILA) INFERRED FROM THIRTY GENE TREES* , 2005, Evolution; international journal of organic evolution.

[47]  J. Baines,et al.  Reduced X-Linked Diversity in Derived Populations of House Mice , 2007, Genetics.

[48]  Robert C. Edgar,et al.  MUSCLE: a multiple sequence alignment method with reduced time and space complexity , 2004, BMC Bioinformatics.

[49]  M. Kimura,et al.  Average time until fixation of a mutant allele in a finite population under continued mutation pressure: Studies by analytical, numerical, and pseudo-sampling methods. , 1980, Proceedings of the National Academy of Sciences of the United States of America.

[50]  Lindell Bromham,et al.  Population size and molecular evolution on islands , 2005, Proceedings of the Royal Society B: Biological Sciences.

[51]  W. Stephan,et al.  Inferring the Demographic History and Rate of Adaptive Substitution in Drosophila , 2006, PLoS genetics.

[52]  Etsuko N. Moriyama,et al.  Codon Usage Bias and tRNA Abundance in Drosophila , 1997, Journal of Molecular Evolution.

[53]  Laurent Duret,et al.  The Impact of Recombination on Nucleotide Substitutions in the Human Genome , 2008, PLoS genetics.

[54]  J. Parsch,et al.  Experimentally Increased Codon Bias in the Drosophila Adh Gene Leads to an Increase in Larval, But Not Adult, Alcohol Dehydrogenase Activity , 2010, Genetics.

[55]  W. Stephan,et al.  Analyzing the evolution of RNA secondary structures in vertebrate introns using Kimura's model of compensatory fitness interactions. , 2008, Molecular biology and evolution.

[56]  David Osumi-Sutherland,et al.  FlyBase: enhancing Drosophila Gene Ontology annotations , 2008, Nucleic Acids Res..

[57]  Gillespie,et al.  Development of Neutral and Nearly Neutral Theories , 1996, Theoretical population biology.

[58]  Elena R. Lozovsky,et al.  Compensatory mutations restore fitness during the evolution of dihydrofolate reductase. , 2010, Molecular biology and evolution.

[59]  William H. Majoros,et al.  A Comparison of Whole-Genome Shotgun-Derived Mouse Chromosome 16 and the Human Genome , 2002, Science.

[60]  M. Takefu,et al.  Compensatory evolution in diploid populations. , 2008, Theoretical population biology.

[61]  M. Kimura.  The role of compensatory neutral mutations in molecular evolution , 1985, Journal of Genetics.

[62]  L. Bromham,et al.  Increased rates of sequence evolution in endosymbiotic bacteria and fungi with small effective population sizes. , 2003, Molecular biology and evolution.

[63]  Casey M. Bergman,et al.  The Evolution of tRNA Genes in Drosophila , 2010, Genome biology and evolution.

[64]  D. A. Kirby,et al.  Maintenance of pre-mRNA secondary structure by epistatic selection. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[65]  Laurent Duret,et al.  Expected Relationship Between the Silent Substitution Rate and the GC Content: Implications for the Evolution of Isochores , 2002, Journal of Molecular Evolution.

[66]  T. Ohta.  Population size and rate of evolution , 1972, Journal of Molecular Evolution.

[67]  Wolfgang Stephan,et al.  Selective constraints in conserved folded RNAs of drosophilid and hominid genomes. , 2011, Molecular biology and evolution.

[68]  Masatoshi Nei,et al.  The neutral theory of molecular evolution in the genomic era. , 2010, Annual review of genomics and human genetics.

[69]  A. Eyre-Walker,et al.  Evidence for Variation in the Effective Population Size of Animal Mitochondrial DNA , 2009, PloS one.

[70]  A. Rich,et al.  Transfer RNA: molecular structure, sequence, and properties. , 1976, Annual review of biochemistry.

[71]  S. Wright,et al.  Sex-Linked Inheritance in Macaque Monkeys: Implications for Effective Population Size and Dispersal to Sulawesi , 2010, Genetics.