Marginal likelihoods in phylogenetics: a review of methods and applications

By providing a framework of accounting for the shared ancestry inherent to all life, phylogenetics is becoming the statistical foundation of biology. The importance of model choice continues to grow as phylogenetic models continue to increase in complexity to better capture micro and macroevolutionary processes. In a Bayesian framework, the marginal likelihood is how data update our prior beliefs about models, which gives us an intuitive measure of comparing model fit that is grounded in probability theory. Given the rapid increase in the number and complexity of phylogenetic models, methods for approximating marginal likelihoods are increasingly important. Here we try to provide an intuitive description of marginal likelihoods and why they are important in Bayesian model testing. We also categorize and review methods for estimating marginal likelihoods of phylogenetic models, highlighting several recent methods that provide well-behaved estimates. Furthermore, we review some empirical studies that demonstrate how marginal likelihoods can be used to learn about models of evolution from biological data. We discuss promising alternatives that can complement marginal likelihoods for Bayesian model choice, including posterior-predictive methods. Using simulations, we find one alternative method based on approximate-Bayesian computation (ABC) to be biased. We conclude by discussing the challenges of Bayesian model choice and future directions that promise to improve the approximation of marginal likelihoods and Bayesian phylogenetics as a whole.Comment: 33 pages, 3 figure

Coming to America: Multiple Origins of New World Geckos

Geckos in the Western Hemisphere provide an excellent model to study faunal assembly at a continental scale. We generated a time-calibrated phylogeny, including exemplars of all New World gecko genera, to produce a biogeographic scenario for the New World geckos. Patterns of New World gecko origins are consistent with almost every biogeographic scenario utilized by a terrestrial vertebrate with different New World lineages showing evidence of vicariance, dispersal via temporary land bridge, overseas dispersal, or anthropogenic introductions. We also recovered a strong relationship between clade age and species diversity, with older New World lineages having more species than more recently arrived lineages. Our data provide the first phylogenetic hypothesis for all New World geckos and highlight the intricate origins and ongoing organization of continental faunas. The phylogenetic and biogeographical hypotheses presented here provide an historical framework to further pursue research on the diversification and assembly of the New World herpetofauna

Molecular phylogenetics of Melastomataceae and Memecylaceae

Melastomataceae are among the most abundant and diversified groups of plants throughout the tropics, but their intrafamily relationships and morphological evolution are poorly understood. Here we report the results of parsimony and maximum likelihood (ML) analyses of cpDNA sequences from the rbcL and ndhF genes and the rpl16 intron, generated for eight outgroups (Crypteroniaceae, Alzateaceae, Rhynchocalycaceae, Oliniaceae, Penaeaceae, Myrtaceae, and Onagraceae) and 54 species of melastomes. The sample represents 42 of the family’s currently recognized ~150 genera, the 13 traditional tribes, and the three subfamilies, Astronioideae, Melastomatoideae, and Memecyloideae (= Memecylaceae DC.). Parsimony and ML yield congruent topologies that place Memecylaceae as sister to Melastomataceae. Pternandra, a Southeast Asian genus of 15 species of which five were sampled, is the firstbranching Melastomataceae. This placement has low bootstrap support (72%), but agrees with morphological treatments that placed Pternandra in Melastomatacaeae because of its acrodromal leaf venation, usually ranked as a tribe or subfamily. The interxylary phloem islands found in Memecylaceae and Pternandra, but not most other Melastomataceae, likely evolved in parallel because Pternandra resembles Melastomataceae in its other wood characters. A newly discovered plesiomorphic character in Pternandra, also present in Memecylaceae, is a fibrous anther endothecium. Higher Melastomataceae lack an endothecium as do the closest relatives of Melastomataceae and Memecylaceae. The next deepest split is between Astronieae, with anthers opening by slits, and all remaining Melastomataceae, which have anthers opening by pores. Within the latter, several generic groups, corresponding to traditional tribes, receive solid statistical support, but relationships among them, with one exception, are different from anything predicted on the basis of morphological data. Thus, Miconieae and Merianieae are sister groups, and both are sister to a trichotomy of Bertolonieae, Microlicieae + Melastomeae, and Dissochaeteae + Blakeeae. Sonerileae/Oxysporeae are nested within Dissochaeteae, Rhexieae within Melastomeae, and African and Asian Melastomeae within neotropical Melastomeae. These findings have profound implications for our understanding of melastome morphological evolution (and biogeography), implying, for example, that berries evolved from capsules minimally four times, stamen connectives went from dorsally enlarged to basal/ventrally enlarged, and loss of an endothecium preceded poricidal dehiscence

Exploring deep phylogenies using protein structure : a dissertation submitted in partial fulfilment of the requirements for the degree of Doctor of Philosophy in Biochemistry, Institute of Natural and Mathematical Sciences, Massey University, Auckland, New Zealand

