Abstract
Substitution models of evolution describe the process of genetic variation through fixed mutations and constitute the basis of the evolutionary analysis at the molecular level. Almost 40 years after the development of first substitution models, highly sophisticated, and data-specific substitution models continue emerging with the aim of better mimicking real evolutionary processes. Here I describe current trends in substitution models of DNA, codon and amino acid sequence evolution, including advantages and pitfalls of the most popular models. The perspective concludes that despite the large number of currently available substitution models, further research is required for more realistic modeling, especially for DNA coding and amino acid data. Additionally, the development of more accurate complex models should be coupled with new implementations and improvements of methods and frameworks for substitution model selection and downstream evolutionary analysis.
Introduction
Substitution models of evolution describe the rates of change of fixed mutations among sequences and constitute the basis of the evolutionary analysis of genetic data at the molecular level. Although, the first substitution models that corrected the effects of multiple replacements were documented a long time ago (Jukes and Cantor, 1969; Dayhoff et al., 1978; Kimura, 1980), they were not applied in evolutionary analysis until much later when they were implemented in phylogenetic software packages based on maximum-likelihood (ML) methods (Felsenstein, 1991; Swofford, 1993). Then, it was noted that substitution model misspecification might bias phylogenetic inferences (Posada and Crandall, 2001; Minin et al., 2003; Lemmon and Moriarty, 2004). As a consequence, the selection of the best-fit substitution model became an essential stage in the pipeline of phylogenetic inference (Posada and Crandall, 1998, 2001), which also increased the popularity of substitution models in phylogenetics.
Nowadays, substitution models are routinely used in diverse areas of evolutionary biology. A large number of substitution models exist and new models are emerging to mimic the evolution of particularly complex real datasets. Nevertheless, recent studies suggest that there is still room for improving the commonly used substitution models (e.g., Keane et al., 2006).
This perspective provides an overview on the trends of substitution models of DNA, codon and amino acid evolution, including goals, and pitfalls of well-established substitution models. It also discusses directions for future research in substitution models of evolution and the establishment of complex substitution models as an essential step for phylogenetic inference.
Trends in DNA substitution models
The first and simplest model to mimic the DNA substitution process was described by Jukes and Cantor (JC) (1969). This model considers one rate of change between all nucleotides and equal nucleotide frequencies. However, changes between bases with equal chemical nature (transitions) are more common than changes between bases with different chemical structure (transversions) because the replacement of a similar structure is more likely in terms of molecular energy. Moreover, the genetic code allows for more transitions than transversions without amino acid replacement (Kimura, 1980; Collins and Jukes, 1994). Motivated by this evidence, Kimura (1980) presented the two-parameter model (K80), where the rates of change differ between transitions and transversions. Similarly, Felsenstein (1981) (F81) extended the JC model to include different nucleotide frequencies, which can also appear as a consequence of the physicochemical properties of nucleotides and the operation of natural selection. A number of models were later developed by incorporating extensions to those original models [i.e., HKY (Hasegawa et al., 1985) and SYM (Zharkikh, 1994)]. Following this trend, the most complex model, the general time-reversible model (GTR; Tavaré, 1986), incorporates different rates for every change and different nucleotide frequencies. In addition, a proportion of invariable sites (+I) (Shoemaker and Fitch, 1989) and/or rate of variation across sites (+G) (Yang, 1994) can be incorporated into any model. For technical details about these parametric DNA substitution models the reader is referred to the following reviews (Liò and Goldman, 1998; Felsenstein, 2004).
The stationary, reversible and homogeneous DNA substitution models derived from all possible combinations (equal/different rates of change and frequencies, with/without +G and/or +I; more than 1600 models) have already been defined, and are currently implemented in some phylogenetic programs (e.g., Darriba et al., 2012; Arenas and Posada, 2014b), see also Table 1. However, despite the large number of available DNA substitution models, the GTR +G +I model usually fits real data better than the other (simpler) alternative models (Sumner et al., 2012b), suggesting that the evolutionary process is very complex. Importantly, the GTR model presents undesirable mathematical properties [i.e., the multiplication of two GTR matrices does not return another GTR matrix (Sumner et al., 2012a)] for its application as a Markov model in phylogenetics (see Gatto et al., 2006; Sumner et al., 2012a,b). Indeed, the frequent selection of the most complex substitution model may suggest that more complex models could improve the fitting to real data (e.g., Jayaswal et al., 2011).
Table 1
| Program | Task | Substitution model | Genome | Rate variation | References | ||
|---|---|---|---|---|---|---|---|
| Non-coding | Coding | Amino acid | |||||
| PhyML | Inference | All | – | Blosum62, CpRev, Dayhoff, FLU, HIVb, HIVw, JTT, LG, Mtart, Mtmam, Mtrev, RtRev, VT, WAG +F | – | +I +G | Guindon et al., 2010 |
| CodonPhyML | Inference | – | GY94b, MG94, YAP, ECMs | – | – | +G | Gil et al., 2013 |
| RAxML | Inference | JC, K80, HKY, GTR | Nta | Blosum62, CpRev, Dayhoff, DUMMY, FLU, HIVb, HIVw, JTT, JonesDCMUT, LG, Mtart, Mtmam, Mtrev, Mtzoa, PMB, RtRev, STMREV, VT, WAG +F | Partitioned models can be specified | +I +G | Stamatakis, 2006 |
| MEGA | Inference | JC, K2P, HKY, TN93, GTR | NG86 | CpRev, Dayhoff, JTT, LG, MtRev, RtRev, WAG | – | +I +G | Tamura et al., 2013 |
| Hyphy | Inference | All | GY94b, MG94, ECM | Dayhoff, HIVb, HIVw, JTT, Mtmam, Mtrev, RtRev, WAG +F | Partitioned models can be specified | +I +G | Pond et al., 2005 |
| PAML | Inference | All | GY94b,c, ECMs | CpRev, Dayhoff, DayhoffDCMUT, Grantham, JTT, JonesDCMUT, LG, Miyata, Mtart, Mtmam, Mtrev24, Mtzoa, WAG +F | – | +I +G | Yang, 2007 |
| MrBayes and BEST | Inference | All | GY94b, MG94 | Blosum62, CpRev, Dayhoff, Mtmam, Mtrev, RtRev, VT, WAG +F | Partitioned models can be specified | +I +G | Ronquist et al., 2012 |
| BEAST | Inference | JC, HKY, TN93, GTR | Nta | Blosum62, CpRev, Dayhoff, FLU, JTT, LG, Mtrev, WAG | – | +I +G | Bouckaert et al., 2014 |
| OmegaMap | Inference | – | NY98b | – | – | – | Wilson and McVean, 2006 |
| Lamarc | Inference | JC, K2P, F84, GTR | – | – | – | +G | Kuhner, 2006 |
| CodABC | Inference | – | GY94 | – | – | +I +G | Arenas et al., 2015a |
| MySSP | Simulation | All | – | – | – | +G | Rosenberg, 2005 |
| Seq-Gen | Simulation | All | Nta | Blosum62, CpRev, JTT, mtREV, PAM, and WAG +F | – | +I +G | Rambaut and Grassly, 1997 |
| indel-Seq-Gen | Simulation | All | Nta | Blosum62, CpRev, JTT, mtREV, PAM, WAG +F | – | +I +G | Strope et al., 2009 |
| INDELible | Simulation | All | GY94b,c, ECMs | Blosum62, CpRev, Dayhoff, DayhoffDCMUT, HIVb, HIVw, JTT, JonesDCMUT, LG, Mtart, Mtmam, Mtrev24, RtRev, VT, WAG +F | – | +I +G | Fletcher and Yang, 2009 |
| EVOLVER | Simulation | All | GY94b,c, ECMs | CpRev, Dayhoff, DayhoffDCMUT, Grantham, JTT, JonesDCMUT, LG, Miyata, Mtart, Mtmam, Mtrev24, Mtzoa, WAG +F | – | +I +G | Yang, 2007 |
| Recodon and NetRecodon | Simulation | All | GY94b | – | – | +I +G | Arenas and Posada, 2007, 2010 |
| ProteinEvolver | Simulation | All | – | 3D structural constraints, Blosum62, CpRev, Dayhoff, DayhoffDCMUT, HIVb, HIVw, JTT, JonesDCMUT, LG, Mtart, Mtmam, Mtrev24, RtRev, VT, WAG +F | – | +I +G | Arenas et al., 2013 |
| CoalEvol | Simulation | All | GY94b,c, MG94, ECMs, HB | Blosum62, CpRev, Dayhoff, DayhoffDCMUT, HIVb, HIVw, JTT, JonesDCMUT, LG, Mtart, Mtmam, Mtrev24, RtRev, VT, WAG +F | – | +I +G | Arenas and Posada, 2014b |
| SIMPROT | Simulation | – | – | PAM, JTT, PMB | – | +G | Pang et al., 2005 |
| PhyloSim | Simulation | All | GY94b,c, ECMs | CpRev, JTT, JonesDCMUT, LG, Mtart, Mtmam, Mtrev24, Mtzoa, WAG +F | – | +I +G | Sipos et al., 2011 |
| π BUSS | Simulation | HKY, TN93, GTR | GY94, MG94 | BLOSUM, Dayhoff, LG, JTT, WAG +F | – | +I +G | Bielejec et al., 2014 |
| SGWE | Simulation | All | GY94b,c, MG94, ECMs, HB | Blosum62, CpRev, Dayhoff, DayhoffDCMUT, HIVb, HIVw, JTT, JonesDCMUT, LG, Mtart, Mtmam, Mtrev24, RtRev, VT, WAG +F | User-specified regions | +I +G | Arenas and Posada, 2014b |
| ALF | Simulation | F84, HKY, TN93, GTR | GY94b, ECMs | Gonnet, JTT, LG, PAM, WAG +F | User-specified regions | +I +G | Dalquen et al., 2012 |
Substitution models implemented in relevant software for phylogenetic inference.
“Task” indicates if the program is oriented to perform evolutionary inference or simulation of molecular evolution. “Substitution model” indicates the implemented models of non-coding DNA evolution [“All” means that all the reversible nucleotide substitution models are considered, JC, …, GTR], codon models [“ECMs” indicates empirical codon models, MG94 model refers to Muse and Gaut (1994) and HB model refers to Halpern and Bruno (1998)] and amino acid models [“+F” indicates that amino acid frequencies can be modeled]. “Rate variation” includes proportion of invariable sites “+I” and substitution rate heterogeneity across sites according to a gamma distribution “+G.” “Genome” indicates the consideration of genome evolution through region-specific substitution models.
Coding sequences are simulated by nucleotide substitution models, avoiding stop codons.
dN/dS can vary across codons.
dN/dS can vary across branches.
In order to improve the fitting to real data, current trends in the development of DNA substitution models involve non-reversible (asymmetric) and non-stationary (nucleotide composition can change over time) matrices (e.g., Boussau and Gouy, 2006; Jayaswal et al., 2011), or even consider neighbor interactions (Lunter and Hein, 2004), that can lead to more accurate phylogenetic inferences (Boussau and Gouy, 2006; Kaehler et al., 2015). However, these models have not been implemented yet in the most popular phylogenetic software packages due to their implicit complexity. In addition to the proposal of novel substitution models, future research should also address the implementation of these complex models in methods and programs for substitution model choice and phylogenetic inference.
Trends in codon substitution models
Sites within codons evolve at different rates and, consequently, they should not be equally treated (Shapiro et al., 2006). Indeed, considering molecular evolution at the codon level allows us to incorporate more realistic evolutionary patterns for each codon position. Actually, evolutionary inferences based on the modeling of codon evolution are more robust than those derived from empirical amino acid models (Benner et al., 1994; Seo and Kishino, 2008).
Interestingly, mutations at the codon level can be classified as synonymous (silent) and non-synonymous (amino-acid replacing), which provides a measure of selective pressure at the molecular level (molecular adaptation). Then, the first codon substitution models included the non-synonymous and synonymous substitution rates, dN and dS respectively (Muse and Gaut, 1994), or the dN/dS ratio (Goldman and Yang, 1994). Despite the large number of (frequently ignored) considerations that should be made for the analysis and interpretation of dN/dS [i.e., dN/dS at the codon level can be biased if recombination is ignored (Anisimova et al., 2003; Arenas and Posada, 2010, 2014a) or if nucleotide frequencies vary across sites (Arenas and Posada, 2014b) and dN/dS should be estimated from samples of different populations (Kryazhimskiy and Plotkin, 2008; Pellissier, 2015)], dN/dS is commonly estimated in evolutionary biology for testing hypothesis related to selective pressure (e.g., Yang and Bielawski, 2000; Perez-Losada et al., 2009, 2011; Lopes et al., 2014; Arenas, 2015b; Arenas et al., 2015b; Lopez-Bueno et al., 2015). dN/dS >1 indicates that substitutions in the protein-coding gene were enriched for those that altered the amino acids states, suggesting diversifying (positive) selection. By contrast, dN/dS < 1 and dN/dS = 1 can be interpreted as purifying (negative) selection and neutral evolution, respectively.
Additional codon models have been proposed in order to better fit particular codon datasets. The GY94 codon model (Goldman and Yang, 1994) was extended to consider different nucleotide models (e.g., GTR; Kosakovsky Pond and Frost, 2005; Pond and Frost, 2005; Pond and Muse, 2005; Arenas and Posada, 2014b), dN/dS variation across sites (Yang et al., 2000; Anisimova et al., 2001; Kosakovsky Pond and Frost, 2005; Yang, 2007) and across branches (e.g., Pond and Frost, 2005; Yang, 2007; Dutheil et al., 2012). Other codon models implement information about the physicochemical properties of the encoded amino acids (Wong et al., 2006). Additionally, several empirical codon models—based on large databases of coding data—have been proposed (Schneider et al., 2005; Doron-Faigenboim and Pupko, 2007; Kosiol et al., 2007) and also codon models that consider codon bias (McVean and Charlesworth, 1999; Nielsen et al., 2007; Yang and Nielsen, 2008) or the effects of different GC contents (Misawa, 2011). Here, a promising trend is the emergence of mutation-selection models (e.g., Halpern and Bruno, 1998; Yang and Nielsen, 2008; Rodrigue et al., 2010) that attempt to integrate complex selection patterns into the mutation process and that outperform simple neutral models (Lawrie et al., 2011). Note that the interplay of mutational biases and weak selection can be complex, where constrained sites can evolve faster than neutral sites (McVean and Charlesworth, 1999; Lawrie et al., 2011). The consideration of these effects in substitution models of evolution can improve the identification of functional regions and the fitting to the data (Lawrie et al., 2011). For technical aspects on codon models the reader is referred to Anisimova and Kosiol (2009) and Cannarozzi and Schneider (2012).
Therefore, codon models allow us to perform accurate evolutionary analysis and to explore signatures of molecular adaptation. However, a problematic technical aspect regarding the implementation of these models is their large exchangeability matrices (61 × 61, note that stop codons are excluded). As a consequence, large amounts of data are needed to generate statistically well-supported empirical codon matrices and the computational burden is heavy. Fortunately, research on codon-based algorithm optimization is already generating new evolutionary tools to simulate (e.g., Fletcher and Yang, 2009; Arenas, 2012; Arenas and Posada, 2012) and analyze (e.g., Pond et al., 2005; Gil et al., 2013; Arenas et al., 2015a; Zoller et al., 2015) codon evolution (see Table 1), although further work is still required in this regard (e.g., the implementation of parallel computing of probability matrices).
The future of codon models will probably be related to the development of more complex models to better fit real data. In this concern, in addition to the development of new empirical models, codon models may follow two interesting trends. First, the consideration of heterogeneity along the sequence and over time, where different sites/regions and time periods could evolve under different models (Arenas, 2015a; Zoller et al., 2015). Note that these partition schemes can be very realistic, for example by considering different models for coding and non-coding regions. Moreover, it is known that codon models based on different codon frequencies across sites can bias dN/dS estimates (Arenas and Posada, 2014b), and Zoller et al. (2015) recently found that a two-partition codon model resolves a phylogeny better than a one-partition codon model. Consequently, methods and software to identify the best codon substitution model for particular codon regions (i.e., following Bao et al., 2008; Delport et al., 2010) and time periods are demanded. A second trend may be the consideration of structural information of proteins in codon models. Information derived from the protein function and from the protein folding stability [i.e., considering energy functions (Grahnen et al., 2011; Liberles et al., 2012; Arenas et al., 2015c)] of the encoded proteins could be considered in codon models. However, these implementations would require large computational costs if the protein structure varies with time or if more than one protein structure is needed to represent the encoded proteins of the dataset.
Trends in amino acid substitution models
Substitution models of amino acid evolution intend to mimic the evolution of protein data, which is crucial for testing a variety of hypotheses such as selection toward novel proteins (e.g., Fares et al., 2002), rate of protein evolution (e.g., Alvarez-Ponce and Fares, 2012), role of protein function on protein evolution (e.g., Liberles et al., 2011), evolutionary aspects of protein-protein interaction networks (e.g., Alvarez-Ponce and Fares, 2012), multiple sequence alignment based on protein evolution (Zhao and Sacan, 2015), or phylogenetic tree and ancestral protein reconstructions (e.g., Perez-Jimenez et al., 2011). Substitution models of amino acid evolution can be classified in two major groups: (i) empirical models—based on large protein databases—[e.g., CpRev (Adachi et al., 2000), Dayhoff (Dayhoff et al., 1978), DayhoffDCMUT (Kosiol and Goldman, 2005), HIVb (Nickle et al., 2007), HIVw (Nickle et al., 2007), JTT (Jones et al., 1992), JonesDCMUT (Kosiol and Goldman, 2005), LG (Le and Gascuel, 2008), Mtart (Abascal et al., 2007), Mtmam (Yang et al., 1998), Mtrev24 (Adachi and Hasegawa, 1996), RtRev (Dimmic et al., 2002), VT (Müller and Vingron, 2000), WAG (Whelan and Goldman, 2001), see Table 1 for their implementation in phylogenetic software], and (ii) parametric models—based on parameters that describe protein evolution—(Rastogi et al., 2006; Liberles et al., 2012).
An empirical amino acid substitution model consists of a 20 × 20 matrix of exchangeability rates and 20 amino acid frequencies. This simplicity (compared to other amino acid models such as those based on structural constraints) leads to advantages but also pitfalls. Empirical models can be incorporated into the commonly used likelihood functions implemented in the standard phylogenetic software by assuming site-independence (all sites evolve under the same model). Heterogeneous evolution can also be modeled by specifying different empirical models for different partitions (i.e., sites or domains) of the protein sequence (Halpern and Bruno, 1998; Lartillot and Philippe, 2004; Zoller and Schneider, 2013). However, a given real dataset may not be properly represented by any of the currently available empirical models. For example, Keane et al. (2006) found that the best fitting empirical model for two large proteobacteria and archaea protein datasets was originally derived from retroviral Pol proteins.
In order to provide an alternative to empirical models, constraints on the protein folding have been considered to generate parametric amino acid substitution models that have led to significant improvements (with respect to empirical models) when fitting real data (e.g., Taverna and Goldstein, 2000; Parisi and Echave, 2005; Goldstein, 2011; Grahnen et al., 2011; Wilke, 2012; Arenas et al., 2013, 2015c; Bordner and Mittelmann, 2013). However, these models are not well-established yet in the evolutionary analysis of protein data because the commonly used evolutionary frameworks implement likelihood functions that cannot deal with site-dependence. Consequently, a current trend in structurally constrained substitution models is to generate site-specific matrices that can be incorporated into common phylogenetic frameworks (Arenas et al., 2015c).
New empirical models are emerging to represent protein families (e.g., Cox and Foster, 2013); similarly, new parametric models are appearing to account for different evolutionary processes across sites and over time—covarion models—(Usmanova et al., 2015), epistatic fitness landscapes (Usmanova et al., 2015), or complex structural constraints (Bordner and Mittelmann, 2013; Arenas et al., 2015c). Nevertheless, the large variety and complexity of current amino acid substitution models cause problems when implemented. As noted, complex amino acid models (such as those based on protein folding stability) are not established yet in popular phylogenetic software due to their implicit complexity. In addition, the user often restricts the candidate substitution models to those empirical models implemented in common substitution model choice programs (e.g., Abascal et al., 2005), which may lead to severe incongruences (Keane et al., 2006).
Therefore, although current research on amino acid substitution models is providing more sophisticated models, these models are usually not applied by evolutionary biologists because they are not implemented in evolutionary frameworks and, consequently, these models are often forgotten. In order to consider complex substitution models in model selection and in evolutionary analysis, an alternative strategy can be the approximate Bayesian computation (ABC) approach (Beaumont, 2010; Csilléry et al., 2010; Sunnaker et al., 2013; Lopes et al., 2014; Arenas, 2015a). Basically, simulated data under different complex models are contrasted with real data through multiple regression adjustments to identify the model that best fits the real data, and to estimate the parameter values of the model corresponding to the studied dataset (see Lopes et al., 2014 for an example of ABC using complex codon models).
Conclusions
The modeling of genome evolution is nowadays important in population genomics and phylogenomics (Kumar et al., 2012; Librado et al., 2014). In that regard, as noted above, a variety of substitution models are known to mimic the evolution of coding and non-coding data. Importantly, the evolutionary process may differ between genomic regions because molecular evolution is often highly heterogeneous (Shapiro et al., 2006; Bofkin and Goldman, 2007; Arbiza et al., 2011) and, consequently, genome evolution should be mimicked with specific substitution models for each genomic region (Arbiza et al., 2011; Dalquen et al., 2012; Arenas and Posada, 2014b), see software implementation of partition models in Table 1. However, there is a need for methods to identify the regions that can better fit with a single substitution model. A strategy to perform this task could emulate genetic algorithms for the detection of recombination breakpoints based on phylogenetic tree incongruence (Fitch and Goodman, 1991; Grassly and Holmes, 1997; Kosakovsky Pond et al., 2006) assuming that different substitution models can also lead to significant signatures of phylogenetic tree incongruence (Minin et al., 2003; Lemmon and Moriarty, 2004). Then, these partitioning schemes could be evaluated with tools such as PartitionFinder (Lanfear et al., 2012). Additionally, a realistic modeling may involve not only region-specific models, but also branch-specific models (Ho, 2009). However, in practice only large datasets could provide sufficient statistical support for identifying the best models at those levels.
It seems clear that future research on substitution models of evolution will involve the development of more sophisticated and realistic substitution models. For example, the increasing amount of genomic data will probably require, in addition to partitioning, the development of complex mixture models where each site/region can be modeled with more than one substitution model. Concomitantly, efforts are needed to design and implement more robust methods to evaluate, compare and apply these complex substitution models.
Funding
This work was supported by the Portuguese Government through the FCT Starting Grant IF/00955/2014.
Conflict of interest statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Acknowledgments
I want to thank Frontiers in Evolutionary and Population Genetics for the invitation to contribute with this inaugural perspective article derived from my new role in the journal as Associate Editor. I thank Marcos Pérez-Losada, at The George Washington University (USA), for helpful comments. I also thank three reviewers for insightful comments.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AbascalF.PosadaD.ZardoyaR. (2007). MtArt: a new model of amino acid replacement for arthropoda. Mol. Biol. Evol.24, 1–5. 10.1093/molbev/msl136
2
AbascalF.ZardoyaR.PosadaD. (2005). ProtTest: selection of best-fit models of protein evolution. Bioinformatics21, 2104–2105. 10.1093/bioinformatics/bti263
3
AdachiJ.HasegawaM. (1996). MOLPHY version 2.3: programs for molecular phylogenetics based in maximum likelihood. Comp. Sci. Monogr.28, 1–150.
4
AdachiJ.WaddellP. J.MartinW.HasegawaM. (2000). Plastid genome phylogeny and a model of amino acid substitution for proteins encoded by chloroplast DNA. J. Mol. Evol.50, 348–358. 10.1007/s002399910038
5
Alvarez-PonceD.FaresM. A. (2012). Evolutionary rate and duplicability in the Arabidopsis thaliana protein-protein interaction network. Genome Biol. Evol.4, 1263–1274. 10.1093/gbe/evs101
6
AnisimovaM.BielawskiJ. P.YangZ. (2001). Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution. Mol. Biol. Evol.18, 1585–1592. 10.1093/oxfordjournals.molbev.a003945
7
AnisimovaM.KosiolC. (2009). Investigating protein-coding sequence evolution with probabilistic codon substitution models. Mol. Biol. Evol.26, 255–271. 10.1093/molbev/msn232
8
AnisimovaM.NielsenR.YangZ. (2003). Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites. Genetics164, 1229–1236. Available online at: http://www.genetics.org/content/164/3/1229.long
9
ArbizaL.PatricioM.DopazoH.PosadaD. (2011). Genome-wide heterogeneity of nucleotide substitution model fit. Genome Biol. Evol.3, 896–908. 10.1093/gbe/evr080
10
ArenasM.Dos SantosH. G.PosadaD.BastollaU. (2013). Protein evolution along phylogenetic histories under structurally constrained substitution models. Bioinformatics29, 3020–3028. 10.1093/bioinformatics/btt530
11
ArenasM.LopesJ. S.BeaumontM. A.PosadaD. (2015a). CodABC: a computational framework to coestimate recombination, substitution, and molecular adaptation rates by approximate bayesian computation. Mol. Biol. Evol.32, 1109–1112. 10.1093/molbev/msu411
12
ArenasM.Lorenzo-RedondoR.Lopez-GalindezC. (2015b). Influence of mutation and recombination on HIV-1 in vitro fitness recovery. Mol. Phylogenet. Evol.94, 264–270. 10.1016/j.ympev.2015.09.001
13
ArenasM.PosadaD. (2007). Recodon: coalescent simulation of coding DNA sequences with recombination, migration and demography. BMC Bioinformatics8:458. 10.1186/1471-2105-8-458
14
ArenasM.PosadaD. (2010). Coalescent simulation of intracodon recombination. Genetics184, 429–437. 10.1534/genetics.109.109736
15
ArenasM.PosadaD. (2012). Simulation of coding sequence evolution, in Codon Evolution, eds CannarozziG. M.SchneiderA. (Oxford: Oxford University Press), 126–132.
16
ArenasM.PosadaD. (2014a). The influence of recombination on the estimation of selection from coding sequence alignments, in Natural Selection: Methods and Applications, ed FaresM. A. (Boca Raton, FL: CRC Press/Taylor & Francis), 112–125.
17
ArenasM.PosadaD. (2014b). Simulation of genome-wide evolution under heterogeneous substitution models and complex multispecies coalescent Histories. Mol. Biol. Evol.31, 1295–1301. 10.1093/molbev/msu078
18
ArenasM.Sánchez-CobosA.BastollaU. (2015c). Maximum likelihood phylogenetic inference with selection on protein folding stability. Mol. Biol. Evol.32, 2195–2207. 10.1093/molbev/msv085
19
ArenasM. (2012). Simulation of molecular data under diverse evolutionary scenarios. PLoS Comput. Biol.8:e1002495. 10.1371/journal.pcbi.1002495
20
ArenasM. (2015a). Advances in computer simulation of genome evolution: toward more realistic evolutionary genomics analysis by approximate bayesian computation. J. Mol. Evol.80, 189–192. 10.1007/s00239-015-9673-0
21
ArenasM. (2015b). Genetic consequences of antiviral therapy on HIV-1. Comput. Math. Methods Med.2015, 9. 10.1155/2015/395826
22
BaoL.GuH.DunnK. A.BielawskiJ. P. (2008). Likelihood-based clustering (LiBaC) for codon models, a method for grouping sites according to similarities in the underlying process of evolution. Mol. Biol. Evol.25, 1995–2007. 10.1093/molbev/msn145
23
BeaumontM. A. (2010). Approximate Bayesian computation in evolution and ecology. Annu. Rev. Ecol. Evol. Syst.41, 379–405. 10.1146/annurev-ecolsys-102209-144621
24
BennerS. A.CohenM. A.GonnetG. H. (1994). Amino acid substitution during functionally constrained divergent evolution of protein sequences. Protein Eng.7, 1323–1332. 10.1093/protein/7.11.1323
25
BielejecF.LemeyP.CarvalhoL. M.BaeleG.RambautA.SuchardM. A. (2014). piBUSS: a parallel BEAST/BEAGLE utility for sequence simulation under complex evolutionary scenarios. BMC Bioinformatics15:133. 10.1186/1471-2105-15-133
26
BofkinL.GoldmanN. (2007). Variation in evolutionary processes at different codon positions. Mol. Biol. Evol.24, 513–521. 10.1093/molbev/msl178
27
BordnerA. J.MittelmannH. D. (2013). A new formulation of protein evolutionary models that account for structural constraints. Mol. Biol. Evol.31, 736–749. 10.1093/molbev/mst240
28
BouckaertR.HeledJ.KühnertD.VaughanT.WuC. H.XieD.et al. (2014). BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput. Biol.10:e1003537. 10.1371/journal.pcbi.1003537
29
BoussauB.GouyM. (2006). Efficient likelihood computations with nonreversible models of evolution. Syst. Biol.55, 756–768. 10.1080/10635150600975218
30
CannarozziG. M.SchneiderA. (2012). Codon Evolution.Oxford: Oxford University Press.
31
CollinsD. W.JukesT. H. (1994). Rates of transition and transversion in coding sequences since the human-rodent divergence. Genomics20, 386–396. 10.1006/geno.1994.1192
32
CoxC. J.FosterP. G. (2013). A 20-state empirical amino-acid substitution model for green plant chloroplasts. Mol. Phylogenet. Evol.68, 218–220. 10.1016/j.ympev.2013.03.030
33
CsilléryK.BlumM. G.GaggiottiO. E.FrançoisO. (2010). Approximate Bayesian Computation (ABC) in practice. Trends Ecol. Evol.25, 410–418. 10.1016/j.tree.2010.04.001
34
DalquenD. A.AnisimovaM.GonnetG. H.DessimozC. (2012). ALF–a simulation framework for genome evolution. Mol. Biol. Evol.29, 1115–1123. 10.1093/molbev/msr268
35
DarribaD.TaboadaG. L.DoalloR.PosadaD. (2012). jModelTest 2: more models, new heuristics and parallel computing. Nat. Methods9, 772. 10.1038/nmeth.2109
36
DayhoffM. O.SchwartzR. M.OrcuttB. C. (1978). A model of evolutionary change in proteins, in Atlas of Protein Sequence and Structure, ed DayhoffM. O. (Washington, DC: National Biomedical Research Foundation), 345–352.
37
DelportW.SchefflerK.BothaG.GravenorM. B.MuseS. V.Kosakovsky PondS. L. (2010). CodonTest: modeling amino acid substitution preferences in coding sequences. PLoS Comput Biol6:e1000885. 10.1371/journal.pcbi.1000885
38
DimmicM. W.RestJ. S.MindellD. P.GoldsteinR. A. (2002). rtREV: an amino acid substitution matrix for inference of retrovirus and reverse transcriptase phylogeny. J. Mol. Evol.55, 65–73. 10.1007/s00239-001-2304-y
39
Doron-FaigenboimA.PupkoT. (2007). A combined empirical and mechanistic codon model. Mol. Biol. Evol.24, 388–397. 10.1093/molbev/msl175
40
DutheilJ. Y.GaltierN.RomiguierJ.DouzeryE. J.RanwezV.BoussauB. (2012). Efficient selection of branch-specific models of sequence evolution. Mol. Biol. Evol.29, 1861–1874. 10.1093/molbev/mss059
41
FaresM. A.BarrioE.Sabater-MuñozB.MoyaA. (2002). The evolution of the heat-shock protein GroEL from Buchnera, the primary endosymbiont of aphids, is governed by positive selection. Mol. Biol. Evol.19, 1162–1170. 10.1093/oxfordjournals.molbev.a004174
42
FelsensteinJ. (1981). Evolutionary trees from DNA sequences: a maximum likelihood approach. J. Mol. Evol.17, 368–376. 10.1007/BF01734359
43
FelsensteinJ. (1991). PHYLIP: Phylogenetic Inference Package, 3.4 Edn. Seattle, WA: University of Washington.
44
FelsensteinJ. (2004). Inferring Phylogenies.Sunderland, MA: Sinauer Associates.
45
FitchD. H. A.GoodmanM. (1991). Phylogenetic scanning: a computer assisted algorithm for mapping gene conversions and othe recombinational events. CABIOS7, 207–215. 10.1093/bioinformatics/7.2.207
46
FletcherW.YangZ. (2009). INDELible: a flexible simulator of biological sequence evolution. Mol. Biol. Evol.26, 1879–1888. 10.1093/molbev/msp098
47
GattoL.CatanzaroD.MilinkovitchM. C. (2006). Assessing the applicability of the GTR nucleotide substitution model through simulations. Evol. Bioinform. Online2, 145–155. Available online at: http://www.la-press.com/assessing-the-applicability-of-the-gtr-nucleotide-substitution-model-t-article-a145
48
GilM.ZanettiM. S.ZollerS.AnisimovaM. (2013). CodonPhyML: fast maximum likelihood phylogeny estimation under codon substitution models. Mol. Biol. Evol.30, 1270–1280. 10.1093/molbev/mst034
49
GoldmanN.YangZ. (1994). A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol. Biol. Evol.11, 725–736.
50
GoldsteinR. A. (2011). The evolution and evolutionary consequences of marginal thermostability in proteins. Proteins79, 1396–1407. 10.1002/prot.22964
51
GrahnenJ. A.NandakumarP.KubelkaJ.LiberlesD. A. (2011). Biophysical and structural considerations for protein sequence evolution. BMC Evol. Biol.11:361. 10.1186/1471-2148-11-361
52
GrasslyN. C.HolmesE. C. (1997). A likelihood method for the detection of selection and recombination using nucleotide sequences. Mol. Biol. Evol.14, 239–247. 10.1093/oxfordjournals.molbev.a025760
53
GuindonS.DufayardJ. F.LefortV.AnisimovaM.HordijkW.GascuelO. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol.59, 307–321. 10.1093/sysbio/syq010
54
HalpernA. L.BrunoW. J. (1998). Evolutionary distances for protein-coding sequences: modeling site-specific residue frequencies. Mol. Biol. Evol.15, 910–917. 10.1093/oxfordjournals.molbev.a025995
55
HasegawaM.KishinoH.YanoT. (1985). Dating the human-ape splitting by a molecular clock of mitochondrial DNA. J. Mol. Evol.22, 160–174. 10.1007/BF02101694
56
HoS. Y. (2009). An examination of phylogenetic models of substitution rate variation among lineages. Biol. Lett.5, 421–424. 10.1098/rsbl.2008.0729
57
JayaswalV.JermiinL. S.PoladianL.RobinsonJ. (2011). Two stationary nonhomogeneous Markov models of nucleotide sequence evolution. Syst. Biol.60, 74–86. 10.1093/sysbio/syq076
58
JonesD. T.TaylorW. R.ThorntonJ. M. (1992). The rapid generation of mutation data matrices from protein sequences. Comput. Appl. Biosci.8, 275–282. 10.1093/bioinformatics/8.3.275
59
JukesT. H.CantorC. R. (1969). Evolution of protein molecules, in Mammalian Protein Metabolism, ed MunroH. M. (New York, NY: Academic Press), 21–132.
60
KaehlerB. D.YapV. B.ZhangR.HuttleyG. A. (2015). Genetic distance for a general non-stationary markov substitution process. Syst. Biol.64, 281–293. 10.1093/sysbio/syu106
61
KeaneT. M.CreeveyC. J.PentonyM. M.NaughtonT. J.MclnerneyJ. O. (2006). Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol. Biol.6:29. 10.1186/1471-2148-6-29
62
KimuraM. (1980). A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J. Mol. Evol.16, 111–120. 10.1007/BF01731581
63
Kosakovsky PondS. L.FrostS. D. (2005). Not so different after all: a comparison of methods for detecting amino Acid sites under selection. Mol. Biol. Evol.22, 1208–1222. 10.1093/molbev/msi105
64
Kosakovsky PondS. L.PosadaD.GravenorM. B.WoelkC. H.FrostS. D. (2006). Automated phylogenetic detection of recombination using a genetic algorithm. Mol. Biol. Evol.23, 1891–1901. 10.1093/molbev/msl051
65
KosiolC.GoldmanN. (2005). Different versions of the Dayhoff rate matrix. Mol. Biol. Evol.22, 193–199. 10.1093/molbev/msi005
66
KosiolC.HolmesI.GoldmanN. (2007). An empirical codon model for protein sequence evolution. Mol. Biol. Evol.24, 1464–1479. 10.1093/molbev/msm064
67
KryazhimskiyS.PlotkinJ. B. (2008). The population genetics of dN/dS. PLoS Genet.4:e1000304. 10.1371/journal.pgen.1000304
68
KuhnerM. K. (2006). LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters. Bioinformatics22, 768–770. 10.1093/bioinformatics/btk051
69
KumarS.FilipskiA. J.BattistuzziF. U.Kosakovsky PondS. L.TamuraK. (2012). Statistics and truth in phylogenomics. Mol. Biol. Evol.29, 457–472. 10.1093/molbev/msr202
70
LanfearR.CalcottB.HoS. Y.GuindonS. (2012). Partitionfinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol.29, 1695–1701. 10.1093/molbev/mss020
71
LartillotN.PhilippeH. (2004). A Bayesian mixture model for across-site heterogeneities in the amino-acid replacement process. Mol. Biol. Evol.21, 1095–1109. 10.1093/molbev/msh112
72
LawrieD. S.PetrovD. A.MesserP. W. (2011). Faster than neutral evolution of constrained sequences: the complex interplay of mutational biases and weak selection. Genome Biol. Evol.3, 383–395. 10.1093/gbe/evr032
73
LeS. Q.GascuelO. (2008). An improved general amino acid replacement matrix. Mol. Biol. Evol.25, 1307–1320. 10.1093/molbev/msn067
74
LemmonA. R.MoriartyE. C. (2004). The importance of proper model assumption in bayesian phylogenetics. Syst. Biol.53, 265–277. 10.1080/10635150490423520
75
LiòP.GoldmanN. (1998). Models of molecular evolution and phylogeny. Genome Res.8, 1233–1244.
76
LiberlesD. A.TeichmannS. A.BaharI.BastollaU.BloomJ.Bornberg-BauerE.et al. (2012). The interface of protein structure, protein biophysics, and molecular evolution. Protein Sci.21, 769–785. 10.1002/pro.2071
77
LiberlesD. A.TisdellM. D.GrahnenJ. A. (2011). Binding constraints on the evolution of enzymes and signalling proteins: the important role of negative pleiotropy. Proc. Biol. Sci.278, 1930–1935. 10.1098/rspb.2010.2637
78
LibradoP.VieiraF. G.Sánchez-GraciaA.KolokotronisS. O.RozasJ. (2014). Mycobacterial phylogenomics: an enhanced method for gene turnover analysis reveals uneven levels of gene gain and loss among species and gene families. Genome Biol. Evol.6, 1454–1465. 10.1093/gbe/evu117
79
LopesJ. S.ArenasM.PosadaD.BeaumontM. A. (2014). Coestimation of recombination, substitution and molecular adaptation rates by approximate Bayesian computation. Heredity112, 255–264. 10.1038/hdy.2013.101
80
López-BuenoA.RastrojoA.PeiroR.ArenasM.AlcamiA. (2015). Ecological connectivity shapes quasispecies structure of RNA viruses in an Antarctic lake. Mol. Ecol.24, 4812–4825. 10.1111/mec.13321
81
LunterG.HeinJ. (2004). A nucleotide substitution model with nearest-neighbour interactions.Bioinformatics20(Suppl. 1), i216–i223. 10.1093/bioinformatics/bth901
82
McVeanG. A.CharlesworthB. (1999). A population genetic model for the evolution of synonymous codon usage: patterns and predictions. Genet. Res. Camb.74, 145–158. 10.1017/S0016672399003912
83
MininV.AbdoZ.JoyceP.SullivanJ. (2003). Performance-based selection of likelihood models for phylogeny estimation. Syst. Biol.52, 674–683. 10.1080/10635150390235494
84
MisawaK. (2011). A codon substitution model that incorporates the effect of the GC contents, the gene density and the density of CpG islands of human chromosomes. BMC Genomics12:397. 10.1186/1471-2164-12-397
85
MullerT.VingronM. (2000). Modeling amino acid replacement. J. Comput. Biol.7, 761–776. 10.1089/10665270050514918
86
MuseS. V.GautB. S. (1994). A likelihood approach for comparing synonymous and nonsynonymous nucleotide substitution rates, with application to the chloroplast genome. Mol. Biol. Evol.11, 715–724.
87
NickleD. C.HeathL.JensenM. A.GilbertP. B.MullinsJ. I.Kosakovsky PondS. L. (2007). HIV-specific probabilistic models of protein evolution. PLoS ONE2:e503. 10.1371/journal.pone.0000503
88
NielsenR.Bauer DumontV. L.HubiszM. J.AquadroC. F. (2007). Maximum likelihood estimation of ancestral codon usage bias parameters in Drosophila. Mol. Biol. Evol.24, 228–235. 10.1093/molbev/msl146
89
PangA.SmithA. D.NuinP. A.TillierE. R. (2005). SIMPROT: using an empirically determined indel distribution in simulations of protein evolution. BMC Bioinformatics6:236. 10.1186/1471-2105-6-236
90
ParisiG.EchaveJ. (2005). Generality of the structurally constrained protein evolution model: assessment on representatives of the four main fold classes. Gene345, 45–53. 10.1016/j.gene.2004.11.025
91
PellissierL. (2015). Stability and the competition-dispersal trade-off as drivers of speciation and biodiversity gradients. Front. Ecol. Evolution3:52. 10.3389/fevo.2015.00052
92
Perez-JimenezR.Inglés-PrietoA.ZhaoZ. M.Sanchez-RomeroI.Alegre-CebolladaJ.KosuriP.et al. (2011). Single-molecule paleoenzymology probes the chemistry of resurrected enzymes. Nat. Struct. Mol. Biol.18, 592–596. 10.1038/nsmb.2020
93
Perez-LosadaM.JobesD. V.SinangilF.CrandallK. A.ArenasM.PosadaD.et al. (2011). Phylodynamics of HIV-1 from a phase III AIDS vaccine trial in Bangkok, Thailand. PLoS ONE6:e16902. 10.1371/journal.pone.0016902
94
Perez-LosadaM.PosadaD.ArenasM.JobesD. V.SinangilF.BermanP. W.et al. (2009). Ethnic differences in the adaptation rate of HIV gp120 from a vaccine trial. Retrovirology6:67. 10.1186/1742-4690-6-67
95
PondS. K.MuseS. V. (2005). Site-to-site variation of synonymous substitution rates. Mol. Biol. Evol.22, 2375–2385. 10.1093/molbev/msi232
96
PondS. L.FrostS. D. (2005). A genetic algorithm approach to detecting lineage-specific variation in selection pressure. Mol. Biol. Evol.22, 478–485. 10.1093/molbev/msi031
97
PondS. L.FrostS. D.MuseS. V. (2005). HYPHY: hypothesis testing using phylogenies. Bioinformatics21, 676–679. 10.1093/bioinformatics/bti079
98
PosadaD.CrandallK. A. (1998). Modeltest: testing the model of DNA substitution. Bioinformatics14, 817–818. 10.1093/bioinformatics/14.9.817
99
PosadaD.CrandallK. A. (2001). Selecting the best-fit model of nucleotide substitution. Syst. Biol.50, 580–601. 10.1080/106351501750435121
100
RambautA.GrasslyN. C. (1997). Seq-Gen: an application for the Monte Carlo simulation of DNA sequence evolution along phylogenetic trees. Comput. Appl. Biosciences13, 235–238. 10.1093/bioinformatics/13.3.235
101
RastogiS.ReuterN.LiberlesD. A. (2006). Evaluation of models for the evolution of protein sequences and functions under structural constraint. Biophys. Chem.124, 134–144. 10.1016/j.bpc.2006.06.008
102
RodrigueN.PhilippeH.LartillotN. (2010). Mutation-selection models of coding sequence evolution with site-heterogeneous amino acid fitness profiles. Proc. Natl. Acad. Sci. U.S.A.107, 4629–4634. 10.1073/pnas.0910915107
103
RonquistF.TeslenkoM.Van Der MarkP.AyresD. L.DarlingA.HöhnaS.et al. (2012). MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol.61, 539–542. 10.1093/sysbio/sys029
104
RosenbergM. S. (2005). MySSP: non-stationary evolutionary sequence simulation, including indels. Evol. Bioinform. Online1, 81–83. Available online at: http://www.la-press.com/myssp-non-stationary-evolutionary-sequence-simulation-including-indels-article-a185
105
SchneiderA.CannarozziG. M.GonnetG. H. (2005). Empirical codon substitution matrix. BMC Bioinformatics6:134. 10.1186/1471-2105-6-134
106
SeoT. K.KishinoH. (2008). Synonymous substitutions substantially improve evolutionary inference from highly diverged proteins. Syst. Biol.57, 367–377. 10.1080/10635150802158670
107
ShapiroB.RambautA.DrummondA. J. (2006). Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences. Mol. Biol. Evol.23, 7–9. 10.1093/molbev/msj021
108
ShoemakerJ. S.FitchW. M. (1989). Evidence from nuclear sequences that invariable sites should be considered when sequence divergence is calculated. Mol. Biol. Evol.6, 270–289.
109
SiposB.MassinghamT.JordanG. E.GoldmanN. (2011). PhyloSim - Monte Carlo simulation of sequence evolution in the R statistical computing environment. BMC Bioinformatics12:104. 10.1186/1471-2105-12-104
110
StamatakisA. (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics22, 2688–2690. 10.1093/bioinformatics/btl446
111
StropeC. L.AbelK.ScottS. D.MoriyamaE. N. (2009). Biological sequence simulation for testing complex evolutionary hypotheses: indel-Seq-Gen version 2.0. Mol. Biol. Evol.26, 2581–2593. 10.1093/molbev/msp174
112
SumnerJ. G.Fernandez-SanchezJ.JarvisP. D. (2012a). Lie Markov models. J. Theor. Biol.298, 16–31. 10.1016/j.jtbi.2011.12.017
113
SumnerJ. G.JarvisP. D.Fernández-SánchezJ.KaineB. T.WoodhamsM. D.HollandB. R. (2012b). Is the general time-reversible model bad for molecular phylogenetics?Syst. Biol.61, 1069–1074. 10.1093/sysbio/sys042
114
SunnakerM.BusettoA. G.NumminenE.CoranderJ.FollM.DessimozC. (2013). Approximate Bayesian computation. PLoS Comput. Biol.9:e1002803. 10.1371/journal.pcbi.1002803
115
SwoffordD. L. (1993). PAUP: Phylogenetic Analysis Using Parsimony, 3.1.1 Edn. Washington, DC: Smithsonian Institution.
116
TamuraK.StecherG.PetersonD.FilipskiA.KumarS. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. Mol. Biol. Evol.30, 2725–2729. 10.1093/molbev/mst197
117
TavaréS. (1986). Some probabilistic and statistical problems in the analysis of DNA sequences, in Some Mathematical Questions in Biology - DNA Sequence Analysis, ed MiuraR. M. (Providence, RI: Amer Math Soc), 57–86.
118
TavernaD. M.GoldsteinR. A. (2000). The distribution of structures in evolving protein populations. Biopolymers53, 1–8. 10.1002/(SICI)1097-0282(200001)53:1<1::AID-BIP1>3.0.CO;2-X. Available online at: http://onlinelibrary.wiley.com/doi/10.1002/%28SICI%291097-0282%28200001%2953:1%3C1::AID-BIP1%3E3.0.CO;2-X/abstract
119
UsmanovaD. R.FerrettiL.PovolotskayaI. S.VlasovP. K.KondrashovF. A. (2015). A model of substitution trajectories in sequence space and long-term protein evolution. Mol. Biol. Evol.32, 542–554. 10.1093/molbev/msu318
120
WhelanS.GoldmanN. (2001). A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol. Biol. Evol.18, 691–699. 10.1093/oxfordjournals.molbev.a003851
121
WilkeC. O. (2012). Bringing molecules back into molecular evolution. PLoS Comput. Biol.8:e1002572. 10.1371/journal.pcbi.1002572
122
WilsonD. J.McVeanG. (2006). Estimating diversifying selection and functional constraint in the presence of recombination. Genetics172, 1411–1425. 10.1534/genetics.105.044917
123
WongW. S.SainudiinR.NielsenR. (2006). Identification of physicochemical selective pressure on protein encoding nucleotide sequences. BMC Bioinformatics7:148. 10.1186/1471-2105-7-148
124
YangZ.BielawskyJ. P. (2000). Statistical methods for detecting molecular adaptation. Trends Ecol. Evol.15, 496–503. 10.1016/S0169-5347(00)01994-7
125
YangZ.NielsenR. (2008). Mutation-selection models of codon substitution and their use to estimate selective strengths on codon usage. Mol. Biol. Evol.25, 568–579. 10.1093/molbev/msm284
126
YangZ.NielsenR.GoldmanN.PedersenA.-M. K. (2000). Codon-substitution models for heterogeneous selection pressure at amino acid sites. Genetics155, 431–449. Available online at: http://www.genetics.org/content/155/1/431.long
127
YangZ.NielsenR.HasegawaM. (1998). Models of amino acid substitution and applications to mitochondrial protein evolution. Mol. Biol. Evol.15, 1600–1611. 10.1093/oxfordjournals.molbev.a025888
128
YangZ. (1994). Maximum likelihood phylogenetic estimation from DNA sequences with variable rates over sites: approximate methods. J. Mol. Evol.39, 306–314. 10.1007/BF00160154
129
YangZ. (2007). PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol.24, 1586–1591. 10.1093/molbev/msm088
130
ZhaoC.SacanA. (2015). UniAlign: protein structure alignment meets evolution. Bioinformatics31, 3139–3146. 10.1093/bioinformatics/btv354
131
ZharkikhA. (1994). Estimation of evolutionary distances between nucleotide sequences. J. Mol. Evol.39, 315–329. 10.1007/BF00160155
132
ZollerS.BoskovaV.AnisimovaM. (2015). Maximum-likelihood tree estimation using codon substitution models with multiple partitions. Mol. Biol. Evol.32, 2208–2216. 10.1093/molbev/msv097
133
ZollerS.SchneiderA. (2013). Improving phylogenetic inference with a semiempirical amino acid substitution model. Mol. Biol. Evol.30, 469–479. 10.1093/molbev/mss229
Summary
Keywords
substitution model, molecular evolution, model selection, molecular adaptation, phylogenetics, phylogenomics
Citation
Arenas M (2015) Trends in substitution models of molecular evolution. Front. Genet. 6:319. doi: 10.3389/fgene.2015.00319
Received
27 July 2015
Accepted
09 October 2015
Published
26 October 2015
Volume
6 - 2015
Edited by
James J. Cai, Texas A&M University, USA
Reviewed by
Benjamin Callahan, Stanford University, USA; David S. Lawrie, Stanford University, USA; Russell Alan Hermansen, Temple University, USA
Updates
Copyright
© 2015 Arenas.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Miguel Arenas marenas@ipatimup.pt
This article was submitted to Evolutionary and Population Genetics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.