Looking for Darwin in Genomic Sequences: Validity and Success Depends on the Relationship Between Model and Data

The role of purifying selection in the origin and maintenance of complex function

Fitness contribution alone should not be the criterion of 'function' in molecular biology and genomics. Disagreement over the use of 'function' in molecular biology and genomics is still with us, almost eight years after publicity surrounding the Encyclopedia of DNA Elements project claimed that 80.4% of the human genome comprises "functional elements". Recent approaches attempt to resolve or reformulate this debate by redefining genomic 'function' in terms of current fitness contribution. In its favour, this redefinition for the genomic context is in apparent conformity with predominant experimental practices, especially in biomedical research, and with ascription of function by selective maintenance. We argue against approaches of this kind, however, on the grounds that they could be seen as non-Darwinian, and fail to properly account for the diversity of non-adaptive processes involved in the origin and maintenance of genomic complexity. We examine cases of molecular and organismal complexity that arise neutrally, showing how purifying selection maintains non-adaptive genomic complexity. Rather than lumping different sorts of genomic complexity together by defining 'function' as fitness contribution, we argue that it is best to separate the heterogeneous contributions of preaptation, exaptation and adaptation to the historical processes of origin and maintenance for complex features.

Improved inference of site-specific positive selection under a generalized parametric codon model when there are multinucleotide mutations and multiple nonsynonymous rates

BACKGROUND

An excess of nonsynonymous substitutions, over neutrality, is considered evidence of positive Darwinian selection. Inference for proteins often relies on estimation of the nonsynonymous to synonymous ratio (ω = dN/dS) within a codon model. However, to ease computational difficulties, ω is typically estimated assuming an idealized substitution process where (i) all nonsynonymous substitutions have the same rate (regardless of impact on organism fitness) and (ii) instantaneous double and triple (DT) nucleotide mutations have zero probability (despite evidence that they can occur). It follows that estimates of ω represent an imperfect summary of the intensity of selection, and that tests based on the ω > 1 threshold could be negatively impacted.

RESULTS

We developed a general-purpose parametric (GPP) modelling framework for codons. This novel approach allows specification of all possible instantaneous codon substitutions, including multiple nonsynonymous rates (MNRs) and instantaneous DT nucleotide changes. Existing codon models are specified as special cases of the GPP model. We use GPP models to implement likelihood ratio tests for ω > 1 that accommodate MNRs and DT mutations. Through both simulation and real data analysis, we find that failure to model MNRs and DT mutations reduces power in some cases and inflates false positives in others. False positives under traditional M2a and M8 models were very sensitive to DT changes. This was exacerbated by the choice of frequency parameterization (GY vs. MG), with rates sometimes > 90% under MG. By including MNRs and DT mutations, accuracy and power was greatly improved under the GPP framework. However, we also find that over-parameterized models can perform less well, and this can contribute to degraded performance of LRTs.

CONCLUSIONS

We suggest GPP models should be used alongside traditional codon models. Further, all codon models should be deployed within an experimental design that includes (i) assessing robustness to model assumptions, and (ii) investigation of non-standard behaviour of MLEs. As the goal of every analysis is to avoid false conclusions, more work is needed on model selection methods that consider both the increase in fit engendered by a model parameter and the degree to which that parameter is affected by un-modelled evolutionary processes.