Something for nothing? Reconstruction of ancestral character states in asterinid sea star development

SUMMARY Traits from early development mapped onto phylogenetic trees can potentially offer insight into the evolutionary history of development by inferring the states of those characters among ancestors at nodes in the phylogeny. A key and often-overlooked aspect of such mapping is the underlying model of character evolution. Without a well-supported and realistic model ("nothing"), character mapping of ancestral traits onto phylogenetic trees might often return results ("something") that lack a sound basis. Here we reconsider a challenging case study in this area of evolutionary developmental biology: the inference of ancestral states for ecological and morphological characters in the reproduction and larval development of asterinid sea stars. We apply improved analytical methods to an expanded set of asterinid phylogenetic data and developmental character states. This analysis shows that the new methods might generally offer some independent insight into choice of a model of character evolution, but that in the specific case of asterinid sea stars the quantitative features of the model (especially the relative probabilities of different directions of change) have an important effect on the results. We suggest caution in applying ancestral state reconstructions in the absence of an independently corroborated model of character evolution, and highlight the need for such modeling in evolutionary developmental biology.     

On phylogenetic tests of irreversible evolution

"Dollo's law" states that, following loss, a complex trait cannot reevolve in an identical manner. Although the law has previously fallen into disrepute, it has only recently been challenged with statistical phylogenetic methods. We employ simulation studies of an irreversible binary character to show that rejections of Dollo's law based on likelihood-ratio tests of transition rate constraints or on reconstructions of ancestral states are frequently incorrect. We identify two major causes of errors: incorrect assignment of root state frequencies, and neglect of the effect of the character state on rates of speciation and extinction. Our findings do not necessarily overturn the conclusions of phylogenetic studies claiming reversals, but we demonstrate devastating flaws in the methods that are the foundation of all such studies. Furthermore, we show that false rejections of Dollo's law can be reduced by the use of appropriate existing models and model selection procedures. More powerful tests of irreversibility require data beyond phylogenies and character states of extant taxa, and we highlight empirical work that incorporates additional information.     

Character Coding and Inapplicable Data

Inapplicable character states occur when character complexes are absent or reduced in some of the taxa. Several approaches have been proposed for representing such states in a character matrix so that the inapplicable condition has no effect on the placement of taxa and/or the applicable states are independent and not redundant. Here we examine each of these approaches and demonstrate that all have shortcomings. Coding inapplicables as “?” (reductive coding), although flawed, is currently the best way to analyze data sets that contain inapplicable character states.     

Greedy Selection of Species for Ancestral State Reconstruction on Phylogenies: Elimination Is Better than Insertion

Accurate reconstruction of ancestral character states on a phylogeny is crucial in many genomics studies. We study how to select species to achieve the best reconstruction of ancestral character states on a phylogeny. We first show that the marginal maximum likelihood has the monotonicity property that more taxa give better reconstruction, but the Fitch method does not have it even on an ultrametric phylogeny. We further validate a greedy approach for species selection using simulation. The validation tests indicate that backward greedy selection outperforms forward greedy selection. In addition, by applying our selection strategy, we obtain a set of the ten most informative species for the reconstruction of the genomic sequence of the so-called boreoeutherian ancestor of placental mammals. This study has broad relevance in comparative genomics and paleogenomics since limited research resources do not allow researchers to sequence the large number of descendant species required to reconstruct an ancestral sequence.     

The distribution and the systematic relevance of the androecial characters oligomery and polymery in the Magnoliophytina

In the system of classification sensu Cronquist, the division between the two subclasses Dilleniidae and Rosidae relies on only a few characters. The major distinguishing character is the occurrence of centripetal or centrifugal sequence of stamen development. We question not only the application and the value of this character, but also the maintenance of the Rosidae and Dilleniidae as separately derived natural taxa, because of the lack of solid distinguishing criteria. We have studied the staminal characters in the 192 families of the Dilleniidae and the Rosidae, using Dahlgren's classification diagram.
The study focussed on the following androecial characters: 1) polyandry; i.e. spiral polyandry, complex polyandry (viz. centrifugal, centripetal and lateral polyandry) and cyclic polyandry 2) diplostemony and obdiplostemony 3) haplostemony and obhaplostemony.
In our proposal only two characters are maintained, namely oligomery (diplostemony s.l.) and polymery (true polyandry s.l.), both having a number of character states. Polyandry, as it has been described up to now, is not to be regarded as homologous in all Magnoliophytina. Both characters correspond to two different lines within the Magnoliatae. In addition, the division between the Dilleniidae and Rosidae sensu Cronquist is no longer tenable.
Finally, a number of modifications of Dahlgren's diagram are suggested and a division of the Magnoliophytina in a polymerous group and an oligomerous group is proposed. The Dilleniidae-Rosidae complex is included in the oligomerous group.     

