Coding polymorphism for phylogeny reconstruction

The methodology of coding polymorphic taxa has received limited attention to date. A search of the taxonomic literature revealed seven types of coding methods. Apart from ignoring polymorphic characters (sometimes called the fixed-only method), two main categories can be distinguished: methods that identify the start of a new character state with the origin of an evolutionary novelty, and methods that identify the new state with the fixation of a novelty. The methods of the first category introduce soft reversals, yielding signals that support cladograms incompatible with true phylogenies. We conclude that coding the plesiomorphy is the method to be preferred, unless the ancestral state is unknown, in which case coding as ambiguous is recommended. This holds for coding polymorphism in species as well as in supraspecific taxa. In this light we remark on methods proposed by previous authors.     

THE PHYLOGENETIC RELATIONSHIPS AND BIOGEOGRAPHY OF TRUE PORPOISES (MAMMALIA: PHOCOENIDAE) BASED ON MORPHOLOGICAL DATA

Prior studies of phylogenetic relationships among phocoenids based on morphology and molecular sequence data conflict and yield unresolved relationships among species. This study evaluates a comprehensive set of cranial, postcranial, and soft anatomical characters to infer interrelationships among extant species and several well-known fossil phocoenids, using two different methods to analyze polymorphic data: polymorphic coding and frequency step matrix. Our phylogenetic results confirmed phocoenid monophyly. The division of Phocoenidae into two subfamilies previously proposed was rejected, as well as the alliance of the two extinct genera Salumiphocaena and Piscolithax with Phocoena dioptrica and Phocoenoides dalli . Extinct phocoenids are basal to all extant species. We also examined the origin and distribution of porpoises within the context of this phylogenetic framework. Phocoenid phylogeny together with available geologic evidence suggests that the early history of phocoenids was centered in the North Pacific during the middle Miocene, with subsequent dispersal into the southern hemisphere in the middle Pliocene. A cooling period in the Pleistocene allowed dispersal of the southern ancestor of Phocoena sinus into the North Pacific (Gulf of California).     

Developmental characters in phylogenetic inference and their absolute timing information

Despite the recent surge of interest in studying the evolution of development, surprisingly little work has been done to investigate the phylogenetic signal in developmental characters. Yet, both the potential usefulness of developmental characters in phylogenetic reconstruction and the validity of inferences on the evolution of developmental characters depend on the presence of such a phylogenetic signal and on the ability of our coding scheme to capture it. In a recent study, we showed, using simulations, that a new method (called the continuous analysis) using standardized time or ontogenetic sequence data and squared-change parsimony outperformed event pairing and event cracking in analyzing developmental data on a reference phylogeny. Using the same simulated data, we demonstrate that all these coding methods (event pairing and standardized time or ontogenetic sequence data) can be used to produce phylogenetically informative data. Despite some dependence between characters (the position of an event in an ontogenetic sequence is not independent of the position of other events in the same sequence), parsimony analysis of such characters converges on the correct phylogeny as the amount of data increases. In this context, the new coding method (developed for the continuous analysis) outperforms event pairing; it recovers a lower proportion of incorrect clades. This study thus validates the use of ontogenetic data in phylogenetic inference and presents a simple coding scheme that can extract a reliable phylogenetic signal from these data.     

A phylogenetic analysis of the tropidurine lizards (Squamata: Tropiduridae), including new characters of squamation and epidermal microstructure

Characters of scale surface microstructure are combined with 'traditional'morphological characters in a phylogenetic analysis of the Tropidurini. Tropidurid lizards show variation in types of coarse and fine scale surface microstructure, in the anatomical distribution of different scale surface features, and in scale organ morphology and distribution. The morphology of the inner surface zone of scales is here described for the first time using scanning electron microscopy. Our phylogeny differs considerably from those proposed in earlier studies. New characters and frequency coding of polymorphic characters help resolve the problematic relationships of several species. Statistical confidence supports recognition of one large cis-Andean and one large trans-Andean clade of species. Based on our results, we synonymize Plesiomicrolophus with Microlophus and Uranoscodon with Tropidurus. The phylogenetic relationships of newly discovered Tropidurus are resolved: T. caltathelys is the sister species of T. melanopleurus; T. xanthochilus is the sister species of T. spinulosus. Tropidurus spinulosus is found to be more closely related to T. strobilurus and a clade of Amazonian species than to T. melanopleurus. The species previously placed in Uracentron are more closely related to species previously placed in Plica than to 77 strobilurus as previously thought.     

