A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia
Identifieur interne : 000612 ( Pmc/Corpus ); précédent : 000611; suivant : 000613A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia
Auteurs : Jonathan J. Fong ; Jeremy M. Brown ; Matthew K. Fujita ; Bastien BoussauSource :
- PLoS ONE [ 1932-6203 ] ; 2012.
Abstract
In resolving the vertebrate tree of life, two fundamental questions remain: 1) what is the phylogenetic position of turtles within amniotes, and 2) what are the relationships between the three major lissamphibian (extant amphibian) groups? These relationships have historically been difficult to resolve, with five different hypotheses proposed for turtle placement, and four proposed branching patterns within Lissamphibia. We compiled a large cDNA/EST dataset for vertebrates (75 genes for 129 taxa) to address these outstanding questions. Gene-specific phylogenetic analyses revealed a great deal of variation in preferred topology, resulting in topologically ambiguous conclusions from the combined dataset. Due to consistent preferences for the same divergent topologies across genes, we suspected systematic phylogenetic error as a cause of some variation. Accordingly, we developed and tested a novel statistical method that identifies sites that have a high probability of containing biased signal for a specific phylogenetic relationship. After removing putatively biased sites, support emerged for a sister relationship between turtles and either crocodilians or archosaurs, as well as for a caecilian-salamander sister relationship within Lissamphibia, with Lissamphibia potentially paraphyletic.
Url:
DOI: 10.1371/journal.pone.0048990
PubMed: 23145043
PubMed Central: 3492174
Links to Exploration step
PMC:3492174Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia</title>
<author><name sortKey="Fong, Jonathan J" sort="Fong, Jonathan J" uniqKey="Fong J" first="Jonathan J." last="Fong">Jonathan J. Fong</name>
<affiliation><nlm:aff id="aff1"><addr-line>Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>College of Natural Sciences, Seoul National University, Seoul, Republic of Korea</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Brown, Jeremy M" sort="Brown, Jeremy M" uniqKey="Brown J" first="Jeremy M." last="Brown">Jeremy M. Brown</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff4"><addr-line>Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Fujita, Matthew K" sort="Fujita, Matthew K" uniqKey="Fujita M" first="Matthew K." last="Fujita">Matthew K. Fujita</name>
<affiliation><nlm:aff id="aff1"><addr-line>Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff5"><addr-line>Museum of Comparative Zoology & Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff6"><addr-line>Department of Biology, University of Texas-Arlington, Arlington, Texas, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Boussau, Bastien" sort="Boussau, Bastien" uniqKey="Boussau B" first="Bastien" last="Boussau">Bastien Boussau</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff7"><addr-line>Laboratorie de Biométrie et Biologie Evolutive, Université de Lyon, Villeurbanne, France</addr-line>
</nlm:aff>
</affiliation>
</author>
</titleStmt>
<publicationStmt><idno type="wicri:source">PMC</idno>
<idno type="pmid">23145043</idno>
<idno type="pmc">3492174</idno>
<idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3492174</idno>
<idno type="RBID">PMC:3492174</idno>
<idno type="doi">10.1371/journal.pone.0048990</idno>
<date when="2012">2012</date>
<idno type="wicri:Area/Pmc/Corpus">000612</idno>
</publicationStmt>
<sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia</title>
<author><name sortKey="Fong, Jonathan J" sort="Fong, Jonathan J" uniqKey="Fong J" first="Jonathan J." last="Fong">Jonathan J. Fong</name>
<affiliation><nlm:aff id="aff1"><addr-line>Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff3"><addr-line>College of Natural Sciences, Seoul National University, Seoul, Republic of Korea</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Brown, Jeremy M" sort="Brown, Jeremy M" uniqKey="Brown J" first="Jeremy M." last="Brown">Jeremy M. Brown</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff4"><addr-line>Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Fujita, Matthew K" sort="Fujita, Matthew K" uniqKey="Fujita M" first="Matthew K." last="Fujita">Matthew K. Fujita</name>
<affiliation><nlm:aff id="aff1"><addr-line>Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff5"><addr-line>Museum of Comparative Zoology & Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff6"><addr-line>Department of Biology, University of Texas-Arlington, Arlington, Texas, United States of America</addr-line>
</nlm:aff>
</affiliation>
</author>
<author><name sortKey="Boussau, Bastien" sort="Boussau, Bastien" uniqKey="Boussau B" first="Bastien" last="Boussau">Bastien Boussau</name>
<affiliation><nlm:aff id="aff2"><addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</nlm:aff>
</affiliation>
<affiliation><nlm:aff id="aff7"><addr-line>Laboratorie de Biométrie et Biologie Evolutive, Université de Lyon, Villeurbanne, France</addr-line>
</nlm:aff>
</affiliation>
</author>
</analytic>
<series><title level="j">PLoS ONE</title>
<idno type="eISSN">1932-6203</idno>
<imprint><date when="2012">2012</date>
</imprint>
</series>
</biblStruct>
</sourceDesc>
</fileDesc>
<profileDesc><textClass></textClass>
</profileDesc>
</teiHeader>
<front><div type="abstract" xml:lang="en"><p>In resolving the vertebrate tree of life, two fundamental questions remain: 1) what is the phylogenetic position of turtles within amniotes, and 2) what are the relationships between the three major lissamphibian (extant amphibian) groups? These relationships have historically been difficult to resolve, with five different hypotheses proposed for turtle placement, and four proposed branching patterns within Lissamphibia. We compiled a large cDNA/EST dataset for vertebrates (75 genes for 129 taxa) to address these outstanding questions. Gene-specific phylogenetic analyses revealed a great deal of variation in preferred topology, resulting in topologically ambiguous conclusions from the combined dataset. Due to consistent preferences for the same divergent topologies across genes, we suspected systematic phylogenetic error as a cause of some variation. Accordingly, we developed and tested a novel statistical method that identifies sites that have a high probability of containing biased signal for a specific phylogenetic relationship. After removing putatively biased sites, support emerged for a sister relationship between turtles and either crocodilians or archosaurs, as well as for a caecilian-salamander sister relationship within Lissamphibia, with Lissamphibia potentially paraphyletic.</p>
</div>
</front>
<back><div1 type="bibliography"><listBibl><biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Benton, Mj" uniqKey="Benton M">MJ Benton</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Meyer, A" uniqKey="Meyer A">A Meyer</name>
</author>
<author><name sortKey="Zardoya, R" uniqKey="Zardoya R">R Zardoya</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Thomson, Rc" uniqKey="Thomson R">RC Thomson</name>
</author>
<author><name sortKey="Shaffer, Hb" uniqKey="Shaffer H">HB Shaffer</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Lyson, Tr" uniqKey="Lyson T">TR Lyson</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Chiari, Y" uniqKey="Chiari Y">Y Chiari</name>
</author>
<author><name sortKey="Cahais, V" uniqKey="Cahais V">V Cahais</name>
</author>
<author><name sortKey="Galtier, N" uniqKey="Galtier N">N Galtier</name>
</author>
<author><name sortKey="Delsuc, F" uniqKey="Delsuc F">F Delsuc</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Parsons, Ts" uniqKey="Parsons T">TS Parsons</name>
</author>
<author><name sortKey="Williams, Ee" uniqKey="Williams E">EE Williams</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Carroll, Rl" uniqKey="Carroll R">RL Carroll</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Vallin, G" uniqKey="Vallin G">G Vallin</name>
</author>
<author><name sortKey="Laurin, M" uniqKey="Laurin M">M Laurin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hedges, Sb" uniqKey="Hedges S">SB Hedges</name>
</author>
<author><name sortKey="Maxson, Lr" uniqKey="Maxson L">LR Maxson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Feller, Ae" uniqKey="Feller A">AE Feller</name>
</author>
<author><name sortKey="Hedges, Sb" uniqKey="Hedges S">SB Hedges</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Milner, Ar" uniqKey="Milner A">AR Milner</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Ruta, M" uniqKey="Ruta M">M Ruta</name>
</author>
<author><name sortKey="Coates, Mi" uniqKey="Coates M">MI Coates</name>
</author>
<author><name sortKey="Quicke, Dlj" uniqKey="Quicke D">DLJ Quicke</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Zardoya, R" uniqKey="Zardoya R">R Zardoya</name>
</author>
<author><name sortKey="Meyer, A" uniqKey="Meyer A">A Meyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hugall, Af" uniqKey="Hugall A">AF Hugall</name>
</author>
<author><name sortKey="Foster, R" uniqKey="Foster R">R Foster</name>
</author>
<author><name sortKey="Lee, Msy" uniqKey="Lee M">MSY Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Roleants, K" uniqKey="Roleants K">K Roleants</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="San Mauro, D" uniqKey="San Mauro D">D San Mauro</name>
</author>
<author><name sortKey="Vences, M" uniqKey="Vences M">M Vences</name>
</author>
<author><name sortKey="Alcobendas, M" uniqKey="Alcobendas M">M Alcobendas</name>
</author>
<author><name sortKey="Zardoya, R" uniqKey="Zardoya R">R Zardoya</name>
</author>
<author><name sortKey="Meyer, A" uniqKey="Meyer A">A Meyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="San Mauro, D" uniqKey="San Mauro D">D San Mauro</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anderson, Js" uniqKey="Anderson J">JS Anderson</name>
</author>
<author><name sortKey="Reisz, Rr" uniqKey="Reisz R">RR Reisz</name>
</author>
<author><name sortKey="Scott, D" uniqKey="Scott D">D Scott</name>
</author>
<author><name sortKey="Frobisch, Nb" uniqKey="Frobisch N">NB Frobisch</name>
</author>
<author><name sortKey="Sumida, Ss" uniqKey="Sumida S">SS Sumida</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Anderson, Js" uniqKey="Anderson J">JS Anderson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Carroll, Rl" uniqKey="Carroll R">RL Carroll</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Carroll, Rl" uniqKey="Carroll R">RL Carroll</name>
</author>
<author><name sortKey="Holmes, R" uniqKey="Holmes R">R Holmes</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Li, C" uniqKey="Li C">C Li</name>
</author>
<author><name sortKey="Wu, X" uniqKey="Wu X">X Wu</name>
</author>
<author><name sortKey="Rieppel, O" uniqKey="Rieppel O">O Rieppel</name>
</author>
<author><name sortKey="Wang, L" uniqKey="Wang L">L Wang</name>
</author>
<author><name sortKey="Zhao, L" uniqKey="Zhao L">L Zhao</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Philippe, H" uniqKey="Philippe H">H Philippe</name>
</author>
<author><name sortKey="Brinkmann, H" uniqKey="Brinkmann H">H Brinkmann</name>
</author>
<author><name sortKey="Lavrov, Dv" uniqKey="Lavrov D">DV Lavrov</name>
</author>
<author><name sortKey="Littlewood, Tj" uniqKey="Littlewood T">TJ Littlewood</name>
</author>
<author><name sortKey="Manuel, M" uniqKey="Manuel M">M Manuel</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Jeffroy, O" uniqKey="Jeffroy O">O Jeffroy</name>
</author>
<author><name sortKey="Brinkmann, H" uniqKey="Brinkmann H">H Brinkmann</name>
</author>
<author><name sortKey="Delsuc, F" uniqKey="Delsuc F">F Delsuc</name>
</author>
<author><name sortKey="Philippe, H" uniqKey="Philippe H">H Philippe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Pattengale, Nd" uniqKey="Pattengale N">ND Pattengale</name>
</author>
<author><name sortKey="Swenson, Km" uniqKey="Swenson K">KM Swenson</name>
</author>
<author><name sortKey="Moret, Bme" uniqKey="Moret B">BME Moret</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Meyer, A" uniqKey="Meyer A">A Meyer</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Simon, C" uniqKey="Simon C">C Simon</name>
</author>
<author><name sortKey="Buckley, Tr" uniqKey="Buckley T">TR Buckley</name>
</author>
<author><name sortKey="Frati, F" uniqKey="Frati F">F Frati</name>
</author>
<author><name sortKey="Stewart, Jb" uniqKey="Stewart J">JB Stewart</name>
</author>
<author><name sortKey="Beckenbach, At" uniqKey="Beckenbach A">AT Beckenbach</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Edwards, Sv" uniqKey="Edwards S">SV Edwards</name>
</author>
<author><name sortKey="Arctander, P" uniqKey="Arctander P">P Arctander</name>
</author>
<author><name sortKey="Wilson, Ac" uniqKey="Wilson A">AC Wilson</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Blouin, Ms" uniqKey="Blouin M">MS Blouin</name>
</author>
<author><name sortKey="Yowell, Ca" uniqKey="Yowell C">CA Yowell</name>
</author>
<author><name sortKey="Courtnely, Ch" uniqKey="Courtnely C">CH Courtnely</name>
</author>
<author><name sortKey="Dame, Jb" uniqKey="Dame J">JB Dame</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Regier, Jc" uniqKey="Regier J">JC Regier</name>
</author>
<author><name sortKey="Schultz, Jw" uniqKey="Schultz J">JW Schultz</name>
</author>
<author><name sortKey="Zwick, A" uniqKey="Zwick A">A Zwick</name>
</author>
<author><name sortKey="Hussey, A" uniqKey="Hussey A">A Hussey</name>
</author>
<author><name sortKey="Ball, B" uniqKey="Ball B">B Ball</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Brinkmann, H" uniqKey="Brinkmann H">H Brinkmann</name>
</author>
<author><name sortKey="Philippe, H" uniqKey="Philippe H">H Philippe</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Philippe, H" uniqKey="Philippe H">H Philippe</name>
</author>
<author><name sortKey="Lartillot, N" uniqKey="Lartillot N">N Lartillot</name>
</author>
<author><name sortKey="Brinkmann, H" uniqKey="Brinkmann H">H Brinkmann</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Fong, Jj" uniqKey="Fong J">JJ Fong</name>
</author>
<author><name sortKey="Fujita, Mk" uniqKey="Fujita M">MK Fujita</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Edwards, Sv" uniqKey="Edwards S">SV Edwards</name>
</author>
<author><name sortKey="Liu, L" uniqKey="Liu L">L Liu</name>
</author>
<author><name sortKey="Pearl, Dk" uniqKey="Pearl D">DK Pearl</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Cranston, Ka" uniqKey="Cranston K">KA Cranston</name>
</author>
<author><name sortKey="Hurwitz, B" uniqKey="Hurwitz B">B Hurwitz</name>
</author>
<author><name sortKey="Ware, D" uniqKey="Ware D">D Ware</name>
</author>
<author><name sortKey="Stein, L" uniqKey="Stein L">L Stein</name>
</author>
<author><name sortKey="Wing, Ra" uniqKey="Wing R">RA Wing</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Shimodaira, H" uniqKey="Shimodaira H">H Shimodaira</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Rodriguez Ezpeleta, N" uniqKey="Rodriguez Ezpeleta N">N Rodriguez-Ezpeleta</name>
</author>
<author><name sortKey="Brinkmann, H" uniqKey="Brinkmann H">H Brinkmann</name>
</author>
<author><name sortKey="Roure, B" uniqKey="Roure B">B Roure</name>
</author>
<author><name sortKey="Lartillot, N" uniqKey="Lartillot N">N Lartillot</name>
</author>
<author><name sortKey="Lang, Bf" uniqKey="Lang B">BF Lang</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lemmon, Ar" uniqKey="Lemmon A">AR Lemmon</name>
</author>
<author><name sortKey="Brown, Jm" uniqKey="Brown J">JM Brown</name>
</author>
<author><name sortKey="Stanger Hall, K" uniqKey="Stanger Hall K">K Stanger-Hall</name>
</author>
<author><name sortKey="Lemmon, Em" uniqKey="Lemmon E">EM Lemmon</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Wolf, Yi" uniqKey="Wolf Y">YI Wolf</name>
</author>
<author><name sortKey="Rogozin, Ib" uniqKey="Rogozin I">IB Rogozin</name>
</author>
<author><name sortKey="Koonin, Ev" uniqKey="Koonin E">EV Koonin</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Aguinaldo, Am" uniqKey="Aguinaldo A">AM Aguinaldo</name>
</author>
<author><name sortKey="Tubeville, Jm" uniqKey="Tubeville J">JM Tubeville</name>
</author>
<author><name sortKey="Linford, Ls" uniqKey="Linford L">LS Linford</name>
</author>
<author><name sortKey="Rivera, Mc" uniqKey="Rivera M">MC Rivera</name>
</author>
<author><name sortKey="Raff, Ra" uniqKey="Raff R">RA Raff</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Telford, Mj" uniqKey="Telford M">MJ Telford</name>
</author>
<author><name sortKey="Bourlat, Sj" uniqKey="Bourlat S">SJ Bourlat</name>
</author>
<author><name sortKey="Economou, A" uniqKey="Economou A">A Economou</name>
</author>
<author><name sortKey="Papillon, D" uniqKey="Papillon D">D Papillon</name>
</author>
<author><name sortKey="Rota Stabelli, O" uniqKey="Rota Stabelli O">O Rota-Stabelli</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Bhullar, B A" uniqKey="Bhullar B">B-A Bhullar</name>
</author>
<author><name sortKey="Bever, Gs" uniqKey="Bever G">GS Bever</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Laurin, M" uniqKey="Laurin M">M Laurin</name>
</author>
<author><name sortKey="Reisz, R" uniqKey="Reisz R">R Reisz</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lee, Msy" uniqKey="Lee M">MSY Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Lee, Msy" uniqKey="Lee M">MSY Lee</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Nagashima, H" uniqKey="Nagashima H">H Nagashima</name>
</author>
<author><name sortKey="Sugahara, F" uniqKey="Sugahara F">F Sugahara</name>
</author>
<author><name sortKey="Takechi, M" uniqKey="Takechi M">M Takechi</name>
</author>
<author><name sortKey="Ericsson, R" uniqKey="Ericsson R">R Ericsson</name>
</author>
<author><name sortKey="Kawashima Ohya, Y" uniqKey="Kawashima Ohya Y">Y Kawashima-Ohya</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Yates, Am" uniqKey="Yates A">AM Yates</name>
</author>
<author><name sortKey="Warren, A" uniqKey="Warren A">A Warren</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hedges, Sb" uniqKey="Hedges S">SB Hedges</name>
</author>
<author><name sortKey="Poling, Ll" uniqKey="Poling L">LL Poling</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Hubbard, Tjp" uniqKey="Hubbard T">TJP Hubbard</name>
</author>
<author><name sortKey="Aken, Bl" uniqKey="Aken B">BL Aken</name>
</author>
<author><name sortKey="Ayling, S" uniqKey="Ayling S">S Ayling</name>
</author>
<author><name sortKey="Ballester, B" uniqKey="Ballester B">B Ballester</name>
</author>
<author><name sortKey="Beal, K" uniqKey="Beal K">K Beal</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stamatakis, A" uniqKey="Stamatakis A">A Stamatakis</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Edgar, Rc" uniqKey="Edgar R">RC Edgar</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Ronquist, F" uniqKey="Ronquist F">F Ronquist</name>
</author>
<author><name sortKey="Huelsenbeck, Jp" uniqKey="Huelsenbeck J">JP Huelsenbeck</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Abascal, F" uniqKey="Abascal F">F Abascal</name>
</author>
<author><name sortKey="Zardoya, R" uniqKey="Zardoya R">R Zardoya</name>
</author>
<author><name sortKey="Posada, D" uniqKey="Posada D">D Posada</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Stamatakis, A" uniqKey="Stamatakis A">A Stamatakis</name>
</author>
<author><name sortKey="Hoover, P" uniqKey="Hoover P">P Hoover</name>
</author>
<author><name sortKey="Rougemont, J" uniqKey="Rougemont J">J Rougemont</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Shimodaira, H" uniqKey="Shimodaira H">H Shimodaira</name>
</author>
<author><name sortKey="Hasegawa, M" uniqKey="Hasegawa M">M Hasegawa</name>
</author>
</analytic>
</biblStruct>
<biblStruct></biblStruct>
<biblStruct></biblStruct>
<biblStruct><analytic><author><name sortKey="Pond, Slk" uniqKey="Pond S">SLK Pond</name>
</author>
<author><name sortKey="Frost, Sdw" uniqKey="Frost S">SDW Frost</name>
</author>
<author><name sortKey="Muse, Sv" uniqKey="Muse S">SV Muse</name>
</author>
</analytic>
</biblStruct>
<biblStruct><analytic><author><name sortKey="Dutheil, J" uniqKey="Dutheil J">J Dutheil</name>
</author>
<author><name sortKey="Boussau, B" uniqKey="Boussau B">B Boussau</name>
</author>
</analytic>
</biblStruct>
</listBibl>
</div1>
</back>
</TEI>
<pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">PLoS One</journal-id>
<journal-id journal-id-type="iso-abbrev">PLoS ONE</journal-id>
<journal-id journal-id-type="publisher-id">plos</journal-id>
<journal-id journal-id-type="pmc">plosone</journal-id>
<journal-title-group><journal-title>PLoS ONE</journal-title>
</journal-title-group>
<issn pub-type="epub">1932-6203</issn>
<publisher><publisher-name>Public Library of Science</publisher-name>
<publisher-loc>San Francisco, USA</publisher-loc>
</publisher>
</journal-meta>
<article-meta><article-id pub-id-type="pmid">23145043</article-id>
<article-id pub-id-type="pmc">3492174</article-id>
<article-id pub-id-type="publisher-id">PONE-D-12-10241</article-id>
<article-id pub-id-type="doi">10.1371/journal.pone.0048990</article-id>
<article-categories><subj-group subj-group-type="heading"><subject>Research Article</subject>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Biology</subject>
<subj-group><subject>Computational Biology</subject>
<subj-group><subject>Sequence Analysis</subject>
</subj-group>
</subj-group>
<subj-group><subject>Evolutionary Biology</subject>
<subj-group><subject>Evolutionary Systematics</subject>
<subj-group><subject>Molecular Systematics</subject>
<subject>Phylogenetics</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group><subject>Genomics</subject>
<subj-group><subject>Genome Databases</subject>
<subj-group><subject>Sequence Databases</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group><subject>Zoology</subject>
<subj-group><subject>Herpetology</subject>
</subj-group>
</subj-group>
</subj-group>
<subj-group subj-group-type="Discipline-v2"><subject>Mathematics</subject>
<subj-group><subject>Statistics</subject>
<subj-group><subject>Statistical Methods</subject>
</subj-group>
</subj-group>
</subj-group>
</article-categories>
<title-group><article-title>A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia</article-title>
<alt-title alt-title-type="running-head">Phylogenomics of the Vertebrate Phylogeny</alt-title>
</title-group>
<contrib-group><contrib contrib-type="author"><name><surname>Fong</surname>
<given-names>Jonathan J.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3"><sup>3</sup>
</xref>
<xref ref-type="corresp" rid="cor1"><sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Brown</surname>
<given-names>Jeremy M.</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff4"><sup>4</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Fujita</surname>
<given-names>Matthew K.</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff5"><sup>5</sup>
</xref>
<xref ref-type="aff" rid="aff6"><sup>6</sup>
</xref>
</contrib>
<contrib contrib-type="author"><name><surname>Boussau</surname>
<given-names>Bastien</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff7"><sup>7</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1"><label>1</label>
<addr-line>Museum of Vertebrate Zoology, University of California, Berkeley, California, United States of America</addr-line>
</aff>
<aff id="aff2"><label>2</label>
<addr-line>Department of Integrative Biology, University of California, Berkeley, California, United States of America</addr-line>
</aff>
<aff id="aff3"><label>3</label>
<addr-line>College of Natural Sciences, Seoul National University, Seoul, Republic of Korea</addr-line>
</aff>
<aff id="aff4"><label>4</label>
<addr-line>Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, United States of America</addr-line>
</aff>
<aff id="aff5"><label>5</label>
<addr-line>Museum of Comparative Zoology & Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America</addr-line>
</aff>
<aff id="aff6"><label>6</label>
<addr-line>Department of Biology, University of Texas-Arlington, Arlington, Texas, United States of America</addr-line>
</aff>
<aff id="aff7"><label>7</label>
<addr-line>Laboratorie de Biométrie et Biologie Evolutive, Université de Lyon, Villeurbanne, France</addr-line>
</aff>
<contrib-group><contrib contrib-type="editor"><name><surname>Hejnol</surname>
<given-names>Andreas</given-names>
</name>
<role>Editor</role>
<xref ref-type="aff" rid="edit1"></xref>
</contrib>
</contrib-group>
<aff id="edit1"><addr-line>Sars International Centre for Marine Molecular Biology, Norway</addr-line>
</aff>
<author-notes><corresp id="cor1">* E-mail: <email>jfong@snu.ac.kr</email>
</corresp>
<fn fn-type="conflict"><p><bold>Competing Interests: </bold>
The authors have declared that no competing interests exist.</p>
</fn>
<fn fn-type="con"><p>Conceived and designed the experiments: JJF MKF. Performed the experiments: JJF MKF. Analyzed the data: JJF MKF JMB BB. Contributed reagents/materials/analysis tools: JMB BB. Wrote the paper: JJF MKF JMB BB.</p>
</fn>
</author-notes>
<pub-date pub-type="collection"><year>2012</year>
</pub-date>
<pub-date pub-type="epub"><day>7</day>
<month>11</month>
<year>2012</year>
</pub-date>
<volume>7</volume>
<issue>11</issue>
<elocation-id>e48990</elocation-id>
<history><date date-type="received"><day>10</day>
<month>4</month>
<year>2012</year>
</date>
<date date-type="accepted"><day>3</day>
<month>10</month>
<year>2012</year>
</date>
</history>
<permissions><copyright-year>2012</copyright-year>
<copyright-holder>Fong et al</copyright-holder>
<license><license-p>This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.</license-p>
</license>
</permissions>
<abstract><p>In resolving the vertebrate tree of life, two fundamental questions remain: 1) what is the phylogenetic position of turtles within amniotes, and 2) what are the relationships between the three major lissamphibian (extant amphibian) groups? These relationships have historically been difficult to resolve, with five different hypotheses proposed for turtle placement, and four proposed branching patterns within Lissamphibia. We compiled a large cDNA/EST dataset for vertebrates (75 genes for 129 taxa) to address these outstanding questions. Gene-specific phylogenetic analyses revealed a great deal of variation in preferred topology, resulting in topologically ambiguous conclusions from the combined dataset. Due to consistent preferences for the same divergent topologies across genes, we suspected systematic phylogenetic error as a cause of some variation. Accordingly, we developed and tested a novel statistical method that identifies sites that have a high probability of containing biased signal for a specific phylogenetic relationship. After removing putatively biased sites, support emerged for a sister relationship between turtles and either crocodilians or archosaurs, as well as for a caecilian-salamander sister relationship within Lissamphibia, with Lissamphibia potentially paraphyletic.</p>
</abstract>
<funding-group><funding-statement>Funding for this work came from the National Science Foundation (NSF) Doctoral Dissertation Improvement Grant (DEB-0909811 [JJF]), NSF Postdoctoral Fellowship in Biology (DBI-0905867 [JMB], DBI-0905714 [MKF]), Human Frontiers Science Program Postdoctoral Fellowship (BB), Centre National de la Recherche Scientifique (BB), and Museum of Vertebrate Zoology (JJF). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</funding-statement>
</funding-group>
<counts><page-count count="14"></page-count>
</counts>
</article-meta>
</front>
<body><sec id="s1"><title>Introduction</title>
<p>“The Origin of Species,” and in particular its singular figure, transformed our thinking of biological diversity from the “great chain of being” to the “tree of life” <xref ref-type="bibr" rid="pone.0048990-Darwin1">[1]</xref>
. Resolving the tree of life is crucial to understand organismal evolution and adaptation, but also has far-reaching benefits to diverse fields such as medicine, conservation, and economics <xref ref-type="bibr" rid="pone.0048990-Yates1">[2]</xref>
. While vertebrates have been the focus of intense phylogenetic research <xref ref-type="bibr" rid="pone.0048990-Benton1">[3]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Thomson1">[5]</xref>
, two fundamental questions in vertebrate systematics remain unanswered: 1) What is the phylogenetic position of turtles within amniotes, and 2) what are the relationships between the three major lissamphibian (extant amphibian) groups–frogs, salamanders, and caecilians?</p>
<p>For more than 150 years, biologists have debated the phylogenetic position of turtles, resulting in no fewer than five different hypotheses (<xref ref-type="fig" rid="pone-0048990-g001">Figure 1A</xref>
) <xref ref-type="bibr" rid="pone.0048990-Meyer1">[4]</xref>
. Earlier studies used the number of temporal skull openings for classification, with the anapsid condition (no openings) found in turtles, the synapsid condition (single opening) found in mammals, and the diapsid condition (two openings) found in birds and non-turtle reptiles <xref ref-type="bibr" rid="pone.0048990-Gaffney1">[6]</xref>
. Morphological and molecular data have suggested four additional hypotheses: turtles as basal sauropsids (reptiles and birds), a turtle-lepidosaur (lizards, snakes, amphisbaenians, and tuatara) sister relationship, a turtle-archosaur (birds and crocodilians) sister relationship, and a turtle-crocodilian sister relationship (<xref ref-type="fig" rid="pone-0048990-g001">Figure 1A</xref>
) (see <xref ref-type="bibr" rid="pone.0048990-Meyer1">[4]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Lyson1">[7]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Chiari1">[9]</xref>
for summary of references). Although recent studies have found strong results supporting specific hypotheses, there is no consensus as different datasets support different hypotheses <xref ref-type="bibr" rid="pone.0048990-Lyson1">[7]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Chiari1">[9]</xref>
.</p>
<fig id="pone-0048990-g001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0048990.g001</object-id>
<label>Figure 1</label>
<caption><title>Alternative hypotheses in the vertebrate phylogeny.</title>
<p>Uncertainties in the vertebrate phylogeny examined in this study. (A) The five alternative hypotheses for the placement of turtles within amniotes 1) turtles as basal amniotes, 2) turtles as basal sauropsids, 3) turtle-lepidosaur sister group, 4) turtle-archosaur sister group, and 5) turtle-crocodilian sister group. (B) monophyletic and (C) paraphyletic alternative hypotheses for lissamphibian (extant amphibians) relationships.</p>
</caption>
<graphic xlink:href="pone.0048990.g001"></graphic>
</fig>
<p>For amphibians, several morphological and physiological characters, including pedicellate teeth and cutaneous respiration, suggest frogs, salamanders, and caecilians share a common origin <xref ref-type="bibr" rid="pone.0048990-Parsons1">[10]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Carroll1">[11]</xref>
. However, the monophyly of Lissamphibia is still under debate, as some paleontological studies have inferred a paraphyletic Lissamphibia <xref ref-type="bibr" rid="pone.0048990-Carroll2">[12]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Cannatella1">[13]</xref>
. There are four proposed branching patterns within Lissamphibia (<xref ref-type="fig" rid="pone-0048990-g001">Figure 1B,C</xref>
). Two hypotheses, Procera and Batrachia, exhibit a monophyletic Lissamphibia, but differ in the interrelationships among frogs, salamanders, and caecilians. The Procera hypothesis proposes a salamander-caecilian sister relationship (morphology: <xref ref-type="bibr" rid="pone.0048990-Vallin1">[14]</xref>
; mitochondrial DNA: <xref ref-type="bibr" rid="pone.0048990-Hedges1">[15]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Feller1">[16]</xref>
), while the Batrachia hypothesis proposes a frog-salamander sister relationship (morphology: <xref ref-type="bibr" rid="pone.0048990-Milner1">[17]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Ruta1">[19]</xref>
; nuclear and combined DNA: <xref ref-type="bibr" rid="pone.0048990-Zardoya1">[20]</xref>
–<xref ref-type="bibr" rid="pone.0048990-SanMauro2">[24]</xref>
) (<xref ref-type="fig" rid="pone-0048990-g001">Figure 1B</xref>
). Conversely, two hypotheses based primarily on paleontological data suggest that Lissamphibia is paraphyletic because of an affinity between caecilians and amniotes (<xref ref-type="fig" rid="pone-0048990-g001">Figure 1C</xref>
) <xref ref-type="bibr" rid="pone.0048990-Carroll2">[12]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Cannatella1">[13]</xref>
, with salamanders sister to either frogs <xref ref-type="bibr" rid="pone.0048990-Anderson1">[25]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Carroll3">[27]</xref>
or caecilians <xref ref-type="bibr" rid="pone.0048990-Carroll4">[28]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Carroll5">[29]</xref>
. In general, paleontological data support a paraphyletic Lissamphibia, while molecular data support the Batrachia hypothesis.</p>
<p>Both turtles and lissamphibians have ancient divergences within vertebrates (>200 Ma for turtles, frogs, salamanders, and caecilians) <xref ref-type="bibr" rid="pone.0048990-Carroll1">[11]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Li1">[30]</xref>
and highly modified morphologies. The lack of intermediate forms, either fossil or extant, obscures any obvious morphological evidence of their respective ancestries. Therefore, molecular studies are the best option for uncovering the information necessary to resolve the enigmatic phylogenetic positions of these groups. However, molecular data are not perfect and exhibit several potential pitfalls, especially when trying to resolve difficult phylogenetic questions <xref ref-type="bibr" rid="pone.0048990-Philippe1">[31]</xref>
. For instance, stochastic error (from insufficient data) and/or systematic error (from inadequate models of substitution), can lead to erroneous inferences <xref ref-type="bibr" rid="pone.0048990-Jeffroy1">[32]</xref>
. Rogue taxa (i.e., taxa with strong support for multiple phylogenetic positions due to either variation across genes or systematic error) can also impede phylogenetic inference by appearing to reduce confidence in other relationships <xref ref-type="bibr" rid="pone.0048990-Pattengale1">[33]</xref>
.</p>
<p>In this study, we address two difficult phylogenetic questions in the vertebrate phylogeny: the placement of turtles among amniotes and the relationships among frogs, salamanders, and caecilians. Minimizing stochastic error requires acquiring a sizeable dataset suitable for testing the hypotheses of interest. We do this in our study by compiling one of the largest datasets for vertebrate phylogenetics to date (75 genes for 129 taxa). Systematic error is more difficult to address and not solvable by acquiring additional data <xref ref-type="bibr" rid="pone.0048990-Philippe1">[31]</xref>
. Some generalized approaches to address systematic error involve transforming data <xref ref-type="bibr" rid="pone.0048990-Meyer2">[34]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Regier1">[38]</xref>
or removing fast evolving sites <xref ref-type="bibr" rid="pone.0048990-Brinkmann1">[39]</xref>
or genes <xref ref-type="bibr" rid="pone.0048990-Philippe2">[40]</xref>
to reduce homoplasy. We test such methods with limited levels of success. We therefore develop a new method to remove those data that are most likely to harbor non-phylogenetic signal. Instead of a generalized removal of fast sites or genes, we take a different approach that identifies and removes sites that have a high probability of containing biased signal for a specific phylogenetic relationship. By minimizing non-phylogenetic signal and removing rogue taxa, ambiguity regarding the preferred hypotheses was greatly reduced, allowing us to infer a sister relationship between turtles and either crocodilians or archosaurs, as well as a sister relationship between caecilians and salamanders. Some support was also found for a paraphyletic Lissamphibia, with the caecilian-salamander clade more closely related to amniotes than to frogs.</p>
</sec>
<sec id="s2"><title>Results</title>
<sec id="s2a"><title>Dataset Characteristics</title>
<p>We obtained DNA sequences of 75 protein-coding genes for 129 taxa from 1) online genomic resources and 2) targeted sequencing of new samples <xref ref-type="bibr" rid="pone.0048990-Fong1">[41]</xref>
. Taxon sampling spans all vertebrates, but is skewed towards mammals (36, available data online) and turtles (45, for future studies within turtles). The concatenated alignment of all genes includes 33,938 base positions, and the overall matrix completeness is 41.6% for a total of 4,378,002 bp of sequence data. On average, taxa in the dataset include 31 of the 75 genes. Of the 3,989 sequences, 878 are new (Genbank #s: JF264630-264720, JN864096-864759, JN885182-885183), while the remaining 3,111 are from online resources. In addition to a standard nucleotide dataset [NUCL], we used three transformations of our data to minimize homoplasy for deeper evolutionary divergences: 1) amino acids (AA) <xref ref-type="bibr" rid="pone.0048990-Meyer2">[34]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Simon1">[35]</xref>
, 2) first and second codon positions (N12) <xref ref-type="bibr" rid="pone.0048990-Edwards1">[36]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Blouin1">[37]</xref>
, and 3) sequence adjustment to account for codon degeneracy (DEGEN1) <xref ref-type="bibr" rid="pone.0048990-Regier1">[38]</xref>
. In addition, for each of these data transformations, we applied four alternative taxonomic and gene sampling strategies: 1) all taxa for all 75 genes (All taxa-75 genes), 2) a subset of 16 taxa (see below) for all 75 genes (16 taxa-75genes), 3) all taxa for a reduced set of genes with sequences from all major taxonomic groups for the turtle (All taxa-31 genes) and lissamphibian (All taxa-26 genes) questions, and 4) 16 taxa for each reduced gene set (Turtle: 16 taxa-31 genes, Lissamphibia: 16 taxa-26 genes). For the 16-taxon datasets, only 2–4 of the most data-complete taxa were included from each major taxonomic group to explore the impacts of missing data, as this dataset is 81.9% complete. Lastly, to investigate the effect of fast evolving genes on phylogenetic reconstruction, we calculated the rate of evolution of each gene and removed the fastest genes from our analyses. Based on the shape of the frequency histogram (<xref ref-type="supplementary-material" rid="pone.0048990.s002">Fig. S2</xref>
), we removed the fastest 25% of genes and concatenated the remaining genes (dataset named ‘slow genes’) for analysis. All data files have been deposited in the Dryad Repository: <ext-link ext-link-type="uri" xlink:href="http://dx.doi.org/10.5061/dryad.25j6h">http://dx.doi.org/10.5061/dryad.25j6h</ext-link>
.</p>
</sec>
<sec id="s2b"><title>Phylogenetic Analyses and Topology Tests</title>
<p>We inferred gene trees for each gene using maximum likelihood (RAxML); each hypothesis for both the turtle placement and lissamphibian relationships was supported by subsets of these individual gene trees (<xref ref-type="fig" rid="pone-0048990-g002">Figure 2</xref>
). For the phylogenetic placement of turtles, when all major amniote groups were included, the turtle-crocodilian topology was most frequently supported (29%), while the turtle-bird topology was most common (49%) when no crocodilian was available (<xref ref-type="fig" rid="pone-0048990-g002">Figure 2A,B</xref>
). For lissamphibian relationships, the Batrachia hypothesis was most frequently supported by individual gene analyses (18%), but support was more evenly spread across each of the hypotheses compared to the turtle question (<xref ref-type="fig" rid="pone-0048990-g002">Figure 2C</xref>
). Coalescent-based analyses such as BEST <xref ref-type="bibr" rid="pone.0048990-Edwards2">[42]</xref>
can be used to deal with incongruent gene trees, but were not possible for our datasets, as analyses failed to converge (potentially due to the large amount of missing data <xref ref-type="bibr" rid="pone.0048990-Cranston1">[43]</xref>
). So we focus on the concatenated analyses and their results for the remainder of this paper. Topology tests (approximately unbiased topology tests [AU tests]) <xref ref-type="bibr" rid="pone.0048990-Shimodaira1">[44]</xref>
were run alongside phylogenetic analyses to rigorously test whether the maximum likelihood tree is statistically better than alternative topologies.</p>
<fig id="pone-0048990-g002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0048990.g002</object-id>
<label>Figure 2</label>
<caption><title>Phylogenetic results from individual gene analyses.</title>
<p>(A) The phylogenetic position of turtles within amniotes when all major groups were present and (B) when no crocodilians were present. (C) The relationships between major lissamphibian groups. The “other” category includes topologies that do not match any of the previously proposed hypotheses, usually with a major amniote group being paraphyletic.</p>
</caption>
<graphic xlink:href="pone.0048990.g002"></graphic>
</fig>
<p>For the phylogenetic placement of turtles, results from concatenated datasets were generally consistent within, but differed between, each data transformation (<xref ref-type="supplementary-material" rid="pone.0048990.s001">Figure S1</xref>
, <xref ref-type="table" rid="pone-0048990-t001">Table 1</xref>
). The NUCL data-type recovered a turtle-crocodilian relationship, which was statistically indistinguishable from a turtle-archosaur sister relationship based on AU tests in all datasets (<xref ref-type="table" rid="pone-0048990-t001">Table 1</xref>
). For the N12 data-type, a turtle-archosaur relationship was recovered in all but one dataset (16 taxa-31 genes: turtle-lepidosaur), and AU tests statistically excluded turtles as basal amniotes. For the DEGEN1 and AA data-types, the turtle-lepidosaur relationship was recovered for all datasets, but results from AU tests were not consistent.</p>
<table-wrap id="pone-0048990-t001" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0048990.t001</object-id>
<label>Table 1</label>
<caption><title>Summary of phylogenetic results from different datasets.</title>
</caption>
<alternatives><graphic id="pone-0048990-t001-1" xlink:href="pone.0048990.t001"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td colspan="4" align="left" rowspan="1">NUCL</td>
<td colspan="4" align="left" rowspan="1">N12</td>
<td colspan="4" align="left" rowspan="1">DEGEN1</td>
<td colspan="4" align="left" rowspan="1">AA</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">TURTLE</td>
<td align="left" rowspan="1" colspan="1">Individual</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">All taxa</td>
<td align="left" rowspan="1" colspan="1">16 taxa</td>
<td align="left" rowspan="1" colspan="1">Slow</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1">Genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">75 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">31 genes</td>
<td align="left" rowspan="1" colspan="1">Genes</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">Basal Amniote</td>
<td align="left" rowspan="1" colspan="1">0%</td>
<td align="left" rowspan="1" colspan="1">0.003</td>
<td align="left" rowspan="1" colspan="1">2e-40</td>
<td align="left" rowspan="1" colspan="1">3e-7</td>
<td align="left" rowspan="1" colspan="1">1e-59</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">1e-5</td>
<td align="left" rowspan="1" colspan="1">8e-5</td>
<td align="left" rowspan="1" colspan="1">0.005</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">0.005</td>
<td align="left" rowspan="1" colspan="1"><bold>0.274</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.237</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">6e-6</td>
<td align="left" rowspan="1" colspan="1"><bold>0.477</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.32</bold>
</td>
<td align="left" rowspan="1" colspan="1">5e-6</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Basal Sauropsid</td>
<td align="left" rowspan="1" colspan="1">13%</td>
<td align="left" rowspan="1" colspan="1">3e-4</td>
<td align="left" rowspan="1" colspan="1">1e-4</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1"><bold>0.056</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.009</td>
<td align="left" rowspan="1" colspan="1"><bold>0.188</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.32</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.383</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.07</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.351</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.424</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.454</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.117</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.137</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.398</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.586</bold>
</td>
<td align="left" rowspan="1" colspan="1">5e-4</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Turtle-Lepidosaur</td>
<td align="left" rowspan="1" colspan="1">10%</td>
<td align="left" rowspan="1" colspan="1">2e-5</td>
<td align="left" rowspan="1" colspan="1">1e-4</td>
<td align="left" rowspan="1" colspan="1">2e-4</td>
<td align="left" rowspan="1" colspan="1">0.016</td>
<td align="left" rowspan="1" colspan="1"><bold>0.144</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.339</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.412</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1">1e-6</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Turtle-Archosaur</td>
<td align="left" rowspan="1" colspan="1">19%</td>
<td align="left" rowspan="1" colspan="1"><bold>0.143</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.177</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.327</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.471</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.505</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.601</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.151</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.526</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.041</td>
<td align="left" rowspan="1" colspan="1"><bold>0.345</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.295</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.283</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.027</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.477</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Turtle-Crocodilian</td>
<td align="left" rowspan="1" colspan="1">29%</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.428</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.489</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.448</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.31</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.094</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.062</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.152</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.007</td>
<td align="left" rowspan="1" colspan="1"><bold>0.094</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.038</td>
<td align="left" rowspan="1" colspan="1"><bold>0.167</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.006</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Other</td>
<td align="left" rowspan="1" colspan="1">29%</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td colspan="4" align="left" rowspan="1"><bold>NUCL</bold>
</td>
<td colspan="4" align="left" rowspan="1"><bold>N12</bold>
</td>
<td colspan="4" align="left" rowspan="1"><bold>DEGEN1</bold>
</td>
<td colspan="4" align="left" rowspan="1"><bold>AA</bold>
</td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold>LISSAMPHIBIA</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>Individual</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>All taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>16 taxa</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>Slow</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"><bold>Genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>75 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>26 genes</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>Genes</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Batrachia</td>
<td align="left" rowspan="1" colspan="1">18%</td>
<td align="left" rowspan="1" colspan="1">0.009</td>
<td align="left" rowspan="1" colspan="1">0.008</td>
<td align="left" rowspan="1" colspan="1"><bold>0.204</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.402</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.643</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.539</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.644</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.234</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Procera</td>
<td align="left" rowspan="1" colspan="1">12%</td>
<td align="left" rowspan="1" colspan="1">0.002</td>
<td align="left" rowspan="1" colspan="1">0.003</td>
<td align="left" rowspan="1" colspan="1">0.037</td>
<td align="left" rowspan="1" colspan="1"><bold>0.069</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.311</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.339</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.29</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.183</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.394</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.419</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.215</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.334</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.343</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.169</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.094</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.087</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Paraphyletic</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Frog-Salamander</td>
<td align="left" rowspan="1" colspan="1">9%</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">0.004</td>
<td align="left" rowspan="1" colspan="1"><bold>0.153</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.512</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.587</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.273</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.403</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.386</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.087</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.269</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.074</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.337</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.361</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.235</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.188</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Paraphyletic</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Caecilian-Salamander</td>
<td align="left" rowspan="1" colspan="1">12%</td>
<td align="left" rowspan="1" colspan="1"><bold>0.103</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.209</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.184</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.093</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.185</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.198</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.061</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.013</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.166</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.093</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.074</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.012</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.23</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.224</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.042</td>
<td align="left" rowspan="1" colspan="1"><bold>0.018</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.873</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Other</td>
<td align="left" rowspan="1" colspan="1">48%</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt101"><label></label>
<p>Summary of results from phylogenetic analyses for turtle placement within the amniote phylogeny, and relationships of Lissamphibia (extant amphibians). Individual gene and concatenated analyses were performed. For individual gene analyses, percentages denote the proportion of genes supporting the hypothesis. For the concatenated analyses, the cell with “X” for each column denote the most likely topology based on RAxML likelihood scores, while the numbers in cells represent p-values based on approximately unbiased topology tests (Shimodaira 2002). Cells in bold font denote statistically indistinguishable topologies (p-value>0.05) from the most likely topology. Sixteen total concatenated analyses were performed (4 data transformations × 4 datasets). Data were transformed in an attempt to reduce the rate of evolution: NUCL = complete nucleotide dataset, N12 = 1<sup>st</sup>
and 2<sup>nd</sup>
codon positions only (Edwards et al. 1991, Blouin et al. 1998), DEGEN1 = codon degeneracy (Regier et al. 2010), AA = amino acids (Meyer 1994). Three different datasets were compiled in attempts to minimize the amount of missing data: 16 taxa (reduced taxon set to include the taxa with the most complete data for each major vertebrate group), 31 genes (for turtle question, genes with representatives from all the major groups in question), 26 genes (for lissamphibian question, genes with representatives from all the major groups in question). The last column summarizes the results from phylogenetic analyses when removing the 19 most quickly evolving genes (25% of total genes).</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>For the lissamphibian question, phylogenetic analyses for N12, DEGEN1, and AA often recovered different relationships, but most AU tests did not exclude any of the four major hypotheses. The NUCL data-type was unique in that the recovered topology had no two lissamphibian groups monophyletic (<xref ref-type="supplementary-material" rid="pone.0048990.s001">Figure S1</xref>
), but of the four major hypotheses, provided the most support for a paraphyletic Lissamphibia, with a caecilian-salamander clade (<xref ref-type="table" rid="pone-0048990-t001">Table 1</xref>
).</p>
<p>Results from phylogenetic analyses and AU tests performed on the ‘slow genes’ dataset are summarized in <xref ref-type="table" rid="pone-0048990-t001">Table 1</xref>
. For the turtle question, a turtle-crocodile relationship was recovered, which was statistically indistinguishable from a turtle-archosaur relationship. For the lissamphibian question, results were identical to the NUCL data-type where no two lissamphibian groups were monophyletic, and AU tests could not statistically reject any of the four hypotheses.</p>
</sec>
<sec id="s2c"><title>Rogue Taxa Analyses</title>
<p>Unstable (rogue) taxa in a phylogeny can affect phylogenetic inference. Removal of these taxa can improve phylogenetic results by increasing resolution and/or support values <xref ref-type="bibr" rid="pone.0048990-Pattengale1">[33]</xref>
. We identified 19–39 rogue taxa for each of the four data-types, with much overlap between data-types. Although phylogenetic relationships of major groups were the same, removal of rogue taxa improved analyses by increasing bootstrap support values of clades.</p>
</sec>
<sec id="s2d"><title>Statistical Analyses</title>
<p>Initial phylogenetic results were inconclusive, possibly due to conflicts between phylogenetic and non-phylogenetic signal. Features of the data that may be correlated with biases in phylogenetic reconstruction include site-specific rates of evolution (site-rates), as well as heterogeneities between clades in GC content (%GC) and amount of missing data (%missing) <xref ref-type="bibr" rid="pone.0048990-RodriguezEzpeleta1">[45]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Lemmon1">[46]</xref>
. We reason that if these correlates of non-phylogenetic signal alone can do a good job of predicting the phylogeny favored by a site in the alignment, this site is likely to be biased and cannot be trusted. A diagram of our methodology to identify biased sites can be found in <xref ref-type="fig" rid="pone-0048990-g003">Figure 3</xref>
. First, we compute site-rates for each site in the alignment, and %GC and %missing per site for major clades relevant to turtle placement and lissamphibian relationships. In addition, we compute site-wise likelihoods for all competing hypotheses regarding the phylogenetic positions of turtles and Lissamphibia and recorded the topology with the highest likelihood for each site. Next, we use Discriminant Function Analysis (DFA; employing a quadratic discriminant function) to predict the favored topology based solely on descriptive statistics (site-rates, %GC, and %missing). Based on the strength with which the DFA was able to predict the topology preferred by any site, we designated sites as putatively biased and progressively removed them from the analysis.</p>
<fig id="pone-0048990-g003" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0048990.g003</object-id>
<label>Figure 3</label>
<caption><title>Flow diagram of data filtering method.</title>
<p>Steps of the new statistical methodology to identify and filter out sites that contain putative non-phylogenetic signal (i.e. biased sites). Analyses pertaining to the phylogenetic position of turtles are used in this example.</p>
</caption>
<graphic xlink:href="pone.0048990.g003"></graphic>
</fig>
<p>We validated our approach on simulated data and a previously published biological dataset <xref ref-type="bibr" rid="pone.0048990-Wolf1">[47]</xref>
. We simulated sequences along a tree with 8 leaves under strong heterogeneities in rates of evolution among sites, relative branch lengths among sites, and equilibrium %GC among taxa (see methods). Phylogenetic reconstruction using all sites without filtering resulted in an artifactual topology where species with similar %GC and high rates of evolution clustered together. We computed %GC and site-rates (simulated dataset was complete, without missing data) on the simulated sequences, and used our procedure to filter sites using site-rates only or both site-rates and %GC. Although the sequences had been simulated under strong compositional heterogeneity, filtering based on site-rates only resulted in better overall results. In fact, the removal of putatively biased sites resulted in the recovery of the correct topology at all thresholds tested. In contrast, filtering based on both %GC and site-rates resulted in the recovery of the correct topology only when removing the largest proportion of sites (<xref ref-type="supplementary-material" rid="pone.0048990.s004">Table S1</xref>
). Contrary to site-rates, %GC contains a complex mixture of phylogenetic and biased signal, which may confuse the method, as shown by the following toy example. If one considers 3 clades A, B, C, with the correct topology ((A,B),C) and convergence towards higher GC content in clades B and C leads to the artifactual topology (A,(B,C)). High GC contents in clade B, clade C, or even in both clades B and C are not by themselves sufficient for predicting that a site is likely to provide biased signal. Only in the case where A is GC poor and both B and C are GC rich can this site be safely assumed to likely provide biased signal. All seven other configurations (all three clades GC rich; the two other configurations with two clades GC rich, the three configurations where two clades are GC poor, and three clades GC poor) are not indicative of a compositional artifact. Consequently, to predict putatively biased sites using compositional statistics for clades, a complex interaction between three variables has to be uncovered by the method. As our DFA does not consider interaction terms between two or more variables, it cannot perform well with %GC. Other predictor variables (e.g., site-rates or %missing) may not require interactions between two or more variables for predicting putatively biased sites, and are thus more amenable to our analysis through DFA. For instance, the rate of a site or the percent of missing data in a particular clade could be enough to predict that a site has the potential for providing biased signal.</p>
<p>To further validate our approach, we used a dataset of eight gene concatenates addressing the Ecdysozoa-Coelomata controversy <xref ref-type="bibr" rid="pone.0048990-Aguinaldo1">[48]</xref>
. In their paper, Wolf et al. (2003) <xref ref-type="bibr" rid="pone.0048990-Wolf1">[47]</xref>
concluded in favor of the Coelomata hypothesis, as analyses of the datasets resulted in 5/8 topologies strongly supporting Coelomata. However, most recent studies support the Ecdysozoa hypothesis and suggest that the Coelomata hypothesis is an artifactual result linked to fast-evolving taxa and inadequate taxonomic sampling <xref ref-type="bibr" rid="pone.0048990-Telford1">[49]</xref>
. The original dataset was a complete, amino acid dataset, so we are unable to calculate %GC and %missing. Therefore, we only computed site-rates and applied our filtering procedure on the eight datasets, comparing it to random removal of sites as a control. After filtering, 6/8 alignments support the Ecdysozoa hypothesis (<xref ref-type="supplementary-material" rid="pone.0048990.s005">Table S2</xref>
), changing the support of three genes from Coelomata to Ecdysozoa. These results suggest that our approach had successfully filtered out biased signal from the alignments.</p>
<p>We computed site-wise descriptive statistics and most likely topologies for the NUCL dataset. As our method focuses on specific phylogenetic questions, we performed filtering of biased sites twice, once for the turtle question and once for the Lissamphibia question, producing two different sets of alignments.</p>
<p>We find that DFA accurately predicts the most likely topology for 47% of the sites for Lissamphibia, and 36% of the sites for turtles. DFA is able to predict the topology with the highest site likelihood more accurately than the control (see methods; sites are correctly predicted by the DFA analysis 1.55× and 1.65× more often than random expectations for Lissamphibia and turtles, respectively) (<xref ref-type="supplementary-material" rid="pone.0048990.s006">Table S3</xref>
). The predictive ability of DFA is significantly better than expected at random, based on the results of permutation tests (<xref ref-type="supplementary-material" rid="pone.0048990.s003">Figure S3</xref>
).</p>
<p>Interestingly, the ability of DFA to predict the preferred topology at a site varies by topology. In lissamphibians, DFA is most able to predict the Procera topology (1.98× more accurately than the control predictor) and least able to predict the Batrachia topology (1.43×). In turtles, DFA is most able to predict the Lepidosaur topology (3.62×) and least able to predict the Archosaur topology (0.66×) (<xref ref-type="supplementary-material" rid="pone.0048990.s006">Table S3</xref>
). For each site, DFA can also be used to calculate a support value corresponding to the strength of its prediction. For instance, regarding Lissamphibia, the 1% most confidently predicted sites based on DFA all support the Procera hypothesis, and for turtles, the 1% most confidently predicted sites all support the Sauropsid topology. This shows that the Procera topology for lissamphibian relationships, and the Lepidosaur and Sauropsid topologies for turtle placement can be predicted by characteristics of the sites that should be unrelated to the site’s preferred topology, and suggests they may be supported in part by non-phylogenetic signal in the alignment. We note that all four candidate topologies for Lissamphibia are predicted with similar accuracies by the DFA analysis, in contrast with the turtle analysis. This may imply that the biased signal we detect is more equally distributed among the different lissamphibian hypotheses than for the turtle hypotheses.</p>
<p>Based on the performance of DFA-filtering when analyzing simulated as well as empirical data, we performed two DFA analyses: three types of descriptive statistics (site-rates, %GC, and %missing) or two types (excluding %GC). We generated several alignments by removing the 10%, 20%, 30%, 40%, or 50% most confidently predicted (i.e. most suspect) sites from the alignment for the turtle and Lissamphibia analyses, and generated phylogenies from these sub-sampled alignments as well as alignments of the discarded sites. For turtles, all phylogenetic analyses and topology tests based on DFA-filtering using all three descriptive statistics support turtles as the sister group to crocodilians (<xref ref-type="supplementary-material" rid="pone.0048990.s007">Table S4</xref>
). Filtered datasets generated without the use of clade-specific %GC as a predictor supported either turtle-crocodilian or turtle-archosaur relationships (<xref ref-type="supplementary-material" rid="pone.0048990.s007">Table S4</xref>
). For Lissamphibia, all analyses using all three descriptive statistics support the same topology in which Lissamphibians are paraphyletic and a caecilian-salamander clade forms the sister group to amniotes (<xref ref-type="supplementary-material" rid="pone.0048990.s007">Table S4</xref>
). However, for analyses excluding %GC, two hypotheses (Procera and Paraphyletic Caecilian-Salamander) are often statistically indistinguishable. Additionally, when excluding %GC and removing 50% and 40% of the data, supported topologies do not match any of the four proposed hypotheses (<xref ref-type="supplementary-material" rid="pone.0048990.s007">Table S4</xref>
). The low bootstrap support values suggest these highly unlikely topologies come from an absence of a clear phylogenetic signal in the remaining sites.</p>
<p>From the four alignments with the 10% most suspect data removed, one for each combination of taxonomic question and number of DFA predictor types (2 or 3), we can exclude all but four possible topologies relating major vertebrate groups. We combine these trees to produce a consensus phylogeny, with relationships within amniotes from the turtle datasets and deeper vertebrate relationships from the lissamphibian datasets. The consensus phylogeny of higher-level vertebrate relationships from our study is in <xref ref-type="fig" rid="pone-0048990-g004">Figure 4</xref>
.</p>
<fig id="pone-0048990-g004" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0048990.g004</object-id>
<label>Figure 4</label>
<caption><title>Consensus vertebrate phylogeny.</title>
<p>Consensus phylogeny from datasets with the 10% most putatively biased sites removed. (A) Turtles are either the sister group to Crocodilians or Archosauria. (B) Lissamphibia: salamanders (Caudata) and caecilians (Gymnophiona) are sister groups, and this group is either the sister group to frogs (Procera hypothesis) or Amniota (rendering Lissamphibia paraphyletic). RAxML bootstrap values are at nodes, with “*” representing support ≥95.</p>
</caption>
<graphic xlink:href="pone.0048990.g004"></graphic>
</fig>
</sec>
</sec>
<sec id="s3"><title>Discussion</title>
<p>Previous studies of the vertebrate phylogeny have resulted in ambiguity regarding the phylogenetic placement of turtles within amniotes and the interrelationships within Lissamphibia (<xref ref-type="fig" rid="pone-0048990-g001">Figure 1</xref>
), in part because the short internodes and long branches that characterize these groups are notoriously difficult problems in phylogenetic inference. Using standard phylogenetic approaches, past studies – as well as similar efforts with our data – have not yielded consistent results (see <xref ref-type="bibr" rid="pone.0048990-Thomson1">[5]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Carroll2">[12]</xref>
). We believe that difficult phylogenetic problems, such as these, could be due to the presence of conflicting phylogenetic signal in the dataset. In large datasets, the problem may not be the amount of phylogenetic signal, but rather the confounding effects of phylogenetic error. Philippe et al. (2011) <xref ref-type="bibr" rid="pone.0048990-Philippe1">[31]</xref>
outline three primary sources of phylogenetic error: 1) incorrect identification of orthologs, 2) erroneous sequence alignments, and 3) inadequate models of evolution. The first two points are addressed in our dataset by rigorously testing orthology and alignment through the marker development and data analysis stages <xref ref-type="bibr" rid="pone.0048990-Fong1">[41]</xref>
. Some standard methods to address the third point are to reduce homoplasy by transforming data and removing genes. For our study, data transformations were ineffective at removing conflicting signal, while removal of fast evolving genes was partially successful, but conflicting signal remained, especially for the lissamphibian question. Accordingly we developed a new method that predicts and removes potentially biased sites for a specific phylogenetic question. Our method tests the potential for biased inferences to result from a high rate of evolution as well as two other potential contributors of non-phylogenetic signal <xref ref-type="bibr" rid="pone.0048990-RodriguezEzpeleta1">[45]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Lemmon1">[46]</xref>
: GC content (%GC) and proportion of missing data (%missing). Our simulations and tests on empirical data showed that our approach is promising in its ability to remove biased signal, notably when %GC is not included as a predictive variable. When we implemented this statistical procedure to filter our data, we reduced conflicting signal and recovered stronger support for higher-level vertebrate relationships.</p>
<sec id="s3a"><title>Phylogenetic Position of Turtles</title>
<p>Past studies have hypothesized five different phylogenetic positions for turtles in the amniote phylogeny, with the most recent molecular studies debating between the turtle-lepidosaur <xref ref-type="bibr" rid="pone.0048990-Lyson1">[7]</xref>
and turtle-archosaur <xref ref-type="bibr" rid="pone.0048990-Crawford1">[8]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Chiari1">[9]</xref>
relationships. Removing the set of sites identified by DFA to have the greatest chance of contributing biased signal allowed statistical exclusion of three previously proposed hypotheses: a turtle-lepidosauria sister grouping, turtles as basal sauropids (reptiles and birds), and turtles as basal amniotes (<xref ref-type="table" rid="pone-0048990-t002">Table 2</xref>
). Our results show that turtles are closely related to birds and crocodilians, but since results differed when clade-specific %GC content was or was not included in the set of predictor variables (<xref ref-type="table" rid="pone-0048990-t002">Table 2</xref>
), we are not able to distinguish between the turtle-archosaur and turtle-crocodilian topologies.</p>
<table-wrap id="pone-0048990-t002" orientation="portrait" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0048990.t002</object-id>
<label>Table 2</label>
<caption><title>Phylogenetic results from filtered datasets.</title>
</caption>
<alternatives><graphic id="pone-0048990-t002-2" xlink:href="pone.0048990.t002"></graphic>
<table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
<col align="center" span="1"></col>
</colgroup>
<thead><tr><td align="left" rowspan="1" colspan="1"></td>
<td colspan="5" align="left" rowspan="1">3 types (site-rates, %GC, %missing)</td>
<td colspan="5" align="left" rowspan="1">2 types (site-rates, %missing)</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">TURTLE</td>
<td align="left" rowspan="1" colspan="1">50%</td>
<td align="left" rowspan="1" colspan="1">40%</td>
<td align="left" rowspan="1" colspan="1">30%</td>
<td align="left" rowspan="1" colspan="1">20%</td>
<td align="left" rowspan="1" colspan="1">10%</td>
<td align="left" rowspan="1" colspan="1">50%</td>
<td align="left" rowspan="1" colspan="1">40%</td>
<td align="left" rowspan="1" colspan="1">30%</td>
<td align="left" rowspan="1" colspan="1">20%</td>
<td align="left" rowspan="1" colspan="1">10%</td>
</tr>
</thead>
<tbody><tr><td align="left" rowspan="1" colspan="1">Archosaur</td>
<td align="left" rowspan="1" colspan="1">7e-5</td>
<td align="left" rowspan="1" colspan="1">2e-5</td>
<td align="left" rowspan="1" colspan="1">1e-4</td>
<td align="left" rowspan="1" colspan="1">1e-4</td>
<td align="left" rowspan="1" colspan="1">0.031</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.313</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.173</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Crocodilian</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1">0.006</td>
<td align="left" rowspan="1" colspan="1"><bold>0.119</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.32</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Lepidosaur</td>
<td align="left" rowspan="1" colspan="1">6e-79</td>
<td align="left" rowspan="1" colspan="1">5e-33</td>
<td align="left" rowspan="1" colspan="1">8e-98</td>
<td align="left" rowspan="1" colspan="1">2e-54</td>
<td align="left" rowspan="1" colspan="1">6e-6</td>
<td align="left" rowspan="1" colspan="1">1e-68</td>
<td align="left" rowspan="1" colspan="1">1e-7</td>
<td align="left" rowspan="1" colspan="1">2e-59</td>
<td align="left" rowspan="1" colspan="1">1e-46</td>
<td align="left" rowspan="1" colspan="1">1e-48</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Basal Sauropsid</td>
<td align="left" rowspan="1" colspan="1">4e-5</td>
<td align="left" rowspan="1" colspan="1">3e-117</td>
<td align="left" rowspan="1" colspan="1">4e-5</td>
<td align="left" rowspan="1" colspan="1">2e-8</td>
<td align="left" rowspan="1" colspan="1">2e-4</td>
<td align="left" rowspan="1" colspan="1">8e-65</td>
<td align="left" rowspan="1" colspan="1">2e-50</td>
<td align="left" rowspan="1" colspan="1">2e-70</td>
<td align="left" rowspan="1" colspan="1">1e-56</td>
<td align="left" rowspan="1" colspan="1">1e-44</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Basal Amniote</td>
<td align="left" rowspan="1" colspan="1">6e-5</td>
<td align="left" rowspan="1" colspan="1">7e-10</td>
<td align="left" rowspan="1" colspan="1">1e-6</td>
<td align="left" rowspan="1" colspan="1">1e-84</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
<td align="left" rowspan="1" colspan="1">5e-26</td>
<td align="left" rowspan="1" colspan="1">1e-4</td>
<td align="left" rowspan="1" colspan="1">5e-7</td>
<td align="left" rowspan="1" colspan="1">1e-6</td>
<td align="left" rowspan="1" colspan="1">12e-5</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"></td>
<td colspan="5" align="left" rowspan="1"><bold>3 types (site-rates, %GC, %missing)</bold>
</td>
<td colspan="5" align="left" rowspan="1"><bold>2 types (site-rates, %missing)</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1"><bold>LISSAMPHIBIA</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>50%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>40%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>30%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>20%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>10%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>50%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>40%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>30%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>20%</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>10%</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Batrachia</td>
<td align="left" rowspan="1" colspan="1">7e-36</td>
<td align="left" rowspan="1" colspan="1">7e-34</td>
<td align="left" rowspan="1" colspan="1">3e-6</td>
<td align="left" rowspan="1" colspan="1">2e-65</td>
<td align="left" rowspan="1" colspan="1">4e-52</td>
<td align="left" rowspan="1" colspan="1">0.004</td>
<td align="left" rowspan="1" colspan="1">0.003</td>
<td align="left" rowspan="1" colspan="1">0.019</td>
<td align="left" rowspan="1" colspan="1">0.015</td>
<td align="left" rowspan="1" colspan="1">0.003</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Procera</td>
<td align="left" rowspan="1" colspan="1">7e-11</td>
<td align="left" rowspan="1" colspan="1">5e-61</td>
<td align="left" rowspan="1" colspan="1">5e-5</td>
<td align="left" rowspan="1" colspan="1">2e-66</td>
<td align="left" rowspan="1" colspan="1">93–94</td>
<td align="left" rowspan="1" colspan="1">0.029</td>
<td align="left" rowspan="1" colspan="1"><bold>0.405</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.375</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>0.131</bold>
</td>
<td align="left" rowspan="1" colspan="1">3e-4</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Paraphyletic</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Frog-Salamander</td>
<td align="left" rowspan="1" colspan="1">3e-76</td>
<td align="left" rowspan="1" colspan="1">2e-7</td>
<td align="left" rowspan="1" colspan="1">1e-6</td>
<td align="left" rowspan="1" colspan="1">6e-67</td>
<td align="left" rowspan="1" colspan="1">5e-6</td>
<td align="left" rowspan="1" colspan="1">2e-7</td>
<td align="left" rowspan="1" colspan="1">0.004</td>
<td align="left" rowspan="1" colspan="1">1e-4</td>
<td align="left" rowspan="1" colspan="1">0.002</td>
<td align="left" rowspan="1" colspan="1">0.001</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Paraphyletic</td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
<td align="left" rowspan="1" colspan="1"></td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Caecilian-Salamander</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1">7e-11</td>
<td align="left" rowspan="1" colspan="1">0.036</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
</tr>
<tr><td align="left" rowspan="1" colspan="1">Other</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1"><bold>X</bold>
</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
<td align="left" rowspan="1" colspan="1">–</td>
</tr>
</tbody>
</table>
</alternatives>
<table-wrap-foot><fn id="nt102"><label></label>
<p>Varying amounts of suspect sites were removed and tested. A) Position of turtles in the amniote phylogeny using three descriptive statistics (site-rates, %GC, and %missing), B) position of turtles in the amniote phylogeny using two descriptive statistics (excluding %GC), C) interrelationships of Lissamphibian groups using three descriptive statistics (site-rates, %GC, and %missing), D) interrelationships of Lissamphibian groups using two descriptive statistics (excluding %GC). The percentage in each column represents the percentage of sites removed from the dataset. Values in cells represent p-values, “X” denotes the best tree, and trees statistically indistinguishable from the best tree are in bold font (Approximately Unbiased topology test p-value >5%).</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Recent results <xref ref-type="bibr" rid="pone.0048990-Crawford1">[8]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Chiari1">[9]</xref>
, as well as our findings, placing turtles as close relatives of crocodilians and birds, necessitates changing the traditional view of turtle evolution, as it prevailed until recently. First, Archosauria is defined as the crown group including the most recent common ancestor of birds and crocodilians <xref ref-type="bibr" rid="pone.0048990-Gauthier1">[50]</xref>
; turtles are either the sister group to or a member of Archosauria. A recent paleontological study supports this relationship between archosaurs and turtles, discovering a unique skull ossification (laterosphenoid) found only in turtles and Archosauriformes <xref ref-type="bibr" rid="pone.0048990-Bhullar1">[51]</xref>
.</p>
<p>An important question in turtle biology is how and when its unique, shelled body plan evolved. Previous work suggested parareptilian (Anapsida) groups as the extinct ancestor of turtles <xref ref-type="bibr" rid="pone.0048990-Laurin1">[52]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Lee1">[53]</xref>
, with one hypothesis pointing towards the elaboration of dermal armor as a precursor to formation of the shell <xref ref-type="bibr" rid="pone.0048990-Lee2">[54]</xref>
. With our results nesting turtles among diapsids, hypotheses of turtle shell evolution from parareptilian ancestors are no longer possible. Turtles are unique in that their ribs develop by encapsulating the shoulder blades and embed within the dermis, sending developmental signals to the dermis to form bone and therefore the carapace <xref ref-type="bibr" rid="pone.0048990-Nagashima1">[55]</xref>
. Understanding how and when the turtle shell arose will come only from studying extinct archosaurian lineages.</p>
</sec>
<sec id="s3b"><title>Relationships within Lissamphibia</title>
<p>Ancestral amphibians appear in the fossil record starting in the late Devonian and are extremely diverse in the Palaeozoic. However, a large gap in the fossil record exists between Palaeozoic amphibians and lissamphibians, with the exception of Stereospondyls extending into the Mesozoic <xref ref-type="bibr" rid="pone.0048990-Yates2">[56]</xref>
and a possible frog ancestor found in the Lower Triassic <xref ref-type="bibr" rid="pone.0048990-Carroll2">[12]</xref>
. It is this gap in the fossil record paired with significant morphological change that has made it difficult to determine the ancestors of and relationships among modern amphibians from paleontological data.</p>
<p>The most recent molecular study based on mitochondrial genomes and eight nuclear genes <xref ref-type="bibr" rid="pone.0048990-SanMauro2">[24]</xref>
supports the Batrachia hypothesis, which is in contrast to most paleontological studies supporting a paraphyletic Lissamphibia <xref ref-type="bibr" rid="pone.0048990-Carroll2">[12]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Cannatella1">[13]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Anderson1">[25]</xref>
–<xref ref-type="bibr" rid="pone.0048990-Carroll5">[29]</xref>
. This raises the question of whether morphological or molecular data are correct <xref ref-type="bibr" rid="pone.0048990-Anderson2">[26]</xref>
. The results of our molecular study supports a caecilian-salamander sister relationship, but cannot distinguish between a monophyletic and paraphyletic Lissamphibia (<xref ref-type="fig" rid="pone-0048990-g004">Figure 4</xref>
). Although our results do not resolve the lissamphibian origin problem, we resurrect an often-overlooked hypothesis of a paraphyletic Lissamphibia. This is the first molecular study to have signal supporting a paraphyletic Lissamphibia, allowing for the possibility of concordance between morphological and molecular data. If the paraphyletic hypothesis is true, caecilians and salamanders would be more closely related to humans and other amniotes than to frogs. To discriminate between the alternative hypotheses of a monophyletic and paraphyletic Lissamphibia, deeper taxon and gene sampling for lissamphibian groups will probably be needed.</p>
</sec>
<sec id="s3c"><title>Removal of Biased Signal</title>
<p>A common notion in molecular systematics is that the solution to resolving difficult relationships is to include ever increasing amounts of data. This belief is based on the idea that true phylogenetic signal will eventually dominate and drive the results of an analysis, circumventing any methodological problems. However, our results suggest that inclusion of more data can introduce biased signal into a dataset, resulting in a lack of resolution or even misleading inferences, a possibility also raised by others <xref ref-type="bibr" rid="pone.0048990-Philippe1">[31]</xref>
. Methods to filter data based on the rate of evolution have been previously used to increase resolution by removing non-phylogenetic signal <xref ref-type="bibr" rid="pone.0048990-Brinkmann1">[39]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Philippe2">[40]</xref>
. In an effort to remove non-phylogenetic signal in a more targeted way, we develop a new DFA-based filtering method that attempts to identify sites contributing biased signal based on several data set characteristics known to cause inference problems in certain contexts (site-rates, %GC, %missing). In addition to using more information for predicting non-phylogenetic signal, our approach is different from other methods because it targets a specific phylogenetic relationship (i.e. sites identified as biased for the turtle question are specific to that question). We anticipate and encourage future studies that will more fully assess which characteristics of sites can be used with this method to accurately predict the presence of phylogenetic bias (e.g. whether or not %GC can improve the prediction of bias and, if so, when). Future studies could also assess potential performance gains from including interaction terms among variables in the DFA or including information about the strength with which a particular site supports its preferred phylogenetic hypothesis.</p>
<p>Proper modeling of molecular evolution and evaluation of the fit between data and model seem to be just as important as the amount of data present in a study. With the advent of new technologies that produce sequence data faster and more cheaply than ever before, datasets will only become larger, and issues relating to signal quality will become even more important in molecular systematics. We view our DFA approach as an important step towards the goal of objectively identifying non-phylogenetic signal in large datasets.</p>
</sec>
<sec id="s3d"><title>Conclusion</title>
<p>With increasingly large datasets being gathered for phylogenetics, many relationships have been confidently resolved. What remains are controversial, difficult to resolve phylogenetic questions, probably arising from conflicting and biased phylogenetic signal in the data. We developed a method for identifying and minimizing biases from molecular data to tackle two persistent yet fundamental problems in vertebrate phylogenetics: the placement of turtles within amniotes and the interrelationships within Lissamphibia. Based on tests of our filtering method on simulated and empirical datasets, we believe that we are able to reduce the amount of conflicting signal in datasets. For the vertebrate phylogeny, the application of this filtering method results in analyses that support turtles being closely related to archosaurs, as either the sister group to crocodilians or archosaurs, and a caecilian-salamander sister relationship, with the possible paraphyly of Lissamphibia. Because of our use of a new statistical approach, we view our results to be tentative and encourage more work from paleontologists and molecular biologists alike to further evaluate these hypotheses and methodology. Given the importance of the historical framework provided by phylogenetic systematics in fields as diverse as developmental biology, genomics, conservation biology and paleontology, we believe approaches like ours will be useful to resolve major phylogenetic questions and advance modern biological thought.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="s4"><title>Materials and Methods</title>
<sec id="s4a"><title>Ethics Statement</title>
<p>This research was conducted under and approved by UC Berkeley's Animal Care and Use Committee (protocol #R279-0211). Tissue samples used in this study were obtained from the Museum of Vertebrate Zoology (MVZ), an institution that serves as a specimen and tissue repository for researchers. The MVZ has a strict policy for researchers when depositing specimens and tissues into the museum, requiring local collecting permits and import permits when necessary.</p>
</sec>
<sec id="s4b"><title>Taxon Sampling</title>
<p>Our study included sampling for all major vertebrate groups except Tuatara (<italic>Sphenodon</italic>
). Omission of the Tuatara is inconsequential to our investigation due to its uncontroversial affinity with Squamata ( = Lepidosauria) <xref ref-type="bibr" rid="pone.0048990-Hedges2">[57]</xref>
. Of the 129 taxa included in this study, data were available for 46 taxa from complete genomes or ESTs (GenBank and Ensembl <xref ref-type="bibr" rid="pone.0048990-Hubbard1">[58]</xref>
), and the remaining 82 samples were newly sequenced. Available data were skewed heavily towards mammals (36 taxa), but also include fish (5 taxa), frog (1 taxon), salamander (1 taxon), lizard (1 taxon), bird (2 taxa), and crocodile (1 taxon). New taxon sampling expanded representation within amphibians, reptiles, and birds. A detailed list of all taxa used in this study, along with Genbank numbers, can be found in <xref ref-type="supplementary-material" rid="pone.0048990.s008">Table S5</xref>
.</p>
</sec>
<sec id="s4c"><title>Marker Sampling</title>
<p>Markers used in this study are single-copy, orthologous, protein-coding genes <xref ref-type="bibr" rid="pone.0048990-Fong1">[41]</xref>
. The single-copy nature of the markers was checked both during marker development and after data collection. During marker development, markers were pre-screened using BLAST to compare with the high-coverage, well-annotated chicken genome and gene families were identified in Ensembl and Metazome. After data collection, to identify and remove paralogous genes, preliminary phylogenetic analyses using RAxML <xref ref-type="bibr" rid="pone.0048990-Stamatakis1">[59]</xref>
were performed for each gene. Each of the trees was analyzed individually by eye for erroneous phylogenetic relationships (e.g. mammal species more closely related to fish) and signatures of gene duplications (i.e. replicated tree topologies within the larger tree). Sequences with erroneous positions were removed, and when gene duplications were detected, the sub-tree that included newly sequenced data was retained for subsequent analyses. We obtained new sequences according to the methods of Fong & Fujita (2011) <xref ref-type="bibr" rid="pone.0048990-Fong1">[41]</xref>
. Briefly, we used cDNA preparations as template to amplify the protein-coding genes using conserved primers. Amplicons were sequenced using ABI3730 chemistry, and sequences were edited using Geneious 5 (Biomatters Ltd.) and aligned using MUSCLE <xref ref-type="bibr" rid="pone.0048990-Edgar1">[60]</xref>
.</p>
</sec>
<sec id="s4d"><title>Datasets and Data-types</title>
<p>Sequences were combined into two main categories of datasets: individual genes and concatenations. Individual datasets for the 75 genes consisted of orthologous sequences from online genomes and the new samples. Combining individual genes using a Perl script (available upon request) produced the concatenated datasets. We compiled seven different concatenated datasets: 1) All taxa-75 genes, 2) 16 taxa-75 genes, 3) All taxa-31 genes (turtle), 4) All taxa-26 genes (Lissamphibia), 5) 16 taxa-31 genes (turtle), 6) 16 taxa-26 genes (Lissamphibia), 7) slow genes (removal of fastest 25% of genes). For the 16-taxon datasets, the vertebrate group Crocodilia is represented by two individuals of the species <italic>Alligator mississippiensis</italic>
(from the EST database and a new sample). We combined the data from both individuals to minimize missing data; this approach is justified, as when there were data from both individuals for a marker, data were identical. The dataset of reduced loci for all taxa was used when evaluating the specific phylogenetic questions (turtle and Lissamphibia). Loci without representatives of all the focal groups were removed, leaving 31 genes for the turtle analysis and 26 genes for the lissamphibian analysis.</p>
<p>The standard nucleotide (NUCL) dataset was transformed to three data-types using the following methods. AA was translated in Geneious 5 (Biomatters Ltd.), 3<sup>rd</sup>
codon positions were removed for N12 using MacClade v4.08 <xref ref-type="bibr" rid="pone.0048990-Maddison1">[61]</xref>
, and DEGEN1 was converted using a Perl script <xref ref-type="bibr" rid="pone.0048990-Regier1">[38]</xref>
.</p>
<p>The rate of evolution for each of the 75 genes was calculated by computing tree length and averaging branch lengths using an R script <xref ref-type="bibr" rid="pone.0048990-R1">[62]</xref>
. Based on the shape of the frequency histogram (<xref ref-type="supplementary-material" rid="pone.0048990.s002">Fig. S2</xref>
), we drew a cut-off of average branch length of 0.04, which denoted the top 25% fastest genes (19 genes) for removal. The names of these genes are: DSCR3, EXOC2, GAPDH, GDE1, GNAS, HPD, H2AFY2, IFT57, MAT2B, OAT, OAZ1, PPM1A, PPTC7, PSAT1, SEC13, SGK1, TAT, UBE2J2, and XPOT (<xref ref-type="supplementary-material" rid="pone.0048990.s008">Table S5</xref>
) <xref ref-type="bibr" rid="pone.0048990-Fong1">[41]</xref>
.</p>
</sec>
<sec id="s4e"><title>Phylogenetic Analyses</title>
<p>All datasets were subject to maximum likelihood analyses using RAxML <xref ref-type="bibr" rid="pone.0048990-Stamatakis1">[59]</xref>
, and a subset of datasets were also subject to Bayesian analyses using MrBayes <xref ref-type="bibr" rid="pone.0048990-Ronquist1">[63]</xref>
.</p>
<p>Since this study deals with a complex, multi-gene dataset, we explored heterogeneous processes of molecular evolution through partitioning the data. Tests of alternative partitioning strategies were performed on the NUCL dataset only, as the N12 dataset is a subset of the NUCL dataset, and the DEGEN1 and AA datasets have information on codon position integrated into gene partitions. For RAxML analyses of the NUCL dataset, three different partitioning strategies were tested: by gene (75 partitions), by gene and 1<sup>st</sup>
+2<sup>nd</sup>
and 3<sup>rd</sup>
codon position (150 partitions), and by gene and codon position (225 partitions). Likelihood ratio tests selected the 150 partitions as the best partitioning strategy.</p>
<p>RAxML nucleotide analyses used the GTRGAMMA model of evolution for tree inference and bootstrapping (1,000 replicates). RAxML amino acid analyses employed the protein gamma model of evolution and the appropriate model of protein evolution selected using ProtTest v2.4 <xref ref-type="bibr" rid="pone.0048990-Abascal1">[64]</xref>
, with empirical amino-acid frequencies for both tree inference and bootstrapping (1,000 replicates). All concatenated datasets were partitioned according to the optimal partitioning strategy. RAxML v7.2.5 and v7.2.6 <xref ref-type="bibr" rid="pone.0048990-Stamatakis1">[59]</xref>
, <xref ref-type="bibr" rid="pone.0048990-Stamatakis2">[65]</xref>
analyses were run locally and on the CIPRES portal v2.2 <xref ref-type="bibr" rid="pone.0048990-Miller1">[66]</xref>
.</p>
<p>For individual gene analyses, the support values of clades are generally very low, since these relatively short genes (average length is ∼450bp) were used to infer the entire vertebrate phylogeny. However, to understand and summarize the phylogenetic signal for each gene, we classified them based on preferred topology (see <xref ref-type="fig" rid="pone-0048990-g002">Figure 2</xref>
) irrespective of nodal support.</p>
<p>Bayesian analyses were only run on individual genes and 16-taxon datasets, as the computational burden for the larger datasets would require extremely long analysis times to achieve stationarity (i.e. >2000 hours). When both RAxML and MrBayes analyses were run, preferred topologies were almost identical, so results should not be compromised by reporting only inferences from RAxML. MrBayes v3.1.2 and v3.2 <xref ref-type="bibr" rid="pone.0048990-Ronquist1">[63]</xref>
analyses were run locally and using the <email>BioHPC@CBSU</email>
resource at Cornell University (<ext-link ext-link-type="uri" xlink:href="http://cbsuapps.tc.cornell.edu">http://cbsuapps.tc.cornell.edu</ext-link>
). All analyses were run with four chains for 10 million generations. Appropriate models of DNA substitution for each partition were selected using MrModeltest v2.3 <xref ref-type="bibr" rid="pone.0048990-Nylander1">[67]</xref>
, and amino acid substitution models were the same as those used in RAxML analyses. Burn-in of MCMC chains was evaluated using the online program AWTY, examining cumulative plots of posterior probabilities of the 20 most variable splits <xref ref-type="bibr" rid="pone.0048990-Wilgenbusch1">[68]</xref>
.</p>
</sec>
<sec id="s4f"><title>Rogue Taxa Analyses</title>
<p>Rogue taxa analyses were performed using RAxML (Stamatakis 2006) and an algorithm described in Pattengale et al. (2010) <xref ref-type="bibr" rid="pone.0048990-Pattengale1">[33]</xref>
. Determining the set of rogue taxa to remove was a multi-step process that was run on the concatenated datasets, separately for each of the four data-types. First, a preliminary RAxML analysis was run with all 129 taxa. Based on this phylogeny, taxa that were obviously in incorrect phylogenetic positions (e.g. turtle placed in the mammal clade) were manually removed. These removed taxa tended to be those with the highest levels of missing data. This was repeated until all remaining taxa were placed in the correct clades. Next, the bootstrap results from the RAxML analysis were used as the input data in rogue taxa analyses. The rogue taxa analysis was run 20 times in total, 10 each for both the strict and majority rules consensus trees, using the “−r” randomization option to select a dropset. To maximize the number of taxa retained and remove the most unstable taxa, taxa were considered rogues if they were identified in ≥5 of the 10 analyses under either strict or majority rules consensus. These taxa were removed and steps repeated until there were no rogue taxa identified or no improvement in bootstrap values.</p>
</sec>
<sec id="s4g"><title>Topology Tests</title>
<p>Approximately unbiased topology tests (AU tests) <xref ref-type="bibr" rid="pone.0048990-Shimodaira1">[44]</xref>
were used to test whether sub-optimal trees were significantly worse than the maximum likelihood tree. AU tests were performed to compare the five turtle and four lissamphibian alternative hypotheses for each of the different datasets. Constrained RAxML analyses were run for each of the different topologies using the GTRGAMMA model of sequence evolution, and per-site log likelihoods calculated. These per-site log likelihoods were then input into the program CONSEL <xref ref-type="bibr" rid="pone.0048990-Shimodaira2">[69]</xref>
.</p>
</sec>
<sec id="s4h"><title>Discriminant Function Analyses (DFA)</title>
<p>Scripts for DFA analyses of the NUCL dataset were written using the R language <xref ref-type="bibr" rid="pone.0048990-R1">[62]</xref>
and rely on the SeqinR <xref ref-type="bibr" rid="pone.0048990-Charif1">[70]</xref>
and MASS <xref ref-type="bibr" rid="pone.0048990-Venables1">[71]</xref>
libraries. These scripts are available from the authors upon request. Per-site log-likelihood scores (LLS) were calculated for each of the constrained phylogenies pertaining to a relevant hypothesis (five turtle positions, four Lissamphibia hypotheses) using RAxML. Site-wise GC content (%GC) and proportion of missing data (%missing) were computed for major clades with potential sister-relationships for turtle placement (turtles, archosaurs, crocodilians, reptiles excluding turtles, and amniotes excluding turtles) and lissamphibian relationships (amniotes, caecilians, caecilians+salamanders, frogs+salamanders, and lissamphibians). Site-specific rates of evolution (site-rates) were calculated for each of the nine constrained phylogenies using HyPhy <xref ref-type="bibr" rid="pone.0048990-Pond1">[72]</xref>
under a GTR model of sequence evolution with model parameters estimated independently for each phylogenetic hypothesis. These rates were then averaged across the five turtle and four lissamphibian topologies. %GC and %missing were both calculated for each site in the NUCL data matrix, averaged across all taxa in each clade of interest.</p>
<p>DFA (from the MASS library) was run with preferred topology as the predicted category, and %GC (for relevant clades), %missing (for relevant clades), and site-rates as predictor variables in one case, and without %GC in another case. Posterior probabilities from the DFA were calculated using leave-one-out cross-validation and normalized with prior probabilities (posterior/prior ratio). The prior probability of assignment to any particular topology was simply the proportion of sites in the alignment preferring that topology. Two different types of DFA were tested to maximize the predictive power of our analysis: linear discriminant analysis (LDA), and quadratic discriminant analysis (QDA). Comparisons of average posterior/prior ratios show that QDA performed best (Turtle: LDA = 1.483, QDA = 1.693; Lissamphibian: LDA: 1.219, QDA = 1.572).</p>
</sec>
<sec id="s4i"><title>DFA Methodology Validation</title>
<p>We validate this new methodology by evaluating its performance in two situations: (1) when analyzing DNA data simulated under conditions known to cause phylogenetic inference problems, and (2) when analyzing empirical amino acid data for a challenging phylogenetic question <xref ref-type="bibr" rid="pone.0048990-Wolf1">[47]</xref>
.</p>
<p>For the simulation study, we simulated two 8-taxon datasets under conditions that cause standard phylogenetic methods to recover the incorrect phylogeny. Tree topologies were balanced and included four groups of two sister species. On each side of the short innermost branch are two sister groups, one of which has a short subtending branch and one of which has a long subtending branch. Equilibrium GC content was set to 80% for long branches and the sister groups that they subtend, while it was set to 50% for all other branches. Other parameters of the substitution matrices were equal among branches. Each simulated dataset was 1,500 bp in length and consisted of one 1000-bp subset, and one 500-bp subset. The 500-bp partition [(((A1∶0.05,A2∶0.05):0.5,(B1∶0.05,B2∶0.05):0.02):0.05,((C1∶0.05,C2∶0.05):0.02,(D1∶0.05,D2∶0.05):0.5):0.05);] showed larger differences in branch lengths than the other one [(((A1∶0.05,A2∶0.05):0.3,(B1∶0.05,B2∶0.05):0.02):0.05,((C1∶0.05,C2∶0.05):0.02,(D1∶0.05,D2∶0.05):0.3):0.05);]. Simulations were performed using bppseqgen <xref ref-type="bibr" rid="pone.0048990-Dutheil1">[73]</xref>
. Precise command lines used for simulation can be obtained from the authors.</p>
<p>The simulated dataset was filtered for biased sites with DFA comparing the true topology used in simulations and the biased topology in which clades with long branches are clustered together. We used either two descriptive statistics (%GC and site-rates) or one descriptive statistic (site-rate), removing 10%, 20%, 30%, 40% and 50% of the sites. These alignments were compared to the random removal of a comparable number of sites. Phylogenetic analyses of the datasets were performed using RAxML with the same parameters as above. Results of these analyses are found in <xref ref-type="supplementary-material" rid="pone.0048990.s004">Table S1</xref>
. Removing sites according to the strength of DFA prediction has different effects on the phylogenetic inference compared to removing sites at random: DFA-based removal allows recovery of the correct topology when enough sites are removed, but random removal does not. In this simulation, filtering based on DFA prediction excluding %GC performs better at escaping the incorrect topology and closely estimating the true phylogeny (<xref ref-type="supplementary-material" rid="pone.0048990.s004">Table S1</xref>
). This result motivated us to try DFA filtering without clade-specific %GC for our empirical data (see above).</p>
<p>For an empirical test of our methodology, we focus on the Coelomata-Ecdysozoa debate regarding metazoan phylogeny <xref ref-type="bibr" rid="pone.0048990-Aguinaldo1">[48]</xref>
. Although some studies support the traditional Coelomata relationship, evidence is mounting in support of a monophyletic Ecdysozoa (see <xref ref-type="bibr" rid="pone.0048990-Telford1">[49]</xref>
for a summary). In the multi-gene dataset of Wolf et al. (2003) <xref ref-type="bibr" rid="pone.0048990-Wolf1">[47]</xref>
, eight macromolecular complex subunits were analyzed, with 5 of 8 genes supporting the Coelomata relationship. We employ our DFA methodology on each of the eight datasets, with some slight modifications; datasets are comprised of amino acid sequences so we only use site-rates as a filtering variable (no missing data or GC content). We generated five filtered alignments varying in the amount of data removed (10%, 20%, 30%, 40%, 50%), and compare these results to random removal of sites. Results dramatically changed compared to those in the original study <xref ref-type="bibr" rid="pone.0048990-Wolf1">[47]</xref>
; originally 5/8 datasets supported the Coelomata relationship, while after data filtering, 6/8 datasets support the Ecdysozoa relationship (<xref ref-type="supplementary-material" rid="pone.0048990.s005">Table S2</xref>
). Three datasets, based on original analyses and random removal of sites, that support the Coelomata relationship (CH, CL, LF) shift support to the Ecdysozoa relationship after DFA data filtering, while three datasets originally supporting the Ecdysozoa (RI, RP) maintain support after data filtering. One of the datasets originally supporting Ecdysozoa (PR) maintains support after some filtering, but show inconsistent results for the 40% and 50% dataset (neither supporting Coelomate nor Ecdysozoa).</p>
<p>While we endeavored to include those predictor variables that we felt were most likely to be correlated with biased signal in the data, we note that these decisions were subjective and we may have left out stronger correlates. Similarly, interactions among predictors were ignored for the sake of tractability. The potential also exists for true phylogenetic signal to manifest itself in %GC in some cases, leading to the exclusion of sites with unbiased signal. However, for both analyses we repeated the QDA procedure after excluding all %GC variables as predictors and report those results as well. Examining the sensitivity of phylogenetic conclusions to these considerations will be an interesting avenue for future work.</p>
</sec>
</sec>
<sec sec-type="supplementary-material" id="s5"><title>Supporting Information</title>
<supplementary-material content-type="local-data" id="pone.0048990.s001"><label>Figure S1</label>
<caption><p><bold>Maximum Likelihood phylogenies of the different data transformations.</bold>
Phylogenies have been simplified to only show higher-level relationships within vertebrates. A) All taxa-NUCL dataset, B) 16 taxa-NUCL dataset, C) All taxa-N12 dataset, D) 16 taxa-N12 dataset (note: the Bayesian analysis recovered a turtle-crocodile relationship), E) All taxa-DEGEN1 dataset, F) 16 taxa-DEGEN1 dataset, G) All taxa-AA dataset, H)16 taxa-AA dataset. Support values for phylogenies with all taxa (A,C,E,G) show RAxML bootstrap values only if ≥50. Support values for phylogenies with 16 taxa show support values in the form of RAxML bootstrap/Bayesian posterior probability. An * indicates full support.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0048990.s001.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s002"><label>Figure S2</label>
<caption><p><bold>Frequency histogram of the rate of evolution for the 75 molecular markers.</bold>
Rate of evolution of the 75 markers was estimated by computing average branch length. The red, vertical line indicates the top 25% fastest genes, which were removed for subsequent phylogenetic analyses.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0048990.s002.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s003"><label>Figure S3</label>
<caption><p><bold>Permutation test results for the predictive ability of discriminant function analysis (DFA).</bold>
Permutation results are compared to random expectations regarding A) lissamphibian relationships and B) the phylogenetic position of turtles. Preferred lissamphibian relationships were permuted among sites (1,000 replicates). Values on the x-axis are the posterior/prior ratio for the preferred topology averaged across all sites for each replicate. The arrow indicates the empirical value, which falls far outside the null distribution.</p>
<p>(PDF)</p>
</caption>
<media xlink:href="pone.0048990.s003.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s004"><label>Table S1</label>
<caption><p><bold>Test of discriminant function analysis (DFA) filtering method on simulated data.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0048990.s004.docx" mimetype="application" mime-subtype="msword"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s005"><label>Table S2</label>
<caption><p><bold>Test of discriminant function analysis filtering method on empirical example.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0048990.s005.docx" mimetype="application" mime-subtype="msword"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s006"><label>Table S3</label>
<caption><p><bold>Predictive power of discriminant function analyses (DFA) for alternative hypotheses.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0048990.s006.docx" mimetype="application" mime-subtype="msword"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s007"><label>Table S4</label>
<caption><p><bold>Phylogenetic results from filtered datasets.</bold>
</p>
<p>(DOCX)</p>
</caption>
<media xlink:href="pone.0048990.s007.docx" mimetype="application" mime-subtype="msword"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
<supplementary-material content-type="local-data" id="pone.0048990.s008"><label>Table S5</label>
<caption><p><bold>List of Taxa and Genbank numbers.</bold>
</p>
<p>(XLS)</p>
</caption>
<media xlink:href="pone.0048990.s008.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for additional data file.</p>
</caption>
</media>
</supplementary-material>
</sec>
</body>
<back><ack><p>This work represents the partial fulfillment of the degree of Doctor of Philosophy for JJF. We would like to thank Drs. Carol Spencer and Jim McGuire, Jerome Fuchs, Chris Witt, Ricardo Pereira, and Adam Leaché for sample collection and/or tissue loans. Drs. Alexandros Stamatakis, Nick Pattengale, and Matt Brandley were invaluable in help with phylogenetic analyses. We would like to thank Jim McGuire, Jeff Boore, Richard Harland, Adam Leaché, Bryan Stuart, Guin Wogan, Shobi Lawalata, Tom Devitt, Sarah Werning, Alison Davis, Sean Reilly, and Dan Portik helpful discussions and revisions of this manuscript. Also, the comments of two anonymous reviewers also helped in improving this manuscript. Lastly, this work would not have been possible without the online computing resources provided by Cyberinfrastructure for Phylogenetic Research (CIPRES).</p>
</ack>
<ref-list><title>References</title>
<ref id="pone.0048990-Darwin1"><label>1</label>
<mixed-citation publication-type="other">Darwin C (1859) On the origin of species. London: John Murray. 502p.</mixed-citation>
</ref>
<ref id="pone.0048990-Yates1"><label>2</label>
<mixed-citation publication-type="other">Yates TL, Salazar-Bravo J, Dragoo JW (2004) The importance of the tree of life to society. In: Cracraft J, Donoghue MJ, editors. Assembling the tree of life. New York: Oxford University Press. 7–17.</mixed-citation>
</ref>
<ref id="pone.0048990-Benton1"><label>3</label>
<mixed-citation publication-type="journal"><name><surname>Benton</surname>
<given-names>MJ</given-names>
</name>
(<year>1990</year>
) <article-title>Phylogeny of the major tetrapod groups: morphological data and divergence dates</article-title>
. <source>J Mol Evol</source>
<volume>30</volume>
: <fpage>409</fpage>
–<lpage>424</lpage>
<pub-id pub-id-type="pmid">2111854</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Meyer1"><label>4</label>
<mixed-citation publication-type="journal"><name><surname>Meyer</surname>
<given-names>A</given-names>
</name>
, <name><surname>Zardoya</surname>
<given-names>R</given-names>
</name>
(<year>2003</year>
) <article-title>Recent advances in the (molecular) phylogeny of vertebrates</article-title>
. <source>Annu Rev Ecol Evol and Syst</source>
<volume>34</volume>
: <fpage>311</fpage>
–<lpage>338</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Thomson1"><label>5</label>
<mixed-citation publication-type="journal"><name><surname>Thomson</surname>
<given-names>RC</given-names>
</name>
, <name><surname>Shaffer</surname>
<given-names>HB</given-names>
</name>
(<year>2010</year>
) <article-title>Rapid progress on the vertebrate tree of life</article-title>
. <source>BMC Biol</source>
<volume>8</volume>
: <fpage>19</fpage>
<pub-id pub-id-type="pmid">20211001</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Gaffney1"><label>6</label>
<mixed-citation publication-type="other">Gaffney ES (1980) Phylogenetic relationships of the major groups of amniotes. In: Panchen AL, editor. The terrestrial environment and the origin of land vertebrates. New York: Academic Press. 593–610.</mixed-citation>
</ref>
<ref id="pone.0048990-Lyson1"><label>7</label>
<mixed-citation publication-type="journal"><name><surname>Lyson</surname>
<given-names>TR</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>MicroRNAs support a turtle + lizard clade</article-title>
. <source>Biol Lett</source>
<volume>8</volume>
: <fpage>104</fpage>
–<lpage>107</lpage>
<pub-id pub-id-type="pmid">21775315</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Crawford1"><label>8</label>
<mixed-citation publication-type="other">Crawford NG, Faircloth BC, McCormack JE, Brumfield RT, Winker K, <etal>et al</etal>
. (2012) More than 1000 ultraconserved elements provide evidence that turtles are the sister group of archosaurs. Biol Lett. doi:10.1098/rsbl.2012.0331.</mixed-citation>
</ref>
<ref id="pone.0048990-Chiari1"><label>9</label>
<mixed-citation publication-type="journal"><name><surname>Chiari</surname>
<given-names>Y</given-names>
</name>
, <name><surname>Cahais</surname>
<given-names>V</given-names>
</name>
, <name><surname>Galtier</surname>
<given-names>N</given-names>
</name>
, <name><surname>Delsuc</surname>
<given-names>F</given-names>
</name>
(<year>2012</year>
) <article-title>Phylogenomic analyses support the position of turtles as the sister group of birds and crocodiles (Archosauria)</article-title>
. <source>BMC Biol</source>
<volume>10</volume>
: <fpage>65</fpage>
<pub-id pub-id-type="pmid">22839781</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Parsons1"><label>10</label>
<mixed-citation publication-type="journal"><name><surname>Parsons</surname>
<given-names>TS</given-names>
</name>
, <name><surname>Williams</surname>
<given-names>EE</given-names>
</name>
(<year>1963</year>
) <article-title>The relationship of the modern Amphibia: a re-examination</article-title>
. <source>Q Rev Biol</source>
<volume>38</volume>
: <fpage>26</fpage>
–<lpage>53</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Carroll1"><label>11</label>
<mixed-citation publication-type="other">Carroll RL (2009) The Rise of Amphibians. Baltimore: The Johns Hopkins University Press. 360p.</mixed-citation>
</ref>
<ref id="pone.0048990-Carroll2"><label>12</label>
<mixed-citation publication-type="journal"><name><surname>Carroll</surname>
<given-names>RL</given-names>
</name>
(<year>2001</year>
) <article-title>The origin and early radiation of terrestrial vertebrates</article-title>
. <source>J Palaeontol</source>
<volume>75</volume>
: <fpage>1202</fpage>
–<lpage>1213</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Cannatella1"><label>13</label>
<mixed-citation publication-type="other">Cannatella DC, Vieites DR, Zhang P, Wake MH, Wake DB (2009) Amphibians (Lissamphibia). In: Hedges SB, Kumar S, editors. The timetree of life. New York: Oxford University Press. 353–356.</mixed-citation>
</ref>
<ref id="pone.0048990-Vallin1"><label>14</label>
<mixed-citation publication-type="journal"><name><surname>Vallin</surname>
<given-names>G</given-names>
</name>
, <name><surname>Laurin</surname>
<given-names>M</given-names>
</name>
(<year>2004</year>
) <article-title>Cranial morphology and affinities of Microbrachis, and a reappraisal of the phylogeny and lifestyle of the first amphibians</article-title>
. <source>J Vert Paleontol</source>
<volume>24</volume>
: <fpage>56</fpage>
–<lpage>72</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Hedges1"><label>15</label>
<mixed-citation publication-type="journal"><name><surname>Hedges</surname>
<given-names>SB</given-names>
</name>
, <name><surname>Maxson</surname>
<given-names>LR</given-names>
</name>
(<year>1993</year>
) <article-title>A molecular perspective on lissamphibian phylogeny</article-title>
. <source>Herpetol Monogr</source>
<volume>7</volume>
: <fpage>27</fpage>
–<lpage>42</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Feller1"><label>16</label>
<mixed-citation publication-type="journal"><name><surname>Feller</surname>
<given-names>AE</given-names>
</name>
, <name><surname>Hedges</surname>
<given-names>SB</given-names>
</name>
(<year>1998</year>
) <article-title>Molecular evidence for the early history of living amphibians</article-title>
. <source>Mol Phylogenet Evol</source>
<volume>9</volume>
: <fpage>509</fpage>
–<lpage>516</lpage>
<pub-id pub-id-type="pmid">9667999</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Milner1"><label>17</label>
<mixed-citation publication-type="journal"><name><surname>Milner</surname>
<given-names>AR</given-names>
</name>
(<year>1993</year>
) <article-title>The Paleozoic relatives of lissamphibians</article-title>
. <source>Herpetol Monogr</source>
<volume>6</volume>
: <fpage>8</fpage>
–<lpage>27</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Trueb1"><label>18</label>
<mixed-citation publication-type="other">Trueb L, Cloutier R (1991) Toward an understanding of the amphibians: two centuries of systematic history. In: Schultze H-P, Trueb L, editors. Origins of the higher groups of tetrapods: controversy and consensus. Ithaca: Cornell University Press. 233–313.</mixed-citation>
</ref>
<ref id="pone.0048990-Ruta1"><label>19</label>
<mixed-citation publication-type="journal"><name><surname>Ruta</surname>
<given-names>M</given-names>
</name>
, <name><surname>Coates</surname>
<given-names>MI</given-names>
</name>
, <name><surname>Quicke</surname>
<given-names>DLJ</given-names>
</name>
(<year>2003</year>
) <article-title>Early tetrapod relationships revisited</article-title>
. <source>Biol Rev</source>
<volume>78</volume>
: <fpage>251</fpage>
–<lpage>345</lpage>
<pub-id pub-id-type="pmid">12803423</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Zardoya1"><label>20</label>
<mixed-citation publication-type="journal"><name><surname>Zardoya</surname>
<given-names>R</given-names>
</name>
, <name><surname>Meyer</surname>
<given-names>A</given-names>
</name>
(<year>2001</year>
) <article-title>On the origin and phylogenetic relationships among living amphibians</article-title>
. <source>Proc Natl Acad Sci USA</source>
<volume>98</volume>
: <fpage>7380</fpage>
–<lpage>7383</lpage>
<pub-id pub-id-type="pmid">11390961</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Hugall1"><label>21</label>
<mixed-citation publication-type="journal"><name><surname>Hugall</surname>
<given-names>AF</given-names>
</name>
, <name><surname>Foster</surname>
<given-names>R</given-names>
</name>
, <name><surname>Lee</surname>
<given-names>MSY</given-names>
</name>
(<year>2007</year>
) <article-title>Calibration choice, rate smoothing, and the pattern of tetrapod diversification according to the long nuclear gene RAG-1</article-title>
. <source>Syst Biol</source>
<volume>56</volume>
: <fpage>543</fpage>
–<lpage>563</lpage>
<pub-id pub-id-type="pmid">17654361</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Roleants1"><label>22</label>
<mixed-citation publication-type="journal"><name><surname>Roleants</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Global patterns of diversification in the history of modern amphibians</article-title>
. <source>Proc Natl Acad Sci USA</source>
<volume>104</volume>
: <fpage>887</fpage>
–<lpage>892</lpage>
<pub-id pub-id-type="pmid">17213318</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-SanMauro1"><label>23</label>
<mixed-citation publication-type="journal"><name><surname>San Mauro</surname>
<given-names>D</given-names>
</name>
, <name><surname>Vences</surname>
<given-names>M</given-names>
</name>
, <name><surname>Alcobendas</surname>
<given-names>M</given-names>
</name>
, <name><surname>Zardoya</surname>
<given-names>R</given-names>
</name>
, <name><surname>Meyer</surname>
<given-names>A</given-names>
</name>
(<year>2005</year>
) <article-title>Initial diversification of living amphibians predated the breakup of Pangaea</article-title>
. <source>Am Nat</source>
<volume>165</volume>
: <fpage>590</fpage>
–<lpage>599</lpage>
<pub-id pub-id-type="pmid">15795855</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-SanMauro2"><label>24</label>
<mixed-citation publication-type="journal"><name><surname>San Mauro</surname>
<given-names>D</given-names>
</name>
(<year>2010</year>
) <article-title>A multilocus timescale for the origin of extant amphibians</article-title>
. <source>Mol Phylogenet Evol</source>
<volume>56</volume>
: <fpage>554</fpage>
–<lpage>561</lpage>
<pub-id pub-id-type="pmid">20399871</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Anderson1"><label>25</label>
<mixed-citation publication-type="journal"><name><surname>Anderson</surname>
<given-names>JS</given-names>
</name>
, <name><surname>Reisz</surname>
<given-names>RR</given-names>
</name>
, <name><surname>Scott</surname>
<given-names>D</given-names>
</name>
, <name><surname>Frobisch</surname>
<given-names>NB</given-names>
</name>
, <name><surname>Sumida</surname>
<given-names>SS</given-names>
</name>
(<year>2008</year>
) <article-title>A stem batrachian from the Early Permian of Texas and the origin of frogs and salamanders</article-title>
. <source>Nature</source>
<volume>453</volume>
: <fpage>515</fpage>
–<lpage>518</lpage>
<pub-id pub-id-type="pmid">18497824</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Anderson2"><label>26</label>
<mixed-citation publication-type="journal"><name><surname>Anderson</surname>
<given-names>JS</given-names>
</name>
(<year>2008</year>
) <article-title>Focal review: the origin(s) of modern amphibians</article-title>
. <source>Evol Biol</source>
<volume>35</volume>
: <fpage>231</fpage>
–<lpage>247</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Carroll3"><label>27</label>
<mixed-citation publication-type="journal"><name><surname>Carroll</surname>
<given-names>RL</given-names>
</name>
(<year>2007</year>
) <article-title>The Palaeozoic ancestry of salamanders, frogs, and caecilians</article-title>
. <source>Zool J Linn Soc</source>
<volume>150</volume>
(suppl. 1)<fpage>1</fpage>
–<lpage>142</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Carroll4"><label>28</label>
<mixed-citation publication-type="journal"><name><surname>Carroll</surname>
<given-names>RL</given-names>
</name>
, <name><surname>Holmes</surname>
<given-names>R</given-names>
</name>
(<year>1980</year>
) <article-title>The skull and jaw musculature as guides to the ancestry of salamanders</article-title>
. <source>Zool J Linn Soc</source>
<volume>68</volume>
: <fpage>1</fpage>
–<lpage>40</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Carroll5"><label>29</label>
<mixed-citation publication-type="other">Carroll RL, Boisvert C, Bolt J, Green DM, Philip N, <etal>et al</etal>
. (2004) Changing patterns of ontogeny from osteopiform fish through Permian tetrapods as a guide to the early evolution of land vertebrates. In: Arratia G, Wilson MHV, Coutier R, editors. Recent advances in the origin and early radiation of vertebrates. Munchen: Verlag Dr. Friedrich Pfeil, Munchen. 321–343.</mixed-citation>
</ref>
<ref id="pone.0048990-Li1"><label>30</label>
<mixed-citation publication-type="journal"><name><surname>Li</surname>
<given-names>C</given-names>
</name>
, <name><surname>Wu</surname>
<given-names>X</given-names>
</name>
, <name><surname>Rieppel</surname>
<given-names>O</given-names>
</name>
, <name><surname>Wang</surname>
<given-names>L</given-names>
</name>
, <name><surname>Zhao</surname>
<given-names>L</given-names>
</name>
(<year>2008</year>
) <article-title>An ancestral turtle from the Late Triassic of southwestern China</article-title>
. <source>Nature</source>
<volume>456</volume>
: <fpage>497</fpage>
–<lpage>501</lpage>
<pub-id pub-id-type="pmid">19037315</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Philippe1"><label>31</label>
<mixed-citation publication-type="journal"><name><surname>Philippe</surname>
<given-names>H</given-names>
</name>
, <name><surname>Brinkmann</surname>
<given-names>H</given-names>
</name>
, <name><surname>Lavrov</surname>
<given-names>DV</given-names>
</name>
, <name><surname>Littlewood</surname>
<given-names>TJ</given-names>
</name>
, <name><surname>Manuel</surname>
<given-names>M</given-names>
</name>
, <etal>et al</etal>
(<year>2011</year>
) <article-title>Solving difficult phylogenetic questions: why more sequences are not enough</article-title>
. <source>PLoS Biol</source>
<volume>9</volume>
: <fpage>e1000602</fpage>
<pub-id pub-id-type="pmid">21423652</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Jeffroy1"><label>32</label>
<mixed-citation publication-type="journal"><name><surname>Jeffroy</surname>
<given-names>O</given-names>
</name>
, <name><surname>Brinkmann</surname>
<given-names>H</given-names>
</name>
, <name><surname>Delsuc</surname>
<given-names>F</given-names>
</name>
, <name><surname>Philippe</surname>
<given-names>H</given-names>
</name>
(<year>2006</year>
) <article-title>Phylogenomics: the beginning of incongruence?</article-title>
<source>Trends Genet</source>
<volume>22</volume>
: <fpage>225</fpage>
–<lpage>231</lpage>
<pub-id pub-id-type="pmid">16490279</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Pattengale1"><label>33</label>
<mixed-citation publication-type="journal"><name><surname>Pattengale</surname>
<given-names>ND</given-names>
</name>
, <name><surname>Swenson</surname>
<given-names>KM</given-names>
</name>
, <name><surname>Moret</surname>
<given-names>BME</given-names>
</name>
(<year>2010</year>
) <article-title>Uncovering hidden phylogenetic consensus</article-title>
. <source>Bioinformat Res Appl</source>
<volume>6053</volume>
: <fpage>128</fpage>
–<lpage>139</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Meyer2"><label>34</label>
<mixed-citation publication-type="journal"><name><surname>Meyer</surname>
<given-names>A</given-names>
</name>
(<year>1994</year>
) <article-title>Shortcomings of the cytochrome b gene as a molecular marker</article-title>
. <source>T Ecol Evol</source>
<volume>9</volume>
: <fpage>278</fpage>
–<lpage>280</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Simon1"><label>35</label>
<mixed-citation publication-type="journal"><name><surname>Simon</surname>
<given-names>C</given-names>
</name>
, <name><surname>Buckley</surname>
<given-names>TR</given-names>
</name>
, <name><surname>Frati</surname>
<given-names>F</given-names>
</name>
, <name><surname>Stewart</surname>
<given-names>JB</given-names>
</name>
, <name><surname>Beckenbach</surname>
<given-names>AT</given-names>
</name>
(<year>1994</year>
) <article-title>Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers</article-title>
. <source>Ann Entomol Soc Am</source>
<volume>87</volume>
: <fpage>651</fpage>
–<lpage>701</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Edwards1"><label>36</label>
<mixed-citation publication-type="journal"><name><surname>Edwards</surname>
<given-names>SV</given-names>
</name>
, <name><surname>Arctander</surname>
<given-names>P</given-names>
</name>
, <name><surname>Wilson</surname>
<given-names>AC</given-names>
</name>
(<year>1991</year>
) <article-title>Mitochondrial resolution of a deep branch in the genealogical tree for perching birds</article-title>
. <source>Proc R Soc Lond B</source>
<volume>243</volume>
: <fpage>99</fpage>
–<lpage>107</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Blouin1"><label>37</label>
<mixed-citation publication-type="journal"><name><surname>Blouin</surname>
<given-names>MS</given-names>
</name>
, <name><surname>Yowell</surname>
<given-names>CA</given-names>
</name>
, <name><surname>Courtnely</surname>
<given-names>CH</given-names>
</name>
, <name><surname>Dame</surname>
<given-names>JB</given-names>
</name>
(<year>1998</year>
) <article-title>Substitution bias, rapid saturation, and the use of mtDNA for nematode systematics</article-title>
. <source>Mol Biol Evol</source>
<volume>15</volume>
: <fpage>1719</fpage>
–<lpage>1727</lpage>
<pub-id pub-id-type="pmid">9866206</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Regier1"><label>38</label>
<mixed-citation publication-type="journal"><name><surname>Regier</surname>
<given-names>JC</given-names>
</name>
, <name><surname>Schultz</surname>
<given-names>JW</given-names>
</name>
, <name><surname>Zwick</surname>
<given-names>A</given-names>
</name>
, <name><surname>Hussey</surname>
<given-names>A</given-names>
</name>
, <name><surname>Ball</surname>
<given-names>B</given-names>
</name>
, <etal>et al</etal>
(<year>2010</year>
) <article-title>Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences</article-title>
. <source>Nature</source>
<volume>463</volume>
: <fpage>1079</fpage>
–<lpage>1083</lpage>
<pub-id pub-id-type="pmid">20147900</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Brinkmann1"><label>39</label>
<mixed-citation publication-type="journal"><name><surname>Brinkmann</surname>
<given-names>H</given-names>
</name>
, <name><surname>Philippe</surname>
<given-names>H</given-names>
</name>
(<year>1999</year>
) <article-title>Archaea sister group of bacteria? Indications from tree reconstruction artifacts in ancient phylogenies</article-title>
. <source>Mol Biol Evol</source>
<volume>16</volume>
: <fpage>817</fpage>
–<lpage>825</lpage>
<pub-id pub-id-type="pmid">10368959</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Philippe2"><label>40</label>
<mixed-citation publication-type="journal"><name><surname>Philippe</surname>
<given-names>H</given-names>
</name>
, <name><surname>Lartillot</surname>
<given-names>N</given-names>
</name>
, <name><surname>Brinkmann</surname>
<given-names>H</given-names>
</name>
(<year>2005</year>
) <article-title>Multigene analyses of bilaterian animals corroborate the monophyly of Ecdysozoa, Lophotrochozoa, and Protostomia</article-title>
. <source>Mol Biol Evol</source>
<volume>22</volume>
: <fpage>1246</fpage>
–<lpage>1253</lpage>
<pub-id pub-id-type="pmid">15703236</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Fong1"><label>41</label>
<mixed-citation publication-type="journal"><name><surname>Fong</surname>
<given-names>JJ</given-names>
</name>
, <name><surname>Fujita</surname>
<given-names>MK</given-names>
</name>
(<year>2011</year>
) <article-title>Evaluating phylogenetic informativeness and data-type usage for new protein-coding genes across Vertebrata</article-title>
. <source>Mol Phylogenet Evol</source>
<volume>61</volume>
: <fpage>300</fpage>
–<lpage>307</lpage>
<pub-id pub-id-type="pmid">21742044</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Edwards2"><label>42</label>
<mixed-citation publication-type="journal"><name><surname>Edwards</surname>
<given-names>SV</given-names>
</name>
, <name><surname>Liu</surname>
<given-names>L</given-names>
</name>
, <name><surname>Pearl</surname>
<given-names>DK</given-names>
</name>
(<year>2007</year>
) <article-title>High-resolution species trees without concatenation</article-title>
. <source>Proc Natl Acad Sci USA</source>
<volume>104</volume>
: <fpage>5936</fpage>
–<lpage>5941</lpage>
<pub-id pub-id-type="pmid">17392434</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Cranston1"><label>43</label>
<mixed-citation publication-type="journal"><name><surname>Cranston</surname>
<given-names>KA</given-names>
</name>
, <name><surname>Hurwitz</surname>
<given-names>B</given-names>
</name>
, <name><surname>Ware</surname>
<given-names>D</given-names>
</name>
, <name><surname>Stein</surname>
<given-names>L</given-names>
</name>
, <name><surname>Wing</surname>
<given-names>RA</given-names>
</name>
(<year>2009</year>
) <article-title>Species trees from highly incongruent gene trees in rice</article-title>
. <source>Syst Biol</source>
<volume>58</volume>
: <fpage>489</fpage>
–<lpage>500</lpage>
<pub-id pub-id-type="pmid">20525603</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Shimodaira1"><label>44</label>
<mixed-citation publication-type="journal"><name><surname>Shimodaira</surname>
<given-names>H</given-names>
</name>
(<year>2002</year>
) <article-title>An approximately unbiased test of phylogenetic tree selection</article-title>
. <source>Syst Biol</source>
<volume>51</volume>
: <fpage>492</fpage>
–<lpage>508</lpage>
<pub-id pub-id-type="pmid">12079646</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-RodriguezEzpeleta1"><label>45</label>
<mixed-citation publication-type="journal"><name><surname>Rodriguez-Ezpeleta</surname>
<given-names>N</given-names>
</name>
, <name><surname>Brinkmann</surname>
<given-names>H</given-names>
</name>
, <name><surname>Roure</surname>
<given-names>B</given-names>
</name>
, <name><surname>Lartillot</surname>
<given-names>N</given-names>
</name>
, <name><surname>Lang</surname>
<given-names>BF</given-names>
</name>
, <etal>et al</etal>
(<year>2007</year>
) <article-title>Detecting and overcoming systematic errors in genome-scale phylogenies</article-title>
. <source>Syst Biol</source>
<volume>56</volume>
: <fpage>389</fpage>
–<lpage>399</lpage>
<pub-id pub-id-type="pmid">17520503</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Lemmon1"><label>46</label>
<mixed-citation publication-type="journal"><name><surname>Lemmon</surname>
<given-names>AR</given-names>
</name>
, <name><surname>Brown</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Stanger-Hall</surname>
<given-names>K</given-names>
</name>
, <name><surname>Lemmon</surname>
<given-names>EM</given-names>
</name>
(<year>2009</year>
) <article-title>The effect of ambiguous data on phylogenetic estimates obtained by maximum likelihood and Bayesian inference</article-title>
. <source>Syst Biol</source>
<volume>58</volume>
: <fpage>130</fpage>
–<lpage>145</lpage>
<pub-id pub-id-type="pmid">20525573</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Wolf1"><label>47</label>
<mixed-citation publication-type="journal"><name><surname>Wolf</surname>
<given-names>YI</given-names>
</name>
, <name><surname>Rogozin</surname>
<given-names>IB</given-names>
</name>
, <name><surname>Koonin</surname>
<given-names>EV</given-names>
</name>
(<year>2003</year>
) <article-title>Coelomata and not Ecdysozoa: evidence from genome-wide phylogenetic analysis</article-title>
. <source>Genome Res</source>
<volume>14</volume>
: <fpage>29</fpage>
–<lpage>36</lpage>
<pub-id pub-id-type="pmid">14707168</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Aguinaldo1"><label>48</label>
<mixed-citation publication-type="journal"><name><surname>Aguinaldo</surname>
<given-names>AM</given-names>
</name>
, <name><surname>Tubeville</surname>
<given-names>JM</given-names>
</name>
, <name><surname>Linford</surname>
<given-names>LS</given-names>
</name>
, <name><surname>Rivera</surname>
<given-names>MC</given-names>
</name>
, <name><surname>Raff</surname>
<given-names>RA</given-names>
</name>
, <etal>et al</etal>
(<year>1997</year>
) <article-title>Evidence for a clade of nematodes, arthropods and other moulting animals</article-title>
. <source>Nature</source>
<volume>387</volume>
: <fpage>489</fpage>
–<lpage>493</lpage>
<pub-id pub-id-type="pmid">9168109</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Telford1"><label>49</label>
<mixed-citation publication-type="journal"><name><surname>Telford</surname>
<given-names>MJ</given-names>
</name>
, <name><surname>Bourlat</surname>
<given-names>SJ</given-names>
</name>
, <name><surname>Economou</surname>
<given-names>A</given-names>
</name>
, <name><surname>Papillon</surname>
<given-names>D</given-names>
</name>
, <name><surname>Rota-Stabelli</surname>
<given-names>O</given-names>
</name>
(<year>2008</year>
) <article-title>The evolution of the Ecdysozoa</article-title>
. <source>Phil Trans R Soc B</source>
<volume>363</volume>
: <fpage>1529</fpage>
–<lpage>1537</lpage>
<pub-id pub-id-type="pmid">18192181</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Gauthier1"><label>50</label>
<mixed-citation publication-type="other">Gauthier J, Padian K (1985) Phylogenetic, functional, and aerodynamic analyses of the origin of birds and their flight. In: Hecht JH, Ostrom GV, Wellnhofer P, editors. The beginnings of birds. Eichstatt: Freunde des Jura-Museum. 185–197.</mixed-citation>
</ref>
<ref id="pone.0048990-Bhullar1"><label>51</label>
<mixed-citation publication-type="journal"><name><surname>Bhullar</surname>
<given-names>B-A</given-names>
</name>
, <name><surname>Bever</surname>
<given-names>GS</given-names>
</name>
(<year>2009</year>
) <article-title>An archosaur-like laterosphenoid in early turtles (Reptilia: Pantestudines)</article-title>
. <source>Breviora</source>
<volume>518</volume>
: <fpage>1</fpage>
–<lpage>11</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Laurin1"><label>52</label>
<mixed-citation publication-type="journal"><name><surname>Laurin</surname>
<given-names>M</given-names>
</name>
, <name><surname>Reisz</surname>
<given-names>R</given-names>
</name>
(<year>1995</year>
) <article-title>A reevaluation of early amniote phylogeny</article-title>
. <source>Zool J Linn Soc</source>
<volume>113</volume>
: <fpage>165</fpage>
–<lpage>223</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Lee1"><label>53</label>
<mixed-citation publication-type="journal"><name><surname>Lee</surname>
<given-names>MSY</given-names>
</name>
(<year>1995</year>
) <article-title>Historical burden in systematics and the interrelationships of ‘parareptiles’</article-title>
. <source>Biol Rev</source>
<volume>70</volume>
: <fpage>459</fpage>
–<lpage>547</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Lee2"><label>54</label>
<mixed-citation publication-type="journal"><name><surname>Lee</surname>
<given-names>MSY</given-names>
</name>
(<year>1996</year>
) <article-title>Correlated progression and the origin of turtles</article-title>
. <source>Nature</source>
<volume>379</volume>
: <fpage>812</fpage>
–<lpage>815</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Nagashima1"><label>55</label>
<mixed-citation publication-type="journal"><name><surname>Nagashima</surname>
<given-names>H</given-names>
</name>
, <name><surname>Sugahara</surname>
<given-names>F</given-names>
</name>
, <name><surname>Takechi</surname>
<given-names>M</given-names>
</name>
, <name><surname>Ericsson</surname>
<given-names>R</given-names>
</name>
, <name><surname>Kawashima-Ohya</surname>
<given-names>Y</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Evolution of the turtle body plan by the folding and creation of new muscle connections</article-title>
. <source>Science</source>
<volume>325</volume>
: <fpage>193</fpage>
–<lpage>196</lpage>
<pub-id pub-id-type="pmid">19590000</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Yates2"><label>56</label>
<mixed-citation publication-type="journal"><name><surname>Yates</surname>
<given-names>AM</given-names>
</name>
, <name><surname>Warren</surname>
<given-names>A</given-names>
</name>
(<year>2000</year>
) <article-title>The phylogeny of ‘higher’ temnospondyls (Vertebrata: Choanata) and its implications for the monophyly and origins of the Stereospondyli</article-title>
. <source>Zool J Linn Soc</source>
<volume>128</volume>
: <fpage>77</fpage>
–<lpage>121</lpage>
</mixed-citation>
</ref>
<ref id="pone.0048990-Hedges2"><label>57</label>
<mixed-citation publication-type="journal"><name><surname>Hedges</surname>
<given-names>SB</given-names>
</name>
, <name><surname>Poling</surname>
<given-names>LL</given-names>
</name>
(<year>1999</year>
) <article-title>A molecular phylogeny of reptiles</article-title>
. <source>Science</source>
<volume>283</volume>
: <fpage>998</fpage>
–<lpage>1001</lpage>
<pub-id pub-id-type="pmid">9974396</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Hubbard1"><label>58</label>
<mixed-citation publication-type="journal"><name><surname>Hubbard</surname>
<given-names>TJP</given-names>
</name>
, <name><surname>Aken</surname>
<given-names>BL</given-names>
</name>
, <name><surname>Ayling</surname>
<given-names>S</given-names>
</name>
, <name><surname>Ballester</surname>
<given-names>B</given-names>
</name>
, <name><surname>Beal</surname>
<given-names>K</given-names>
</name>
, <etal>et al</etal>
(<year>2009</year>
) <article-title>Ensembl 2009</article-title>
. <source>Nucleic Acids Res</source>
<volume>37</volume>
: <fpage>D690</fpage>
–<lpage>D697</lpage>
<pub-id pub-id-type="pmid">19033362</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Stamatakis1"><label>59</label>
<mixed-citation publication-type="journal"><name><surname>Stamatakis</surname>
<given-names>A</given-names>
</name>
(<year>2006</year>
) <article-title>RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models</article-title>
. <source>Bioinformatics</source>
<volume>22</volume>
: <fpage>2688</fpage>
–<lpage>2690</lpage>
<pub-id pub-id-type="pmid">16928733</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Edgar1"><label>60</label>
<mixed-citation publication-type="journal"><name><surname>Edgar</surname>
<given-names>RC</given-names>
</name>
(<year>2004</year>
) <article-title>MUSCLE: multiple sequence alignment with high accuracy and high throughput</article-title>
. <source>Nucleic Acids Res</source>
<volume>32</volume>
: <fpage>1792</fpage>
–<lpage>1797</lpage>
<pub-id pub-id-type="pmid">15034147</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Maddison1"><label>61</label>
<mixed-citation publication-type="other">Maddison DR, Maddison WP (2005) MacClade 4: Analysis of phylogeny and character evolution. Version 4.08a. <ext-link ext-link-type="uri" xlink:href="http://macclade.org">http://macclade.org</ext-link>
</mixed-citation>
</ref>
<ref id="pone.0048990-R1"><label>62</label>
<mixed-citation publication-type="other">R Development Core Team (2011). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0. URL: <ext-link ext-link-type="uri" xlink:href="http://www.R-project.org/">http://www.R-project.org/</ext-link>
</mixed-citation>
</ref>
<ref id="pone.0048990-Ronquist1"><label>63</label>
<mixed-citation publication-type="journal"><name><surname>Ronquist</surname>
<given-names>F</given-names>
</name>
, <name><surname>Huelsenbeck</surname>
<given-names>JP</given-names>
</name>
(<year>2003</year>
) <article-title>MrBayes 3: Bayesian phylogenetic inference under mixed models</article-title>
. <source>Bioinformatics</source>
<volume>19</volume>
: <fpage>1572</fpage>
–<lpage>1574</lpage>
<pub-id pub-id-type="pmid">12912839</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Abascal1"><label>64</label>
<mixed-citation publication-type="journal"><name><surname>Abascal</surname>
<given-names>F</given-names>
</name>
, <name><surname>Zardoya</surname>
<given-names>R</given-names>
</name>
, <name><surname>Posada</surname>
<given-names>D</given-names>
</name>
(<year>2005</year>
) <article-title>ProtTest: Selection of best-fit models of protein evolution</article-title>
. <source>Bioinformatics</source>
<volume>21</volume>
: <fpage>2104</fpage>
–<lpage>2105</lpage>
<pub-id pub-id-type="pmid">15647292</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Stamatakis2"><label>65</label>
<mixed-citation publication-type="journal"><name><surname>Stamatakis</surname>
<given-names>A</given-names>
</name>
, <name><surname>Hoover</surname>
<given-names>P</given-names>
</name>
, <name><surname>Rougemont</surname>
<given-names>J</given-names>
</name>
(<year>2008</year>
) <article-title>A rapid bootstrap algorithm for the RAxML web-servers</article-title>
. <source>Syst Biol</source>
<volume>75</volume>
: <fpage>758</fpage>
–<lpage>771</lpage>
<pub-id pub-id-type="pmid">18853362</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Miller1"><label>66</label>
<mixed-citation publication-type="other">Miller MA, Holder MT, Vos R, Midford PE, Liebowitz T, <etal>et al</etal>
. (2009) The CIPRES Portals. CIPRES. 2009-08-04. URL: <ext-link ext-link-type="uri" xlink:href="http://www.phylo.org/sub_sections/portal">http://www.phylo.org/sub_sections/portal</ext-link>
Accessed: 2009-08-04. (Archived by WebCite(r) at <ext-link ext-link-type="uri" xlink:href="http://www.webcitation.org/5imQlJeQa">http://www.webcitation.org/5imQlJeQa</ext-link>
).</mixed-citation>
</ref>
<ref id="pone.0048990-Nylander1"><label>67</label>
<mixed-citation publication-type="other">Nylander JAA (2004) MrModeltest 2.3. Program distributed by the author. Evolutionary Biology Centre, Uppsala University.</mixed-citation>
</ref>
<ref id="pone.0048990-Wilgenbusch1"><label>68</label>
<mixed-citation publication-type="other">Wilgenbusch JC, Warren DL, Swofford DL (2004) AWTY: A system for graphical exploration of MCMC convergence in Bayesian phylogenetic inference. <ext-link ext-link-type="uri" xlink:href="http://ceb.csit.fsu.edu/awty">http://ceb.csit.fsu.edu/awty</ext-link>
</mixed-citation>
</ref>
<ref id="pone.0048990-Shimodaira2"><label>69</label>
<mixed-citation publication-type="journal"><name><surname>Shimodaira</surname>
<given-names>H</given-names>
</name>
, <name><surname>Hasegawa</surname>
<given-names>M</given-names>
</name>
(<year>2001</year>
) <article-title>CONSEL: for assessing the confidence of phylogenetic tree selection</article-title>
. <source>Bioinformatics</source>
<volume>17</volume>
: <fpage>1246</fpage>
–<lpage>1247</lpage>
<pub-id pub-id-type="pmid">11751242</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Charif1"><label>70</label>
<mixed-citation publication-type="other">Charif D, Lobry JR (2007) SeqinR 1.0–2: a contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis. In: Bastolla U, Porto M, Roman HE, Vendruscolo M, editors. Structural approaches to sequence evolution: molecules, networks, populations. Springer Verlag, New York. 207–232.</mixed-citation>
</ref>
<ref id="pone.0048990-Venables1"><label>71</label>
<mixed-citation publication-type="other">Venables WN, Ripley BD (2002) Modern applied statistics with S, 4<sup>th</sup>
edition. New York: Springer. 512p.</mixed-citation>
</ref>
<ref id="pone.0048990-Pond1"><label>72</label>
<mixed-citation publication-type="journal"><name><surname>Pond</surname>
<given-names>SLK</given-names>
</name>
, <name><surname>Frost</surname>
<given-names>SDW</given-names>
</name>
, <name><surname>Muse</surname>
<given-names>SV</given-names>
</name>
(<year>2005</year>
) <article-title>HyPhy: hypothesis testing using phylogenies</article-title>
. <source>Bioinformatics</source>
<volume>21</volume>
: <fpage>676</fpage>
–<lpage>679</lpage>
<pub-id pub-id-type="pmid">15509596</pub-id>
</mixed-citation>
</ref>
<ref id="pone.0048990-Dutheil1"><label>73</label>
<mixed-citation publication-type="journal"><name><surname>Dutheil</surname>
<given-names>J</given-names>
</name>
, <name><surname>Boussau</surname>
<given-names>B</given-names>
</name>
(<year>2008</year>
) <article-title>Non-homogeneous models of sequence evolution in the Bio++ suite of libraries and programs. BMC Evol Biol</article-title>
. <volume>8</volume>
: <fpage>255</fpage>
</mixed-citation>
</ref>
</ref-list>
</back>
</pmc>
</record>
Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Ticri/CIDE/explor/CyberinfraV1/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000612 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000612 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien |wiki= Ticri/CIDE |area= CyberinfraV1 |flux= Pmc |étape= Corpus |type= RBID |clé= PMC:3492174 |texte= A Phylogenomic Approach to Vertebrate Phylogeny Supports a Turtle-Archosaur Affinity and a Possible Paraphyletic Lissamphibia }}
Pour générer des pages wiki
HfdIndexSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/RBID.i -Sk "pubmed:23145043" \ | HfdSelect -Kh $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd \ | NlmPubMed2Wicri -a CyberinfraV1
![]() | This area was generated with Dilib version V0.6.25. | ![]() |