Identification of apomorphies and the role of groundpatterns in molecular systematics

Putative apomorphic character states are the only relevant phylogenetic signal contained in sets of sequence data. Using the sequence position as a character, a way to identify putative apomorphies prior to phylogenetic analysis is proposed. It is shown that distance-matrix methods use trivial characters. The concept of the asymmetrical split is presented for determination of character polarity. It is furthermore argued that groundpatterns (node sequences) should be reconstructed prior to the study of relationships between taxa of high phylogenetic age. The 'evolutionary noise'contained in groundpatterns can be illustrated with a network of distances using a split-decomposition analysis.     

Fixed Character States and the Optimization of Molecular Sequence Data

A method is proposed to optimize molecular sequence data that does not employ multiple sequence alignment. This method treats entire homologous contiguous stretches of sequence data as individual characters. This sequence is treated as the homologous unit employed in phylogeny reconstruction. The sets of specific sequences exhibited by the terminal taxa constitute the character states. The number of states is then less than or equal to the number of unique sequences (or homologous fragments) exhibited by the data. A matrix of transformation costs is created to relate the states to one another. The cells of this matrix are defined as the minimum transformation cost between each pair of states based on insertion–deletion and base substitution costs. The diagnosis of a topology then follows existing dynamic programming techniques, with the number of states greatly expanded. Since the possible sequences reconstructed at nodes are limited to those exhibited by the terminals, cladograms constructed in this way may be longer than those of other methods in that they require a greater number of weighted evolutionary events. Example data, the effects of missing data, restricted ancestors, and putative long-branch attraction are discussed.     

Comparative and phylogenetic analysis of developmental sequences

Event pairing has been proposed for the optimization of developmental sequences (event sequences) on a given phylogenetic hypothesis (cladogram) to determine instances of sequence heterochrony. Here, we show that event pairing is faulty, leading to the optimization of impossible hypothetical ancestors, the underestimation of the lengths of the developmental sequences on the tree, and the proposition of synapomorphies that are not supported by the data. When used for phylogenetic analysis, event pairing can even produce cladograms that are inconsistent with the data. These errors are caused by the fact that event pairing treats dependent features as if they were independent. We present a new method for comparative and phylogenetic analysis of developmental sequences that does not exhibit these errors. Our method applies Search-based character optimization and treats the entire developmental sequence as a single character that is then analyzed by using an edit cost function, which specifies the transformation cost between pairs of observed and unobserved character states, and dynamic programming. In other words, the developmental sequence is directly optimized on the tree. We used event pairing as an edit cost function, but others are possible.     

Vagaries in the delimitation of character states in quantitative variation--an experimental study

An experimental study on the delimitation of character states in continuous variation indicates that (1) the way data are presented influences the assignment of character states and (2) states in the same data set are delimited in various ways by different individuals. Forty-nine individuals were given a set of graphs denoting variation of 10 characters in the genus Kalmia (Ericaceae) and outgroups, all identification having been removed from the graphs. The variation was represented in one of three ways: as 95% confidence intervals on a linear scale, as 95% confidence intervals on a log10 scale, or with bars showing SD x 2 on a linear scale. No two individuals scored a set of graphs in the same way, and only one character in one representation was scored identically by all individuals; the scoring for this character was completely different when the ordinate was changed from linear to logarithmic. Together, the 49 individuals delimited states within each character between 9 and 16 different ways. In general, variation represented by 2 x SD bars elicited the largest numbers of different scorings, yet with a relatively low number of states; the complexity of the patterns in the graphs in this representation was greatest. Expert knowledge appears to be of dubious value in delimiting states in such variation, and if such characters are to be used in phylogenetic analyses, states could be delimited by people who know nothing of the details of the study being scored; in any case, presentation of data and an explicit protocol to follow when delimiting states are essential. In converting data of this type into character states, psychological factors are particularly likely to come into play. Other implications of our experiments include the severe underdetermination of some phylogenetic hypotheses by observation and the heterogeneous nature of morphological data.     