QUANTITATIVE CHARACTERS IN PHYLOGENETIC ANALYSIS

Abstract— When analysing phylogcnetic relationships at low laxonomic levels it is often the ease that many of the features that can be used to separate taxa show continuous variation. The theoretical and practical problems for the use of such quantitative characters in phylogenetic analysis arc examined. Three methods of coding continuous data into discrete characters are assessed in detail: simple gap-coding, generalised gap-coding and segment-coding, a form of range-coding. The methods are applied to a data set gathered for Eucatyptus L'Hérit. informal subgenus Symphyomyrtus section Maidenana (Myrtaceae). Each method is capable of distorting relative differences between taxa, but segment-coding produces the least amount of distortion, provided the range of variation of the character is divided into a sufficient number of character states.
Continuous quantitative characters provide data for phylogenetic analysis that arc more noisy than those provided by discrete qualitative characters and should, therefore, only be used when the number of qualitative characters is insufficient for resolution of relationships. The results of such analyses should be recognised as provisional pending the discovery of more readily informative characters.     

Robust discriminant analysis and its application to identify protein coding regions of rice genes

Identification of protein coding regions is fundamentally a statistical pattern recognition problem. Discriminant analysis is a statistical technique for classifying a set of observations into predefined classes and it is useful to solve such problems. It is well known that outliers are present in virtually every data set in any application domain, and classical discriminant analysis methods (including linear discriminant analysis (LDA) and quadratic discriminant analysis (QDA)) do not work well if the data set has outliers. In order to overcome the difficulty, the robust statistical method is used in this paper. We choose four different coding characters as discriminant variables and an approving result is presented by the method of robust discriminant analysis.     

Intraspecific variability and timing in ancestral ecology reconstruction: a test case from the cape flora

Thamnochortus (ca. 32 species) is an ecologically diverse genus of Restionaceae. Restionaceae comprise a major component of the southern African Cape flora, wherein eco-diversification might have been important in the generation of high levels of species richness. In an attempt to reconstruct the macroecological history of Thamnochortus, it was found that standard procedures for character state optimization make two inappropriate assumptions. The first is that ancestors are monomorphic (i.e., ecologically uniform) and the second is that eco-diversification follows, or is slower than, lineage diversification. We demonstrate a variety of coding schemes with which the assumption of monomorphy can be avoided. For unordered discrete ecological characters, presence coding and generalized frequency coding (GFC) are suboptimal because they occasionally yield illogical assignments of no state to ancestors. Polymorphism coding or use of the program DIVA are preferable in this respect but are applicable only with parsimony. For continuous eco-characters (e.g., a rainfall gradient, where individual species occur in ranges), GFC and MaxMin coding provide equally valid solutions to optimizing ranges with parsimony. However, MaxMin can be extended to likelihood approaches and is therefore preferable. With respect to rates and timing, all algorithms currently employed for ancestral ecology reconstruction bias toward slow rates of eco-diversification relative to lineage diversification. An alternative to this bias is provided by DIVA, which biases toward accelerated rates of eco-diversification and thus inferences of ecology-driven speciation. We see no way of choosing between these biases; however, phylogeneticists should be aware of them. Applying these methods to Thamnochortus, we find there to be important differences in details, yet general congruence, regarding the historical ecology of this clade. We infer the most recent common ancestor of Thamnochortus to have been a post-fire resprouting species distributed on rocky, well-drained, sandstone-derived soils at lower-middle elevations, in regions of moderate levels of yearly (primarily winter) rainfall. This species would have been distributed in habitats much like those of the southwestern Cape mountains today. Major ecological trends include shifts to lower rainfall regimes and shifts from sandstone to limestone-derived alkaline soils at lower altitudes.     

Adapting generalized frequency coding to use colour spectra in the determination of phylogenetic relationships: an example with hummingbirds

