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.