More is better: the uses of developmental genetic data to reconstruct perianth evolution

The origin and evolution of the perianth remains enigmatic. While it seems likely that an undifferentiated perianth consisting of tepals arose early in angiosperm evolution, it is unclear when and how differentiated perianths consisting of distinct organs, such as petals and sepals, arose. Phylogenetic reconstructions of ancestral perianth states across angiosperms have traditionally relied on morphological data from extant species, but these analyses often produce equivocal results. Here we describe the use of developmental genetic data as an additional strategy to infer the ancestral perianth character state for different angiosperm clades. By assessing functional data in combination with expression data in a maximum likelihood framework, we provide a novel approach for investigating the evolutionary history of the perianth. Results of this analysis provide new insights into perianth evolution and provide a proof of concept for using this strategy to explore the incorporation of developmental genetic data in character state reconstructions. As the assumptions outlined here are tested and more genetic data are generated, we hope that ancestral state reconstructions based on multiple lines of evidence will converge.     

A space-efficient construction of the Burrows-Wheeler transform for genomic data.

Algorithms for exact string matching have substantial application in computational biology. Time-efficient data structures which support a variety of exact string matching queries, such as the suffix tree and the suffix array, have been applied to such problems. As sequence databases grow, more space-efficient approaches to exact matching are becoming more important. One such data structure, the compressed suffix array (CSA), based on the Burrows-Wheeler transform, has been shown to require memory which is nearly equal to the memory requirements of the original database, while supporting common sorts of query problems time efficiently. However, building a CSA from a sequence in efficient space and time is challenging. In 2002, the first space-efficient CSA construction algorithm was presented. That implementation used (1+2 log2 |summation|)(1+epsilon) bits per character (where epsilon is a small fraction). The construction algorithm ran in as much as twice that space, in O(| summation|n log(n)) time. We have created an implementation which can also achieve these asymptotic bounds, but for small alphabets, and only uses 1/2 (1+|summation|)(1+epsilon) bits per character, a factor of 2 less space for nucleotide alphabets. We present time and space results for the CSA construction and querying of our implementation on publicly available genome data which demonstrate the practicality of this approach.     

Two Amyloid States of the Prion Protein Display Significantly Different Folding Patterns

It has been well established that a single amino acid sequence can give rise to several conformationally distinct amyloid states. The extent to which amyloid structures formed within the same sequence are different, however, remains unclear. To address this question, we studied two amyloid states (referred to as R- and S-fibrils) produced in vitro from highly purified full-length recombinant prion protein. Several biophysical techniques including X-ray diffraction, CD, Fourier transform infrared spectroscopy (FTIR), hydrogen-deuterium exchange, proteinase K digestion, and binding of a conformation-sensitive fluorescence dye revealed that R- and S-fibrils have substantially different secondary, tertiary, and quaternary structures. While both states displayed a 4. 8-Å meridional X-ray diffraction typical for amyloid cross-β-spines, they showed markedly different equatorial profiles, suggesting different folding pattern of β-strands. The experiments on hydrogen-deuterium exchange monitored by FTIR revealed that only small fractions of amide protons were protected in R- or S-fibrils, an argument for the dynamic nature of their cross-β-structure. Despite this fact, both amyloid states were found to be very stable conformationally as judged from temperature-induced denaturation monitored by FTIR and the conformation-sensitive dye. Upon heating to 80 °C, only local unfolding was revealed, while individual state-specific cross-β features were preserved. The current studies demonstrated that the two amyloid states formed by the same amino acid sequence exhibited significantly different folding patterns that presumably reflect two different architectures of cross-β-structure. Both S- and R-fibrils, however, shared high conformational stability, arguing that the energy landscape for protein folding and aggregation can contain several deep free-energy minima.     

The phylogeny of amphibian metamorphosis