We introduce a methodology to perform phylogenetic analyses on bird taxa using parsimony based on plumage colour characters obtained as spectra. The generalized frequency coding method was adapted to transform continuous colour spectral data into discrete variables. As an example of the application of the methodology, we analysed two species of the hummingbird genus Campylopterus and all species of the genera Anthracothorax and Eulampis , which occur in South and Central America and the Caribbean islands. The genus Topaza was used as an outgroup because of its basal position within the hummingbirds. The results were mostly in accordance with both traditional morphological as well as molecular analyses indicating that there is phylogenetic information in plumage colours measured in an integral and objective way and that the new method is appropriate for phylogenetic analyses at this lower taxonomic level.     

Phylogeny and biogeography of the extant species of triplespine fishes (Triacanthidae, Tetraodontiformes)

A new phylogenetic hypothesis for the living species of triplespine fishes of the Indo-Western Pacific family Triacanthidae (Tetraodontiformes, Teleostei) is proposed. A data set of 55 morphological characters (34 osteological and 21 morphometric) was constructed. A cladistic analysis of the osteological data set yielded a single most-parsimonious tree. This cladogram does not support the monophyly of one of the four genera, Tripodichthys , but Bremer values for this analysis are low. The osteological data set was then combined with a data set of 21 morphometric characters that had previously been used to diagnose the four genera. The analysis of the combined data set produced the same phylogenetic hypothesis, but with greater nodal support. The biogeographical distribution of the living species is then interpreted with the use of this new phylogenetic information.     

Algebraic regularities of the ratio of heterozygosity with a mean and variance of quantitative character

The relationships of heterozygosity with the mean and variance of quantitative character were considered under neutrality, additivity and overdominance of polymorphic loci. Attention was drawn to dependence of the patterns of relationships on the number of polymorphic loci (which varied from 1 to 10) and on the type of polymorphic loci, both homogeneous (polymorphic loci are of the same type) and heterogeneous (polymorphic loci are of the two types) samples of 10 polymorphic loci and their combination. It is shown that increase in the number of polymorphic loci is accompanied with extension of the limits of corresponding relations, whereas the patterns of these relations depend on the type of connection of separate polymorphic locus with the quantitative character and on the ratio of different loci in the set of polymorphic loci. It is assumed that the relationship of heterozygosity with quantitative characters, displayed in the number of experimental works, may contain a component mediated by similar statistical effects. It is inferred that the discrepancy between different authors' experimental data on existence or lack of some relationships between multiplicative heterozygosity and morphological variability of quantitative characters can be explained by different types of relations of polymorphic loci to quantitative characters encountered in their works.     

Artifacts of Coding Amino Acids and Other Composite Characters for Phylogenetic Analysis

A phylogenetic analysis can be no better than the characters on which it is based. Just as it is inappropriate to code character states of individual characters as separate presence/absence characters, it is inappropriate to combine independent characters because not all information in the data is being utilized. Composite characters link otherwise discernible states from different characters together to form new character states. There are two related problems with this coding. First, there is a loss of hierarchic information between the reductive and composite characters when unordered states are used. Second, the linking of separate characters that occurs during the construction of composite character states can create putative synapomorphies that were not present in the separate characters. For amino acid characters, the problem may occur whenever more than one position of a codon is variable among the terminals sampled. Groups that are resolved as paraphyletic with reductive coding may be resolved as monophyletic with composite coding. The artificial character states indicated by the amino acid characters are unlikely to be congruent with the true gene tree.     

POLYMORPHIC TAXA, MISSING VALUES AND CLADISTIC ANALYSIS