Recent times have seen an exponential growth in protein sequence and structure data. The most popular way of characterising newly determined protein sequences is to compare them to well characterised sequences and predict the function of novel sequences based on homology. This practice has been highly successful for a majority of proteins. However, these sequence based methods struggle with certain deeply diverging proteins and hence cannot always recover evolutionary histories. Another feature of proteins, namely their structures, has been shown to retain evolutionary signals over longer time scales compared to the respective sequences that encode them. The structure therefore presents an opportunity to uncover the evolutionary signal that otherwise escapes conventional sequence-based methods. Structural phylogenetics refers to the comparison of protein structures to extract evolutionary relationships. The area of structural phylogenetics has been around for a number of years and multiple approaches exist to delineate evolutionary relationships from protein structures. However, once the relationships have been recovered from protein structural data, no methods exist, at present, to verify the robustness of these relationships. Because of the nature of the structural data, conventional sequence-based methods, e.g. bootstrapping, cannot be applied. This work introduces the first ever use of a molecular dynamics (MD)-based bootstrap method, which can add a measure of significance to the relationships inferred from the structure-based analysis. This work begins in Chapter 2 by thoroughly investigating the use of a protein structural comparison metric Qscore, which has previously been used to generate structural phylogenies, and highlights its strengths and weaknesses. The mechanistic exploration of the structural comparison metric reveals a size difference limit of no more than 5-10% in the sizes of protein structures being compared for accurate phylogenetic inference to be made. Chapter 2 also explores the MD-based bootstrap method to offer an interpretation of the significance values recovered. Two protein structural datasets, one relatively more conserved at the sequence level than the other and with different levels of structural conservation are used as controls to simplify the interpretation of the statistics recovered from the MD-based bootstrap method. Chapter 3 then sees the application of the Qscore metric to the aminoacyl-tRNA synthetases. The aminoacyl-tRNA synthetases are believed to have been present at the dawn of life, making them one of the most ancient protein families. Due to the important functional role they play, these proteins are conserved at both sequence and structural levels and well-characterised using both sequence and structure-based comparative methods. This family therefore offered inferences which could be informed with structural analysis using an automated method. Successful recovery of known relationships raised confidence in the ability of structural phylogenetic analysis based on Qscore to detect evolutionary signals. In Chapter 4, a structural phylogeny was created for a protein structural dataset presenting either the histone fold or its ancestral precursor. This structural dataset comprised of proteins that were significantly diverged at a sequence level, however shared a common structural motif. The structural phylogeny recovered the split between bacterial and non-bacterial proteins. Furthermore, TATA protein associated factors were found to have multiple points of origin. Moreover, some mismatch was found between the classifications of these proteins between SCOP and PFam, which also did not agree with the results from this work. Using the structural phylogeny a model outlining the evolution of these proteins was proposed. The structural phylogeny of the Ferritin-like superfamily has previously been generated using the Qscore metric and supported qualitatively. Chapter 5 recovers the structural phylogeny of the Ferritin-like superfamily and finds quantitative support for the inferred relationships from the first ever implementation of the MD-based bootstrap method. The use of the MD-based bootstrap method simultaneously allows for the resolution of polytomies in structural databases. Some limitations of the MD-based bootstrap method, highlighted in Chapter 2, are revisited in Chapter 5. This work indicates that evolutionary signals can be successfully extracted from protein structures for deeply diverging proteins and that the MD-based bootstrap method can be used to gauge the robustness of relationships inferred

Genetic Diversity, Phylogenetics and Molecular Systematics of Guizotia Cass. (Asteraceae)

The genus Guizotia belongs to the tribe Heliantheae in the family Asteraceae. It has been placed under different subtribes. The genus has its center of origin, distribution and genetic diversity in Ethiopia, where G. abyssinica (niger) has been domesticated. Amplified Fragment Length Polymorphism (AFLP), Random Amplified Polymorphic DNA (RAPD) and DNA sequencing were applied to study the genetic diversity, phylogenetics, and molecular systematics of this genus. A large number of niger populations, representing all regions in Ethiopia where this crop is grown, was investigated using AFLP and RAPD molecular marker techniques. The extent of genetic variation in niger is distributed throughout its growing regions, regardless of the extent and altitude of cultivation. Despite the fact that most of the variation was within populations, significant population differentiation was obtained (AMOVA; P < 0.001) in all guizotias. It is concluded that both G. abyssinica and its wild and/or weedy relatives have wide genetic bases that need to be conserved and utilized for the improvement of G. abyssinica. Further collection of niger germplasm and exploration and conservation of highly localized guizotias are recommended. Most of the diagnostic markers generated from AFLPs and RAPDs in this study were specific to G. arborescens and G. zavattarii. Phylogenetic analyses of the genus Guizotia were undertaken based on molecular sequence data from the internal transcribed spacers (ITS) and five chloroplast DNA regions. The analyses revealed a close phylogenetic relationship between G. abyssinica and G. scabra ssp. schimperi and support the previous suggestion that the latter is the progenitor of the former. According to this study, G. scabra ssp. scabra, G. scabra ssp. schimperi, and the Chelelu and Ketcha populations are best viewed at present as separate species within the genus Guizotia. Those perennial guizotias with highly localized geographic distribution appears to have evolved first during the evolutionary history of the genus. This study supports the placement of the genus Guizotia within the subtribe Milleriinae. It is suggested that the present species composition of Guizotia and the subtribal placement of the genus need to be redefined