Frogs have one of the most extreme metamorphoses among vertebrates. How did this metamorphosis evolve? By combining the methods previously proposed by Mabee and Humphries (1993) and Velhagen (1997), I develop a phylogenetic method suited for rigorous analysis of this question. In a preliminary analysis using 12 transformation sequence characters and 36 associated event sequence characters, all drawn from the osteology of the skull, the evolution of metamorphosis is traced on an assumed phylogeny. This phylogeny has lissamphibians (frogs, salamanders, and caecilians) monophyletic, with frogs the sister group of salamanders. Successive outgroups used are temnospondyls and discosauriscids, both of which are fossil groups for which ontogenetic data are available. In the reconstruction of character evolution, an unambiguous change (synapomorphy) along the branch leading to lissamphibians is a delay in the lengthening of the maxilla until metamorphosis, in accordance with my previous suggestion (Reiss, 1996). However, widening of the interpterygoid vacuity does not appear as a synapomophy of lissamphibians, due to variation in the character states in the outgroups. From a more theoretical perspective, the reconstructed evolution of amphibian metamorphosis involves examples of heterochrony, through the shift of ancestral premetamorphic events to the metamorphic period, caenogenesis, through the origin of new larval features, and terminal addition, through the origin of new adult features. Other changes don't readily fit these categories. This preliminary study provides evidence that metamorphic changes in frogs arose as further modifications of changes unique to lissamphibians, as well as a new method by which such questions can be examined.     

ON COMBINING PROTEIN SEQUENCES AND NUCLEIC ACID SEQUENCES IN PHYLOGENETIC ANALYSIS: THE HOMEOBOX PROTEIN CASE

Abstract— Amino acid encoding genes contain character state information that may be useful for phylogenetic analysis on at least two levels. The nucleotide sequence and the translated amino acid sequences have both been employed separately as character states for cladistic studies of various taxa, including studies of the genealogy of genes in multigene families. In essence, amino acid sequences and nucleic acid sequences are two different ways of character coding the information in a gene. Silent positions in the nucleotide sequence (first or third positions in codons that can accrue change without changing the identity of the amino acid that the triplet codes for) may accrue change relatively rapidly and become saturated, losing the pattern of historical divergence. On the other hand, non-silent nucleotide alterations and their accompanying amino acid changes may evolve too slowly to reveal relationships among closely related taxa. In general, the dynamics of sequence change in silent and non-silent positions in protein coding genes result in homoplasy and lack of resolution, respectively. We suggest that the combination of nucleic acid and the translated amino acid coded character states into the same data matrix for phylogenetic analysis addresses some of the problems caused by the rapid change of silent nucleotide positions and overall slow rate of change of non-silent nucleotide positions and slowly changing amino acid positions. One major theoretical problem with this approach is the apparent non-independence of the two sources of characters. However, there are at least three possible outcomes when comparing protein coding nucleic acid sequences with their translated amino acids in a phylogenetic context on a codon by codon basis. First, the two character sets for a codon may be entirely congruent with respect to the information they convey about the relationships of a certain set of taxa. Second, one character set may display no information concerning a phylogenetic hypothesis while the other character set may impart information to a hypothesis. These two possibilities are cases of non-independence, however, we argue that congruence in such cases can be thought of as increasing the weight of the particular phylogenetic hypothesis that is supported by those characters. In the third case, the two sources of character information for a particular codon may be entirely incongruent with respect to phylogenetic hypotheses concerning the taxa examined. In this last case the two character sets are independent in that information from neither can predict the character states of the other. Examples of these possibilities are discussed and the general applicability of combining these two sources of information for protein coding genes is presented using sequences from the homeobox region of 46 homeobox genes fromDrosophila melanogasterto develop a hypothesis of genealogical relationship of these genes in this large multigene family.     

THE EFFECT OF ORDERED CHARACTERS ON PHYLOGENETIC RECONSTRUCTION