Abstract Missing values have been used in cladistic analyses when data are unavailable, inapplicable or sometimes when character states are variable within terminal taxa. The practice of scoring taxa as having "missing values" for polymorphic characters introduces errors into the calculation of cladogram lengths and consistency indices because some character change is hidden within terminals. Because these hidden character steps are not counted, the set of most parsimonious cladograms may differ from those that would be found if polymorphic taxa had been broken into monomorphic subunits. In some cases, the trees found when polymorphisms are scored as missing values may not include any of the most parsimonious trees found when the data are scored properly. Additionally, in some cases, polymorphic taxa may be found to be polyphyletic when broken into monomorphic subunits; this is undetected when polymorphisms are treated as missing. Because of these problems, terminal units in cladistic analysis should be based on unique, fixed combinations of characters. Polymorphic taxa should be subdivided into subunits that are monomorphic for each character used in the analysis. Disregarding errors in topology, the additional hidden steps in a cladogram in which polymorphisms are scored as missing can be calculated by a simple formula, based on the observation that if it is assumed that polymorphic terminals include all combinations of character states, 2 ^p − 1 additional steps are required for each taxon in which p polymorphic binary characters are scored as missing values. Thus, when several polymorphisms are scored as missing in the same taxon, very large errors can be introduced into the calculation of tree length.     

Heterosis, marker mutational processes and population inbreeding history.

Genotype-fitness correlations (GFC) have previously been studied using allozyme markers and have often focused on short-term processes such as recent inbreeding. Thus, models of GFC usually neglect marker mutation and only use heterozygosity as a genotypic index. Recently, GFC have also been reported (i) with DNA markers such as microsatellites, characterized by high mutation rates and specific mutational processes and (ii) using new individual genotypic indices assumed to be more precise than heterozygosity. The aim of this article is to evaluate the theoretical impact of marker mutation on GFC. We model GFC due to short-term processes generated by the current breeding system (partial selfing) and to long-term processes generated by past population history (hybridization). Various mutation rates and mutation models corresponding to different kinds of molecular markers are considered. Heterozygosity is compared to other genotypic indices designed for specific marker types. Highly mutable markers (such as microsatellites) are particularly suitable for the detection of GFC that evolve in relation to short-term processes, whereas GFC due to long-term processes are best observed with intermediate mutation rates. Irrespective of the marker type and population scenario, heterozygosity usually provides higher correlations than other genotypic indices under most biologically plausible conditions.     

New methods to identify conserved microsatellite loci and develop primer sets of high cross-species utility - as demonstrated for birds

We have developed a new approach to create microsatellite primer sets that have high utility across a wide range of species. The success of this method was demonstrated using birds. We selected 35 avian EST microsatellite loci that had a high degree of sequence homology between the zebra finch Taeniopygia guttata and the chicken Gallus gallus and designed primer sets in which the primer bind sites were identical in both species. For 33 conserved primer sets, on average, 100% of loci amplified in each of 17 passerine species and 99% of loci in five non-passerine species. The genotyping of four individuals per species revealed that 24-76% (mean 48%) of loci were polymorphic in the passerines and 18-26% (mean 21%) in the non-passerines. When at least 17 individuals were genotyped per species for four Fringillidae finch species, 71-85% of loci were polymorphic, observed heterozygosity was above 0.50 for most loci and no locus deviated significantly from Hardy-Weinberg proportions. This new set of microsatellite markers is of higher cross-species utility than any set previously designed. The loci described are suitable for a range of applications that require polymorphic avian markers, including paternity and population studies. They will facilitate comparisons of bird genome organization, including genome mapping and studies of recombination, and allow comparisons of genetic variability between species whilst avoiding ascertainment bias. The costs and time to develop new loci can now be avoided for many applications in numerous species. Furthermore, our method can be readily used to develop microsatellite markers of high utility across other taxa.     

Spectral partitioning of phylogenetic data sets based on compatibility

We describe two new methods to partition phylogenetic data sets of discrete characters based on pairwise compatibility. The partitioning methods make no assumptions regarding the phylogeny, model of evolution, or characteristics of the data. The methods first build a compatibility graph, in which each node represents a character in the data set. Edges in the compatibility graph may represent strict compatibility of characters or they may be weighted based on a fractional compatibility scoring procedure that measures how close the characters are to being compatible. Given the desired number of partitions, the partitioning methods then seek to cluster the characters with the highest average pairwise compatibility, so that characters in each cluster are more compatible with each other than they are with characters in the other cluster(s). Partitioning according to these criteria is computationally intractable (NP-hard); however, spectral methods can quickly provide high-quality solutions. We demonstrate that the spectral partitioning effectively identifies characters with different evolutionary histories in simulated data sets, and it is better at highlighting phylogenetic conflict within empirical data sets than previously used partitioning methods.     

