TESTING FOR UNEQUAL AMOUNTS OF EVOLUTION IN A CONTINUOUS CHARACTER ON DIFFERENT BRANCHES OF A PHYLOGENETIC TREE USING LINEAR AND SQUARED-CHANGE PARSIMONY: AN EXAMPLE USING LESSER ANTILLEAN ANOLIS LIZARDS

Although a large body of work investigating tests of correlated evolution of two continuous characters exists, hypotheses such as character displacement are really tests of whether substantial evolutionary change has occurred on a particular branch or branches of the phylogenetic tree. In this study, we present a methodology for testing such a hypothesis using ancestral character state reconstruction and simulation. Furthermore, we suggest how to investigate the robustness of the hypothesis test by varying the reconstruction methods or simulation parameters. As a case study, we tested a hypothesis of character displacement in body size of Caribbean Anolis lizards. We compared squared-change, weighted squared-change, and linear parsimony reconstruction methods, gradual Brownian motion and speciational models of evolution, and several resolution methods for linear parsimony. We used ancestor reconstruction methods to infer the amount of body size evolution, and tested whether evolutionary change in body size was greater on branches of the phylogenetic tree in which a transition from occupying a single-species island to a two-species island occurred. Simulations were used to generate null distributions of reconstructed body size change. The hypothesis of character displacement was tested using Wilcoxon Rank-Sums. When tested against simulated null distributions, all of the reconstruction methods resulted in more significant P-values than when standard statistical tables were used. These results confirm that P-values for tests using ancestor reconstruction methods should be assessed via simulation rather than from standard statistical tables. Linear parsimony can produce an infinite number of most parsimonious reconstructions in continuous characters. We present an example of assessing the robustness of our statistical test by exploring the sample space of possible resolutions. We compare ACCTRAN and DELTRAN resolutions of ambiguous character reconstructions in linear parsimony to the most and least conservative resolutions for our particular hypothesis.     

A morphometrics‐based phylogeny of the temperate Gondwanan mite harvestmen (Opiliones,Cyphophthalmi, Pettalidae)

A phylogenetic estimation of the temperate Gondwanan mite harvestman family Pettalidae (Arachnida, Opiliones, Cyphophthalmi) was conducted using 143 morphological variables (59 raw and 84 scaled measurements) from 37 ingroup and 15 outgroup terminals. We used custom algorithms to do pairwise comparisons between characters and identify sets of dependent characters, which were collapsed using principal components analysis. We analysed the resulting data without discretization under the parsimony criterion. Monophyly or paraphyly of most groups suspected from previous molecular and morphological phylogenetic studies were recovered. Trees were optimized for monophyly of 20 different focus clades by varying character phylogenetic independence. This yielded a final tree with monophyly of 15 out of 20 focus clades, including the South African pettalids, which contains the troglomorphic species Speleosiro argasiformis Lawrence, 1931. Two of the remaining five clades were found paraphyletic, with the genera Aoraki, Rakaia, and Siro always being found polyphyeletic.     

Compatibility analysis and its applications

A two-state character is defined as uniquely derived if it has only evolved once in the history of a group, without subsequent reversal. Two independent characters cannot both be uniquely derived if all four possible combinations (or all three excluding that of the two ancestral forms) occur.
A number of ways of choosing compatible sets of uniquely derived characters are discussed and used to derive possible unrooted and rooted trees. Results of these are related to those chosen on parsimony criteria, using data for orthopteroid groups, and the assumptions of both methods are compared. Application of compatibility analysis to the moth genera Teldenia and Argodrepana is also discussed. Compatibility and parsimony methods are complementary rather than exclusive of each other.     

Phylogenetic analysis and intraspecific variation: performance of parsimony,likelihood, and distance methods