Abstract Morphological structures are likely to undergo more than a single change during the course of evolution. As a result, multistate characters are common in systematic studies and must be dealt with. Particularly interesting is the question of whether or not multistate characters should be treated as ordered (additive) or unordered (non-additive). In accepting a particular hypothesis of order, numerous others are necessarily rejected. We review some of the criteria often used to order character states and the underlying assumptions inherent in these criteria.
The effects that ordered multistate characters can have on phylogenetic reconstruction are examined using 27 data sets. It has been suggested that hypotheses of character state order are more informative then hypotheses of unorder and may restrict the number of equally parsimonious trees as well as increase tree resolution. Our results indicate that ordered characters can produce more, equal or less equally parsimonious trees and can increase, decrease or have no effect on tree resolution. The effect on tree resolution can be a simple gain in resolution or a dramatic change in sister-taxa relationships. In cases where several outgroups are included in the data matrix, hypotheses of order can change character polarities by altering outgroup topology. Ordered characters result in a different topology from unordered characters only when the hierarchy of the cladogram disagrees with the investigator's a priori hypothesis of order. If the best criterion for assessing character evolution is congruence with other characters, the practice of ordering multistate characters is inappropriate.     

On the origin of termite workers: weighing up the phylogenetic evidence

Resolving the phylogenetic history of a 'true' worker caste in termites is essential to our understanding of termite eusocial evolution. Whether this caste is ancient and monophyletic or derived and polyphyletic will have a tremendous impact on our interpretation of termite eusocial history and remains an outstanding question in termite biology. Recent work has begun to re-examine this question in light of new phylogenetic information, but new questions have now arisen about how best to model character state changes in termite caste systems. In the present paper, we compare the models of Grandcolas and D'Haese [J. Evol. Biol. 15 (2002) 885] and Thompson et al. [J. Evol Biol. 13 (2000) 8691 and attempt to make explicit how these proposals differ with respect to the number of, and homology between, character states. We highlight the support each model has for the two principal, but competing, evolutionary hypotheses outlined above.     

The evolution of pelvic osteology and soft tissues on the line to extant birds (Neornithes)

Substantial differences in pelvic osteology and soft tissues separate crown group crocodylians (Crocodylia) and birds (Neornithes). A phylogenetic perspective including fossils reveals that these disparities arose in a stepwise pattern along the line to extant birds, with major changes occurring both within and outside Aves. Some character states that preceded the origin of Neornithes are only observable or inferable in extinct taxa. These transitional states are important for recognizing the derived traits of neornithines. Palaeontological and neontological data are vital for reconstructing the sequence of pelvic changes along the line to Neornithes. Soft tissue correlation with osteological structures allows changes in soft tissue anatomy to be traced along a phylogenetic framework, and adds anatomical significance to systematic characters from osteology. Explicitly addressing homologies of bone surfaces reveals many subtleties in pelvic evolution that were previously unrecognized or implicit. I advocate that many anatomical features often treated as independent characters should be interpreted as different character states of the same character. Relatively few pelvic character states are unique to Neornithes. Indeed, many features evolved quite early along the line to Neornithes, blurring the distinction between 'avian' and 'non-avian' anatomy.     

A Neural Computational Model of Incentive Salience

Incentive salience is a motivational property with ‘magnet-like’ qualities. When attributed to reward-predicting stimuli (cues), incentive salience triggers a pulse of ‘wanting’ and an individual is pulled toward the cues and reward. A key computational question is how incentive salience is generated during a cue re-encounter, which combines both learning and the state of limbic brain mechanisms. Learning processes, such as temporal-difference models, provide one way for stimuli to acquire cached predictive values of rewards. However, empirical data show that subsequent incentive values are also modulated on the fly by dynamic fluctuation in physiological states, altering cached values in ways requiring additional motivation mechanisms. Dynamic modulation of incentive salience for a Pavlovian conditioned stimulus (CS or cue) occurs during certain states, without necessarily requiring (re)learning about the cue. In some cases, dynamic modulation of cue value occurs during states that are quite novel, never having been experienced before, and even prior to experience of the associated unconditioned reward in the new state. Such cases can include novel drug-induced mesolimbic activation and addictive incentive-sensitization, as well as natural appetite states such as salt appetite. Dynamic enhancement specifically raises the incentive salience of an appropriate CS, without necessarily changing that of other CSs. Here we suggest a new computational model that modulates incentive salience by integrating changing physiological states with prior learning. We support the model with behavioral and neurobiological data from empirical tests that demonstrate dynamic elevations in cue-triggered motivation (involving natural salt appetite, and drug-induced intoxication and sensitization). Our data call for a dynamic model of incentive salience, such as presented here. Computational models can adequately capture fluctuations in cue-triggered ‘wanting’ only by incorporating modulation of previously learned values by natural appetite and addiction-related states.