"Cryptic" repeating triplets of purines and pyrimidines (cRRY(i)) are frequent and polymorphic: analysis of coding cRRY(i) in the proopiomelanocortin (POMC) and TATA-binding protein (TBP) genes.

Triplets of the form of purine, purine, pyrimidine (RRY(i)) are enhanced in frequency in the genomes of primates, rodents, and bacteria. Some RRY(i) are "cryptic" repeats (cRRY(i)) in which no one tandem run of a trinucleotide predominates. A search of human GenBank sequence revealed that the sequences of cRRY(i) are highly nonrandom. Three randomly chosen human cRRY(i) were sequenced in search of polymorphic alleles. Multiple polymorphic alleles were found in cRRY(i) in the coding regions of the genes for proopiomelanocortin (POMC) and TATA-binding protein (TBP). The highly polymorphic TBP cRRY(i) was characterized in detail. Direct sequencing of 157 unrelated human alleles demonstrated the presence of 20 different alleles which resulted in 29-40 consecutive glutamines in the amino-terminal region of TBP. These alleles are differentially distributed among the races. PCR was used to screen 1,846 additional alleles in order to characterize more fully the range of variation in the population. Three additional alleles were discovered, but there was no example of a substantial sequence amplification as is seen in the repeat sequences associated with X-linked spinal and bulbar muscular atrophy, myotonic dystrophy, or the fragile-X syndrome. The structure of the TBP cRRY(i) is conserved in the five monkey species examined. In the chimpanzee, examination of four individuals revealed that the cRRY(i) was highly polymorphic, but the pattern of polymorphism differed from that in humans. The TBP cRRY(i) displays both similarities with and differences from the previously described RRY(i) in the coding sequence of the androgen receptor. Our data suggest how simple tandem repeats could evolve from cryptic repeats.     

Character analysis in morphological phylogenetics: problems and solutions

Many aspects of morphological phylogenetics are controversial in the theoretical systematics literature and yet are often poorly explained and justified in empirical studies. In this paper, I argue that most morphological characters describe variation that is fundamentally quantitative, regardless of whether they are coded qualitatively or quantitatively by systematists. Given this view, three fundamental problems in morphological character analysis (definition, delimitation, and ordering of character states) may have a common solution: coding morphological characters as continuous quantitative traits. A new parsimony method (step-matrix gap-weighting, a modification of Thiele's approach) is proposed that allows quantitative traits to be analyzed as continuous variables. The problem of scaling or weighting quantitative characters relative to qualitative characters (and to each other) is reviewed, and three possible solutions are described. The new coding method is applied to data from hoplocercid lizards, and the results show the sensitivity of phylogenetic conclusions to different scaling methods. Although some authors reject the use of continuous, overlapping, quantitative characters in phylogenetic analysis, quantitative data from hoplocercid lizards that are coded using the new approach contain significant phylogenetic structure and exhibit levels of homoplasy similar to those seen in data that are coded qualitatively.     

A general approach to proving the minimality of phylogenetic trees illustrated by an example with a set of 23 vertebrates

Summary We have recently described a method of building phylogenetic trees and have outlined an approach for proving whether a particular tree is optimal for the data used. In this paper we describe in detail the method of establishing lower bounds on the length of a minimal tree by partitioning the data set into subsets. All characters that could be involved in duplications in the data are paired with all other such characters. A matching algorithm is then used to obtain the pairing of characters that reveals the most duplications in the data. This matching may still not account for all nucleotide substitutions on the tree. The structure of the tree is then used to help select subsets of three or more. characters until the lower bound found by partitioning is equal to the length of the tree. The tree must then be a minimal tree since no tree can exist with a length less than that of the lower bound.The method is demonstrated using a set of 23 vertebrate cytochrome c sequences with the criterion of minimizing the total number of nucleotide substitutions. There are 131130 7045768798 9603440625 topologically distinct trees that can be constructed from this data set. The method described in this paper does identify 144 minimal tree variants. The method is general in the sense that it can be used for other data and other criteria of length. It need not however always be possible to prove a tree minimal but the method will give an upper and lower bound on the length of minimal trees.