Intraspecific variation is abundant in all types of systematic characters but is rarely addressed in simulation studies of phylogenetic method performance. We compared the accuracy of 15 phylogenetic methods using simulations to (1) determine the most accurate method(s) for analyzing polymorphic data (under simplified conditions) and (2) test if generalizations about the performance of phylogenetic methods based on previous simulations of fixed (nonpolymorphic) characters are robust to a very different evolutionary model that explicitly includes intraspecific variation. Simulated data sets consisted of allele frequencies that evolved by genetic drift. The phylogenetic methods included eight parsimony coding methods, continuous maximum likelihood, and three distance methods (UPGMA, neighbor joining, and Fitch-Margoliash) applied to two genetic distance measures (Nei's and the modified Cavalli-Sforza and Edwards chord distance). Two sets of simulations were performed. The first examined the effects of different branch lengths, sample sizes (individuals sampled per species), numbers of characters, and numbers of alleles per locus in the eight-taxon case. The second examined more extensively the effects of branch length in the four-taxon, two-allele case. Overall, the most accurate methods were likelihood, the additive distance methods (neighbor joining and Fitch-Margoliash), and the frequency parsimony method. Despite the use of a very different evolutionary model in the present article, many of the results are similar to those from simulations of fixed characters. Similarities include the presence of the "Felsenstein zone," where methods often fail, which suggests that long-branch attraction may occur among closely related species through genetic drift. Differences between the results of fixed and polymorphic data simulations include the following: (1) UPGMA is as accurate or more accurate than nonfrequency parsimony methods across nearly all combinations of branch lengths, and (2) likelihood and the additive distance methods are not positively misled under any combination of branch lengths tested (even when the assumptions of the methods are violated and few characters are sampled). We found that sample size is an important determinant of accuracy and affects the relative success of methods (i.e., distance and likelihood methods outperform parsimony at small sample sizes). Attempts to generalize about the behavior of phylogenetic methods should consider the extreme examples offered by fixed-mutation models of DNA sequence data and genetic-drift models of allele frequencies.     

The historical biogeography of Mammalia

Palaeobiogeographic reconstructions are underpinned by phylogenies, divergence times and ancestral area reconstructions, which together yield ancestral area chronograms that provide a basis for proposing and testing hypotheses of dispersal and vicariance. Methods for area coding include multi-state coding with a single character, binary coding with multiple characters and string coding. Ancestral reconstruction methods are divided into parsimony versus Bayesian/likelihood approaches. We compared nine methods for reconstructing ancestral areas for placental mammals. Ambiguous reconstructions were a problem for all methods. Important differences resulted from coding areas based on the geographical ranges of extant species versus the geographical provenance of the oldest fossil for each lineage. Africa and South America were reconstructed as the ancestral areas for Afrotheria and Xenarthra, respectively. Most methods reconstructed Eurasia as the ancestral area for Boreoeutheria, Euarchontoglires and Laurasiatheria. The coincidence of molecular dates for the separation of Afrotheria and Xenarthra at approximately 100 Ma with the plate tectonic sundering of Africa and South America hints at the importance of vicariance in the early history of Placentalia. Dispersal has also been important including the origins of Madagascar's endemic mammal fauna. Further studies will benefit from increased taxon sampling and the application of new ancestral area reconstruction methods.     

Reconstructing ancestral character states: a critical reappraisal

Using parsimony to reconstruct ancestral character states on a phylogenetic tree has become a popular method for testing ecological and evolutionary hypotheses. Despite its popularity, the assumptions and uncertainties of reconstructing the ancestral states of a single character have received less attention than the much less challenging endeavor of reconstructing phylogenetic trees from many characters. Recent research suggests that parsimony reconstructions are often sensitive to violations of the almost universal assumption of equal probabilities of gains and losses. In addition, maximum likelihood has been developed as an alternative to parsimony reconstruction, and has also revealed a surprising amount of uncertainty in ancestral reconstructions.     

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.     

Phylogeny of the arachnid order Opiliones (Arthropoda) inferred from a combined approach of complete 18S and partial 28S ribosomal DNA sequences and morphology

The phylogenetic relationships among the main evolutionary lines of the arachnid order Opiliones were investigated by means of molecular (complete 18S rDNA and the D3 region of the 28S rDNA genes) and morphological data sets. Equally and differentially weighted parsimony analyses of independent and combined data sets provide evidence for the monophyly of the Opiliones. In all the analyses, the internal relationships of the group coincide in the monophyly of the following main groups: Cyphophthalmi, Eupnoi Palpatores, Dyspnoi Palpatores, and Laniatores. The Cyphophthalmi are monophyletic and sister to a clade that includes all the remaining opilionid taxa (=Phalangida). Within the Phalangida the most supported hypothesis suggests that Palpatores are paraphyletic, as follows: (Eupnoi (Dyspnoi + Laniatores)), but the alternative hypothesis (Laniatores (Eupnoi + Dyspnoi)) is more parsimonious in some molecular data analyses. Relationships within the four main clades are also addressed. Evolution of some morphological characters is discussed, and plesiomorphic states of these characters are evaluated using molecular data outgroup polarization. Finally, Martens' hypothesis of opilionid evolution is assessed in relation to our results.     

Systematics and morphological evolution within the moss family Bryaceae: a comparison between parsimony and Bayesian methods for reconstruction of ancestral character states

The Bryaceae are a large cosmopolitan moss family including genera of significant morphological and taxonomic complexity. Phylogenetic relationships within the Bryaceae were reconstructed based on DNA sequence data from all three genomic compartments. In addition, maximum parsimony and Bayesian inference were employed to reconstruct ancestral character states of 38 morphological plus four habitat characters and eight insertion/deletion events. The recovered phylogenetic patterns are generally in accord with previous phylogenies based on chloroplast DNA sequence data and three major clades are identified. The first clade comprises Bryum bornholmense, B. rubens, B. caespiticium, and Plagiobryum. This corroborates the hypothesis suggested by previous studies that several Bryum species are more closely related to Plagiobryum than to the core Bryum species. The second clade includes Acidodontium, Anomobryum, and Haplodontium, while the third clade contains the core Bryum species plus Imbribryum. Within the latter clade, B. subapiculatum and B. tenuisetum form the sister clade to Imbribryum. Reconstructions of ancestral character states under maximum parsimony and Bayesian inference suggest fourteen morphological synapomorphies for the ingroup and synapomorphies are detected for most clades within the ingroup. Maximum parsimony and Bayesian reconstructions of ancestral character states are mostly congruent although Bayesian inference shows that the posterior probability of ancestral character states may decrease dramatically when node support is taken into account. Bayesian inference also indicates that reconstructions may be ambiguous at internal nodes for highly polymorphic characters.     

Inferring the ancestral sexuality and reproductive condition in sponges (Porifera)

Considerable diversity abounds among sponges with respect to reproductive and developmental biology. Their ancestral sexual mode (gonochorism vs. hermaphroditism) and reproductive condition (oviparity vs. viviparity) however remain unclear, and these traits appear to have undergone correlated evolution in the phylum. To infer ancestral traits and investigate this putative correlation, we used DNA sequence data from two loci (18S ribosomal RNA and cytochrome c oxidase subunit I) to explore the phylogenetic relationships of 62 sponges whose reproductive traits have been previously documented. Although the inferred tree topologies, using the limited data available, favoured paraphyly of sponges, we also investigated ancestral character‐state reconstruction on a phylogeny with constrained sponge monophyly. Both parsimony‐ and likelihood‐based ancestral state reconstructions indicate that viviparity (brooding) was the likely reproductive mode of the ancestral sponge. Hermaphroditism is favoured over gonochorism as the sexual condition of the sponge ancestor under parsimony, but the reconstruction is ambiguous under likelihood, rendering the ancestry of sexuality unresolved in our study. These results are insensitive to the constraint of sponge monophyly when tracing the reproductive characters using parsimony methods. However, the maximum likelihood analysis of the monophyletic hypothetical tree rendered gonochorism as ancestral for the phylum. A test of trait correlation unambiguously favours the concerted evolution of sexuality and reproductive mode in sponges (hermaphroditism/viviparity, gonochorism/oviparity). Although testing ecological hypotheses for the pattern of sponge reproduction is beyond the scope of our analyses, we postulate that certain physiological constrains might be key causes for the correlation of reproductive characters.     

Testing the accuracy of methods for reconstructing ancestral states of continuous characters.

Many methods are available for estimating ancestral values of continuous characteristics, but little is known about how well these methods perform. Here we compare six methods: linear parsimony, squared-change parsimony, one-parameter maximum likelihood (Brownian motion), two-parameter maximum likelihood (Ornstein-Uhlenbeck process), and independent comparisons with and without branch-length information. We apply these methods to data from 20 morphospecies of Pleistocene planktic Foraminifera in order to estimate ancestral size and shape variables, and compare these estimates with measurements on fossils close to the phylogenetic position of 13 ancestors. No method produced accurate estimates for any variable: estimates were consistently less good as predictors of the observed values than were the averages of the observed values. The two-parameter maximum-likelihood model consistently produces the most accurate size estimates overall. Estimation of ancestral sizes is confounded by an evolutionary trend towards increasing size. Shape showed no trend but was still estimated very poorly: we consider possible reasons. We discuss the implications of our results for the use of estimates of ancestral characteristics.     

Ancestral state reconstruction of body size in the Caniformia (Carnivora, Mammalia): the effects of incorporating data from the fossil record

A recent molecular phylogeny of the mammalian order Carnivora implied large body size as the ancestral condition for the caniform subclade Arctoidea using the distribution of species mean body sizes among living taxa. "Extant taxa-only" approaches such as these discount character state observations for fossil members of living clades and completely ignore data from extinct lineages. To more rigorously reconstruct body sizes of ancestral forms within the Caniformia, body size and first appearance data were collected for 149 extant and 367 extinct taxa. Body sizes were reconstructed for four ancestral nodes using weighted squared-change parsimony on log-transformed body mass data. Reconstructions based on extant taxa alone favored large body sizes (on the order of 10 to 50 kg) for the last common ancestors of both the Caniformia and Arctoidea. In contrast, reconstructions incorporating fossil data support small body sizes (< 5 kg) for the ancestors of those clades. When the temporal information associated with fossil data was discarded, body size reconstructions became ambiguous, demonstrating that incorporating both character state and temporal information from fossil taxa unambiguously supports a small ancestral body size, thereby falsifying hypotheses derived from extant taxa alone. Body size reconstructions for Caniformia, Arctoidea, and Musteloidea were not sensitive to potential errors introduced by uncertainty in the position of extinct lineages relative to the molecular topology, or to missing body size data for extinct members of an entire major clade (the aquatic Pinnipedia). Incorporating character state observations and temporal information from the fossil record into hypothesis testing has a significant impact on the ability to reconstruct ancestral characters and constrains the range of potential hypotheses of character evolution. Fossil data here provide the evidence to reliably document trends of both increasing and decreasing body size in several caniform clades. More generally, including fossils in such analyses incorporates evidence of directional trends, thereby yielding more reliable ancestral character state reconstructions.     

Ancestral inference and the study of codon bias evolution: implications for molecular evolutionary analyses of the Drosophila melanogaster subgroup

Reliable inference of ancestral sequences can be critical to identifying both patterns and causes of molecular evolution. Robustness of ancestral inference is often assumed among closely related species, but tests of this assumption have been limited. Here, we examine the performance of inference methods for data simulated under scenarios of codon bias evolution within the Drosophila melanogaster subgroup. Genome sequence data for multiple, closely related species within this subgroup make it an important system for studying molecular evolutionary genetics. The effects of asymmetric and lineage-specific substitution rates (i.e., varying levels of codon usage bias and departures from equilibrium) on the reliability of ancestral codon usage was investigated. Maximum parsimony inference, which has been widely employed in analyses of Drosophila codon bias evolution, was compared to an approach that attempts to account for uncertainty in ancestral inference by weighting ancestral reconstructions by their posterior probabilities. The latter approach employs maximum likelihood estimation of rate and base composition parameters. For equilibrium and most non-equilibrium scenarios that were investigated, the probabilistic method appears to generate reliable ancestral codon bias inferences for molecular evolutionary studies within the D. melanogaster subgroup. These reconstructions are more reliable than parsimony inference, especially when codon usage is strongly skewed. However, inference biases are considerable for both methods under particular departures from stationarity (i.e., when adaptive evolution is prevalent). Reliability of inference can be sensitive to branch lengths, asymmetry in substitution rates, and the locations and nature of lineage-specific processes within a gene tree. Inference reliability, even among closely related species, can be strongly affected by (potentially unknown) patterns of molecular evolution in lineages ancestral to those of interest.     

Reconstructing phylogenies from allozyme data: comparing method performance with congruence

Allozyme data are widely used to infer the phylogenies of populations and closely-related species. Numerous parsimony, distance, and likelihood methods have been proposed for phylogenetic analysis of these data; the relative merits of these methods have been debated vigorously, but their accuracy has not been well explored. In this study, I compare the performance of 13 phylogenetic methods (six parsimony, six distance, and continuous maximum likelihood) by applying a congruence approach to eight allozyme data sets from the literature. Clades are identified that are supported by multiple data sets other than allozymes (e.g. morphology, DNA sequences), and the ability of different methods to recover these 'known' clades is compared. The results suggest that (1) distance and likelihood methods generally outperform parsimony methods, (2) methods that utilize frequency data tend to perform well, and (3) continuous maximum likelihood is among the most accurate methods, and appears to be robust to violations of its assumptions. These results are in agreement with those from recent simulation studies, and help provide a basis for empirical workers to choose among the many methods available for analysing allozyme characters.     

Approximate methods for estimating the pattern of nucleotide substitution and the variation of substitution rates among sites

We propose two approximate methods (one based on parsimony and one on pairwise sequence comparison) for estimating the pattern of nucleotide substitution and a parsimony-based method for estimating the gamma parameter for variable substitution rates among sites. The matrix of substitution rates that represents the substitution pattern can be recovered through its relationship with the observable matrix of site pattern frequences in pairwise sequence comparisons. In the parsimony approach, the ancestral sequences reconstructed by the parsimony algorithm were used, and the two sequences compared are those at the ends of a branch in the phylogenetic tree. The method for estimating the gamma parameter was based on a reinterpretation of the numbers of changes at sites inferred by parsimony. Three data sets were analyzed to examine the utility of the approximate methods compared with the more reliable likelihood methods. The new methods for estimating the substitution pattern were found to produce estimates quite similar to those obtained from the likelihood analyses. The new method for estimating the gamma parameter was effective in reducing the bias in conventional parsimony estimates, although it also overestimated the parameter. The approximate methods are computationally very fast and appear useful for analyzing large data sets, for which use of the likelihood method requires excessive computation.      

Weighted parsimony outperforms other methods of phylogenetic inference under models appropriate for morphology

One of the lasting controversies in phylogenetic inference is the degree to which specific evolutionary models should influence the choice of methods. Model‐based approaches to phylogenetic inference (likelihood, Bayesian) are defended on the premise that without explicit statistical models there is no science, and parsimony is defended on the grounds that it provides the best rationalization of the data, while refraining from assigning specific probabilities to trees or character‐state reconstructions. Authors who favour model‐based approaches often focus on the statistical properties of the methods and models themselves, but this is of only limited use in deciding the best method for phylogenetic inference—such decision also requires considering the conditions of evolution that prevail in nature. Another approach is to compare the performance of parsimony and model‐based methods in simulations, which traditionally have been used to defend the use of models of evolution for DNA sequences. Some recent papers, however, have promoted the use of model‐based approaches to phylogenetic inference for discrete morphological data as well. These papers simulated data under models already known to be unfavourable to parsimony, and modelled morphological evolution as if it evolved just like DNA, with probabilities of change for all characters changing in concert along tree branches. The present paper discusses these issues, showing that under reasonable and less restrictive models of evolution for discrete characters, equally weighted parsimony performs as well or better than model‐based methods, and that parsimony under implied weights clearly outperforms all other methods.     

追溯山茱萸科植物的性状进化和生物地理学历史——方法选择的影响?

This study compares results on reconstructing the ancestral state of characters and ancestral areas of distribution in Cornaceae to gain insights into the impact of using different analytical methods. Ancestral character state reconstructions were compared among three methods (parsimony, maximum likelihood, and stochastic character mapping) using MESQUITE and a full Bayesian method in BAYESTRAITS and inferences of ancestral area distribution were compared between the parsimony-based dispersal-vicariance analysis (DIVA) and a newly developed maximum likelihood (ML) method. Results indicated that among the six inflorescence and fruit characters examined, "perfect" binary characters (no homoplasy, no polymorphism within terminals, and no missing data) are little affected by choice of method, while homoplasious characters and missing data are sensitive to methods used. Ancestral areas at deep nodes of the phylogeny are substantially different between DIVA and ML and strikingly different between analyses including and excluding fossils at three deepest nodes. These results, while raising caution in making conclusions on trait evolution and historical biogeography using conventional methods, demonstrate a limitation in our current understanding of character evolution and biogeography. The biogeographic history favored by the ML analyses including fossils suggested the origin and early radiation of Cornus likely occurred in the late Cretaceous and earliest Tertiary in Europe and intercontinental disjunctions in three lineages involved movements across the North Atlantic Land Bridge (BLB) in the early and mid Tertiary. This result is congruent with the role of NALB for post-Eocene migration and in connecting tropical floras in North America and Africa, and in eastern Asia and South America. However, alternative hypotheses with an origin in eastern Asia and early Trans-Beringia migrations of the genus cannot be ruled out.     

Grazers, shredders and filtering carnivores--the evolution of feeding ecology in Drusinae (Trichoptera: Limnephilidae): insights from a molecular phylogeny

We examined the phylogenetic relationships between species and genera within the caddisfly subfamily Drusinae (Trichoptera: Limnephilidae) using sequence data from two mitochondrial loci (cytochrome oxidase 1, large subunit rRNA) and one nuclear gene (wingless). Sequence data were analysed for 28 species from five genera from the subfamily. We analysed individual and combined data sets using a Bayesian Markov Chain Monte Carlo and a maximum parsimony approach and compared the performance of each partition for resolving phylogenetic relationships at this level. In terms of resolution and phylogenetic utility wingless outperformed the two mitochondrial gene partitions. Using both Shimodaira-Hasegawa and expected likelihood weights tests we tested several hypotheses of relationships previously inferred based on adult morphological characters. The data did not support the generic concept, or many previously proposed species groupings, based on adult morphology. In contrast, the molecular data correlated with the morphology and feeding ecology of larvae. Using Bayesian ancestral character state reconstructions we inferred the evolution of feeding ecology and relevant larval morphological characters. Our analyses showed that within the subfamily Drusinae two derived feeding types evolved. One of these--grazing epilithic algae--is otherwise unusual in the Limnephilidae and may have promoted the high degree of diversity in the Drusinae.     

Tracking character evolution and biogeographic history through time in Cornaceae - Does choice of methods matter?

This study compares results on reconstructing the ancestral state of characters and ancestral areas of distribution in Cornaceae to gain insights into the impact of using different analytical methods. Ancestral character state reconstructions were compared among three methods (parsimony, maximum likelihood, and stochastic character mapping) using MESQUITE and a full Bayesian method in BAYESTRAITS and inferences of ancestral area distribution were compared between the parsimony-based dispersal-vicariance analysis (DIVA) and a newly developed maximum likelihood (ML) method. Results indicated that among the six inflorescence and fruit char-acters examined, “perfect” binary characters (no homoplasy, no polymorphism within terminals, and no missing data) are little affected by choice of method, while homoplasious characters and missing data are sensitive to methods used. Ancestral areas at deep nodes of the phylogeny are substantially different between DIVA and ML and strikingly different between analyses including and excluding fossils at three deepest nodes. These results, while raising caution in making conclusions on trait evolution and historical biogeography using conventional methods, demonstrate a limitation in our current understanding of character evolution and biogeography. The biogeographic history favored by the ML analyses including fossils suggested the origin and early radiation of Cornus likely occurred in the late Cretaceous and earliest Tertiary in Europe and intercontinental disjunctions in three lineages involved movements across the North Atlantic Land Bridge (BLB) in the early and mid Tertiary. This result is congruent with the role of NALB for post-Eocene migration and in connecting tropical floras in North America and Africa, and in eastern Asia and South America. However, alternative hypotheses with an origin in eastern Asia and early Trans-Beringia migrations of the genus cannot be ruled out.