Phylogenetic mixture models for proteins

Autor: Si Quang Le, Nicolas Lartillot, Olivier Gascuel
Přispěvatelé: Laboratoire d'Informatique de Robotique et de Microélectronique de Montpellier (LIRMM), Centre National de la Recherche Scientifique (CNRS)-Université de Montpellier (UM), Méthodes et Algorithmes pour la Bioinformatique (MAB), Centre National de la Recherche Scientifique (CNRS)-Université de Montpellier (UM)-Centre National de la Recherche Scientifique (CNRS)-Université de Montpellier (UM)
Rok vydání: 2008
Předmět:
0106 biological sciences
JTT
Biology
Bioinformatics
010603 evolutionary biology
01 natural sciences
General Biochemistry
Genetics and Molecular Biology

Evolution
Molecular

Matrix (chemical analysis)
03 medical and health sciences
Phylogenetics
phylogenetic inference
Single amino acid
Phylogeny
030304 developmental biology
chemistry.chemical_classification
Likelihood Functions
0303 health sciences
maximum-likelihood estimations
Models
Genetic

Phylogenetic tree
Substitution (logic)
Proteins
Amino acid substitution
CAT profile model
Mixture model
[SDV.BIBS]Life Sciences [q-bio]/Quantitative Methods [q-bio.QM]
Amino acid
amino acid replacement matrices
Amino Acid Substitution
chemistry
Biochemistry
[INFO.INFO-BI]Computer Science [cs]/Bioinformatics [q-bio.QM]
General Agricultural and Biological Sciences
WAG and LG
Research Article
Zdroj: Philosophical Transactions of the Royal Society B: Biological Sciences
Philosophical Transactions of the Royal Society B: Biological Sciences, Royal Society, The, 2008, 363, pp.3965-3976. ⟨10.1098/rstb.2008.0180⟩
ISSN: 1471-2970
0962-8436
DOI: 10.1098/rstb.2008.0180
Popis: Standard protein substitution models use a single amino acid replacement rate matrix that summarizes the biological, chemical and physical properties of amino acids. However, site evolution is highly heterogeneous and depends on many factors: genetic code; solvent exposure; secondary and tertiary structure; protein function; etc. These impact the substitution pattern and, in most cases, a single replacement matrix is not enough to represent all the complexity of the evolutionary processes. This paper explores in maximum-likelihood framework phylogenetic mixture models that combine several amino acid replacement matrices to better fit protein evolution. We learn these mixture models from a large alignment database extracted from HSSP, and test the performance using independent alignments from TreeBase . We compare unsupervised learning approaches, where the site categories are unknown, to supervised ones, where in estimations we use the known category of each site, based on its exposure or its secondary structure. All our models are combined with gamma-distributed rates across sites. Results show that highly significant likelihood gains are obtained when using mixture models compared with the best available single replacement matrices. Mixtures of matrices also improve over mixtures of profiles in the manner of the CAT model. The unsupervised approach tends to be better than the supervised one, but it appears difficult to implement and highly sensitive to the starting values of the parameters, meaning that the supervised approach is still of interest for initialization and model comparison. Using an unsupervised model involving three matrices, the average AIC gain per site with TreeBase test alignments is 0.31, 0.49 and 0.61 compared with LG (named after Le & Gascuel 2008 Mol. Biol. Evol. 25 , 1307–1320), WAG and JTT, respectively. This three-matrix model is significantly better than LG for 34 alignments (among 57), and significantly worse for 1 alignment only. Moreover, tree topologies inferred with our mixture models frequently differ from those obtained with single matrices, indicating that using these mixtures impacts not only the likelihood value but also the output tree. All our models and a PhyML implementation are available from http://atgc.lirmm.fr/mixtures .
Databáze: OpenAIRE