Rethinking pathways to the dioecy-polyploidy association: Genera with many dioecious species have fewer polyploids.

PREMISE: Numerous studies have found a positive association between dioecy and polyploidy; however, this association presents a theoretical conflict: While polyploids are predicted to benefit from self-reproduction for successful establishment, dioecious species cannot self-reproduce. We propose a theoretical framework to resolve this apparent conflict. We hypothesize that the inability of dioecious species to self-reproduce hinders their establishment as polyploids. We therefore expect that genera with many dioecious species have fewer polyploids, leading to a negative association between polyploidy and dioecy across genera. METHODS: We used three publicly available databases to determine ploidy and sexual systems for 131 genera and 546 species. We quantified (1) the relationship between the frequency of polyploid species and the frequency of dioecious species across genera, and (2) the proportion of polyploids with hermaphroditism and dioecy across species, adjusting for phylogenetic history. RESULTS: Across genera, we found a negative relationship between the proportion of polyploids and the proportion of dioecious species, a consistent trend across clades. Across all species, we found that sexual system (dioecious or not) was not associated with polyploidy. CONCLUSIONS: Polyploids are rare in genera in which the majority of species are dioecious, consistent with the theory that self-reproduction favors polyploid establishment. The low frequency of polyploidy among dioecious species indicates the association is not as widespread as previously suggested. Our findings are consistent with previous studies identifying a positive relationship between the two traits, but only if polyploidy promotes a transition to dioecy, and not the reverse.

Disproportionate extinction of South American mammals drove the asymmetry of the Great American Biotic Interchange.

The interchange between the previously disconnected faunas of North and South America was a massive experiment in biological invasion. A major gap in our understanding of this invasion is why there was a drastic increase in the proportion of mammals of North American origin found in South America. Four nonmutually exclusive mechanisms may explain this asymmetry: 1) Higher dispersal rate of North American mammals toward the south, 2) higher origination of North American immigrants in South America, 3) higher extinction of mammals with South American origin, and 4) similar dispersal rate but a larger pool of native taxa in North versus South America. We test among these mechanisms by analyzing ∼20,000 fossil occurrences with Bayesian methods to infer dispersal and diversification rates and taxonomic selectivity of immigrants. We find no differences in the dispersal and origination rates of immigrants. In contrast, native South American mammals show higher extinction. We also find that two clades with North American origin (Carnivora and Artiodactyla) had significantly more immigrants in South America than other clades. Altogether, the asymmetry of the interchange was not due to higher origination of immigrants in South America as previously suggested, but resulted from higher extinction of native taxa in southern South America. These results from one of the greatest biological invasions highlight how biogeographic processes and biotic interactions can shape continental diversity.

Target sequence capture of Barnadesioideae (Compositae) demonstrates the utility of low coverage loci in phylogenomic analyses.

Target sequence capture has emerged as a powerful method to sequence hundreds or thousands of genomic regions in a cost- and time-efficient approach. In most cases, however, targeted regions lack full sequence information for certain samples, due to taxonomic, laboratory, or stochastic factors. Loci lacking molecular data for a large number of samples are commonly excluded from downstream analyses, even though they may still contain valuable information. On the other hand, including data-poor loci may bias phylogenetic analyses. Here we use a target sequence capture dataset of an ecologically and taxonomically diverse group of spiny sunflowers (Asteraceae, or Compositae: Barnadesioideae) to test how the inclusion or exclusion of such data-poor loci affects phylogenetic inference. We investigate the sensitivity of concatenation and coalescent approaches to missing data with matrices of varying taxonomic completeness by filtering loci with different proportions of missing samples prior to data analysis. We find that missing data affect both the topology and branch support of the resulting phylogenies. The matrix containing all loci yielded the overall highest node support values, independently of the amount of missing nucleotides. These results provide empirical support to earlier suggestions based on single genes and data simulations that taxa with high amounts of missing data should not be readily dismissed as they can provide essential information for phylogenomic reconstruction.

The origin of modern patterns of continental diversity in Mauritiinae palms: the Neotropical museum and the Afrotropical graveyard.

While the latitudinal diversity gradient has received much attention, biodiversity and species richness also vary between continents across similar latitudes. Fossil information can be used to understand the evolutionary mechanisms that generated such variation between continents of similar latitudes. We integrated fossil data into a phylogenetic analysis of the Mauritiinae palms, whose extant diversity is restricted to the Neotropics, but extended across Africa and India during most of the Cenozoic. Mauritiinae diverged from its sister lineage Raphiinae ca 106 Ma. Using ancestral state estimation and a lineage through time analysis, we found that diversity arose globally during the late Cretaceous and Palaeocene across South America, Africa and India. The Palaeocene-Eocene transition (ca 56 Ma) marked the end of global Mauritiinae expansion, and the beginning of their decline in both Africa and India. Mauritiinae disappeared from the Indian subcontinent and Africa at the end of the Eocene and the Miocene, respectively. By contrast, Neotropical diversity steadily increased over the last 80 Myr. Taken together, our results suggest that the Neotropics functioned as a continental-scale refuge for Mauritiinae palms, where lineages survived and diversified while global climatic changes that drastically reduced rainforests led to their demise on other continents.

A bioinformatic platform to integrate target capture and whole genome sequences of various read depths for phylogenomics.

The increasing availability of short-read whole genome sequencing (WGS) provides unprecedented opportunities to study ecological and evolutionary processes. Although loci of interest can be extracted from WGS data and combined with target sequence data, this requires suitable bioinformatic workflows. Here, we test different assembly and locus extraction strategies and implement them into secapr, a pipeline that processes short-read data into multilocus alignments for phylogenetics and molecular ecology analyses. We integrate the processing of data from low-coverage WGS (<30×) and target sequence capture into a flexible framework, while optimizing de novo contig assembly and loci extraction. Specifically, we test different assembly strategies by contrasting their ability to recover loci from targeted butterfly protein-coding genes, using four data sets: a WGS data set across different average coverages (10×, 5× and 2×) and a data set for which these loci were enriched prior to sequencing via target sequence capture. Using the resulting de novo contigs, we account for potential errors within contigs and infer phylogenetic trees to evaluate the ability of each assembly strategy to recover species relationships. We demonstrate that choosing multiple sizes of kmer simultaneously for assembly results in the highest yield of extracted loci from de novo assembled contigs, while data sets derived from sequencing read depths as low as 5× recovers the expected species relationships in phylogenetic trees. By making the tested assembly approaches available in the secapr pipeline, we hope to inspire future studies to incorporate complementary data and make an informed choice on the optimal assembly strategy.

Biogeographic origins of southern African Silene (Caryophyllaceae).

Silene (Caryophyllaceae) is distributed predominantly in the northern Hemisphere, where it is most diverse around the Mediterranean Basin. The genus is also well represented in North Africa, extending into tropical, sub-Saharan and southern Africa. Eight native species are recognized in southern Africa, taxonomically placed in two sections: Elisanthe and Silene s.l. Although the taxonomy of the southern African taxa has recently been revised, their phylogenetic relationships and biogeographic history remain unclear. This study aims to infer the phylogenetic position and geographic origins of the southern African taxa. We generated DNA sequences of nuclear and plastid loci from several individuals belonging to all eight species of Silene recognized from southern Africa, and combined our DNA sequences with existing data representing species from major clades (i.e. sections) based on the recently revised Silene infrageneric taxonomy. We used a Bayesian coalescent species tree continuous diffusion approach to co-estimate the species tree and the ancestral areas of representative members of the genus. Our results show that the perennial southern African members of section Elisanthe form a strongly-supported clade with the Eurasian annual S. noctiflora and the Central Asian perennial S. turkestanica. The rest of the perennial species form a strongly-supported clade together with the annual S. aethiopica, which is nested in a larger Mediterranean clade comprising mostly annual species classified in section Silene s.l. Estimates of ancestral areas indicate a late Pleistocene dispersal to southern Africa from central and East Africa for the sub-Saharan members of section Silene s.l. The Elisanthe clade is inferred to have colonized southern Africa through long-distance dispersal from Eurasia during the late Pleistocene. Our findings support the hypothesis of a relatively recent colonization into southern Africa resulting from two independent dispersal events during the Pleistocene.

In situ radiation explains the frequency of dioecious palms on islands.

BACKGROUND AND AIMS: Dioecy has evolved up to 5000 times in angiosperms, despite the potentially high intrinsic costs to unisexuality. Dioecy prevents inbreeding, which is especially relevant on isolated islands when gene pools are small. Dioecy is also associated with certain dispersal traits, such as fruit size and type. However, the influence of dioecy on other life history traits and island distribution remains poorly understood. Here, we test the effect of dioecy on palm (Arecaceae) speciation rates, fruit size and frequency on islands. METHODS: We used phylogenetic comparative methods to estimate the ancestral state of the sexual system and its impact on speciation rates and fruit size. Frequency of sexual systems, effect of insularity on the probability of being dioecious, and phylogenetic clustering of island dioecious vs. mainland species were inferred. Lastly, we determined the interplay of insularity and sexual system on speciation rates. KEY RESULTS: Palms repeatedly evolved different sexual systems (dioecy, monoecy and polygamy) from a hermaphrodite origin. Differences in speciation rates and fruit size among the different sexual systems were not identified. An effect of islands on the probability of the palms being dioecious was also not found. However, we found a high frequency and phylogenetic clustering of dioecious palms on islands, which were not correlated with higher speciation rates. CONCLUSIONS: The high frequency and phylogenetic clustering may be the result of in situ radiation and suggest an 'island effect' for dioecious palms, which was not explained by differential speciation rates. This island effect also cannot be attributed to long-distance dispersal due to the lack of fruit size difference among sexual systems, and particularly because palm dispersal to islands is highly constrained by the interaction between the sizes of fruit and frugivores. Taken together, we suggest that trait flexibility in sexual system evolution and the in situ radiation of dioecious lineages are the underlying causes of the outstanding distribution of palms on islands.

Amazonia is the primary source of Neotropical biodiversity.

The American tropics (the Neotropics) are the most species-rich realm on Earth, and for centuries, scientists have attempted to understand the origins and evolution of their biodiversity. It is now clear that different regions and taxonomic groups have responded differently to geological and climatic changes. However, we still lack a basic understanding of how Neotropical biodiversity was assembled over evolutionary timescales. Here we infer the timing and origin of the living biota in all major Neotropical regions by performing a cross-taxonomic biogeographic analysis based on 4,450 species from six major clades across the tree of life (angiosperms, birds, ferns, frogs, mammals, and squamates), and integrate >1.3 million species occurrences with large-scale phylogenies. We report an unprecedented level of biotic interchange among all Neotropical regions, totaling 4,525 dispersal events. About half of these events involved transitions between major environmental types, with a predominant directionality from forested to open biomes. For all taxonomic groups surveyed here, Amazonia is the primary source of Neotropical diversity, providing >2,800 lineages to other regions. Most of these dispersal events were to Mesoamerica (∼1,500 lineages), followed by dispersals into open regions of northern South America and the Cerrado and Chaco biomes. Biotic interchange has taken place for >60 million years and generally increased toward the present. The total amount of time lineages spend in a region appears to be the strongest predictor of migration events. These results demonstrate the complex origin of tropical ecosystems and the key role of biotic interchange for the assembly of regional biotas.

Selective extinction against redundant species buffers functional diversity.

The extinction of species can destabilize ecological processes. A way to assess the ecological consequences of species loss is by examining changes in functional diversity. The preservation of functional diversity depends on the range of ecological roles performed by species, or functional richness, and the number of species per role, or functional redundancy. However, current knowledge is based on short timescales and an understanding of how functional diversity responds to long-term biodiversity dynamics has been limited by the availability of deep-time, trait-based data. Here, we compile an exceptional trait dataset of fossil molluscs from a 23-million-year interval in the Caribbean Sea (34 011 records, 4422 species) and develop a novel Bayesian model of multi-trait-dependent diversification to reconstruct mollusc (i) diversity dynamics, (ii) changes in functional diversity, and (iii) extinction selectivity over the last 23 Myr. Our results identify high diversification between 23-5 Mya, leading to increases in both functional richness and redundancy. Conversely, over the last three million years, a period of high extinction rates resulted in the loss of 49% of species but only 3% of functional richness. Extinction rates were significantly higher in small, functionally redundant species suggesting that competition mediated the response of species to environmental change. Taken together, our results identify long-term diversification and selective extinction against redundant species that allowed functional diversity to grow over time, ultimately buffering the ecological functions of biological communities against extinction.

Soil fertility and flood regime are correlated with phylogenetic structure of Amazonian palm communities.

BACKGROUND AND AIMS: Identifying the processes that generate and maintain biodiversity requires understanding of how evolutionary processes interact with abiotic conditions to structure communities. Edaphic gradients are strongly associated with floristic patterns but, compared with climatic gradients, have received relatively little attention. We asked (1) How does the phylogenetic composition of palm communities vary along edaphic gradients within major habitat types? and (2) To what extent are phylogenetic patterns determined by (a) habitat specialists, (b) small versus large palms, and (c) hyperdiverse genera? METHODS: We paired data on palm community composition from 501 transects of 0.25 ha located in two main habitat types (non-inundated uplands and seasonally inundated floodplains) in western Amazonian rain forests with information on soil chemistry, climate, phylogeny and metrics of plant size. We focused on exchangeable base concentration (cmol+ kg-1) as a metric of soil fertility and a floristic index of inundation intensity. We used a null model approach to quantify the standard effect size of mean phylogenetic distance for each transect (a metric of phylogenetic community composition) and related this value to edaphic variables using generalized linear mixed models, including a term for spatial autocorrelation. KEY RESULTS: Overall, we recorded 112 008 individuals belonging to 110 species. Palm communities in non-inundated upland transects (but not floodplain transects) were more phylogenetically clustered in areas of low soil fertility, measured as exchangeable base concentration. In contrast, floodplain transects with more severe flood regimes (as inferred from floristic structure) tended to be phylogenetically clustered. Nearly half of the species recorded (44 %) were upland specialists while 18 % were floodplain specialists. In both habitat types, phylogenetic clustering was largely due to the co-occurrence of small-sized habitat specialists belonging to two hyperdiverse genera (Bactris and Geonoma). CONCLUSIONS: Edaphic conditions are associated with the phylogenetic community structure of palms across western Amazonia, and different factors (specifically, soil fertility and inundation intensity) appear to underlie diversity patterns in non-inundated upland versus floodplain habitats. By linking edaphic gradients with palm community phylogenetic structure, our study reinforces the need to integrate edaphic conditions in eco-evolutionary studies in order to better understand the processes that generate and maintain tropical forest diversity. Our results suggest a role for edaphic niche conservatism in the evolution and distribution of Amazonian palms, a finding with potential relevance for other clades.

Toward a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa.

Rapidly growing biological data-including molecular sequences and fossils-hold an unprecedented potential to reveal how evolutionary processes generate and maintain biodiversity. However, researchers often have to develop their own idiosyncratic workflows to integrate and analyze these data for reconstructing time-calibrated phylogenies. In addition, divergence times estimated under different methods and assumptions, and based on data of various quality and reliability, should not be combined without proper correction. Here we introduce a modular framework termed SUPERSMART (Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa), and provide a proof of concept for dealing with the moving targets of evolutionary and biogeographical research. This framework assembles comprehensive data sets of molecular and fossil data for any taxa and infers dated phylogenies using robust species tree methods, also allowing for the inclusion of genomic data produced through next-generation sequencing techniques. We exemplify the application of our method by presenting phylogenetic and dating analyses for the mammal order Primates and for the plant family Arecaceae (palms). We believe that this framework will provide a valuable tool for a wide range of hypothesis-driven research questions in systematics, biogeography, and evolution. SUPERSMART will also accelerate the inference of a "Dated Tree of Life" where all node ages are directly comparable. [Bayesian phylogenetics; data mining; divide-and-conquer methods; GenBank; multilocus multispecies coalescent; next-generation sequencing; palms; primates; tree calibration.].

The road to evolutionary success: insights from the demographic history of an Amazonian palm.

Evolutionary success, as demonstrated by high abundance and a wide geographical range, is related to genetic variation and historical demography. Here we assess how climatic change during the Quaternary influenced the demography and distribution of the Neotropical swamp palm Mauritia flexuosa. Using microsatellite loci and coalescent analyses we examined how demographical dynamics affected genetic diversity, effective population size and connectivity through time and space. Mauritia flexuosa presents significant genetic differentiation between the Amazonian and Cerrado biomes and among different river basins. Amazonian lineages are ancient compared to lineages from the Cerrado, a pattern corroborated using the fossil pollen record, where the species was absent from the Cerrado during the cold and dry periods of the last glacial cycles, then returned during the wet, interglacial phases. Coalescent simulations show that the pattern of observed genetic diversity for M. flexuosa is most likely due to a range retraction during the Last Glacial Maximum, leading to multiple refugia and resulting in high differentiation between Amazonian and Cerrado biomes. Isolation-by-distance and by-environment also shaped the distribution and evolutionary success of M. flexuosa. Our study provides new insights into the historical factors that affected geographical distribution and structure genetic diversity, contributing to long-term evolutionary success.

Biological evidence supports an early and complex emergence of the Isthmus of Panama.

The linking of North and South America by the Isthmus of Panama had major impacts on global climate, oceanic and atmospheric currents, and biodiversity, yet the timing of this critical event remains contentious. The Isthmus is traditionally understood to have fully closed by ca. 3.5 million years ago (Ma), and this date has been used as a benchmark for oceanographic, climatic, and evolutionary research, but recent evidence suggests a more complex geological formation. Here, we analyze both molecular and fossil data to evaluate the tempo of biotic exchange across the Americas in light of geological evidence. We demonstrate significant waves of dispersal of terrestrial organisms at approximately ca. 20 and 6 Ma and corresponding events separating marine organisms in the Atlantic and Pacific oceans at ca. 23 and 7 Ma. The direction of dispersal and their rates were symmetrical until the last ca. 6 Ma, when northern migration of South American lineages increased significantly. Variability among taxa in their timing of dispersal or vicariance across the Isthmus is not explained by the ecological factors tested in these analyses, including biome type, dispersal ability, and elevation preference. Migration was therefore not generally regulated by intrinsic traits but more likely reflects the presence of emergent terrain several millions of years earlier than commonly assumed. These results indicate that the dramatic biotic turnover associated with the Great American Biotic Interchange was a long and complex process that began as early as the Oligocene-Miocene transition.

Endemic palm species shed light on habitat shifts and the assembly of the Cerrado and Restinga floras.

Species expansions into new habitats are often associated with physiological adaptations, for instance when rain forest lineages colonize dry habitats. Although such shifts have been documented for the Brazilian savanna (Cerrado), little is known about the biogeographic origin of species occupying an extreme South American habitat type, the coastal dunes (Restinga). We examined the formation of this poorly known, endangered habitat by reconstructing the evolutionary history of two endemic species. Due to the proposed recency and uniqueness of this habitat, we hypothesized that Restinga species of the palm genus Allagoptera to be recently evolved and to present derived morphological characters. To detect habit shifts in absolute time, we used one plastid and nine nuclear genes to reconstruct the phylogenetic and biogeographic history of Allagoptera. We used light microscopy and stable isotope analysis to explore whether morphological adaptations occurred concomitantly with habitat shifts. Phylogenetic relationships were well supported and we found ancestral lineages of Allagoptera to be widely distributed throughout habitats that are currently occupied by extant species. Over the last ca. 7Ma Allagoptera has shifted its preference to increasingly dry habitats. Coincident with the colonization of the Cerrado and Restinga, morphological adaptations also evolved, including subterranean stems that are fire-resistant and long underground stem and root systems that facilitate water access. We did not find differences in metabolic pathway or modifications to pollen morphology when compared to other palm lineages. Assuming that the evolutionary history of Allagoptera is indicative of the habitat in which it occurs, our results infer a recent origin for Cerrado species. Although little is known about the formation of the Restinga habitat, our results also suggest a longer history than currently proposed; with an origin of Restinga habitats dating back to the Late Pliocene.

Biogeography of the Malagasy Celastraceae: Multiple independent origins followed by widespread dispersal of genera from Madagascar.

Of the 97 currently recognized genera of Celastraceae, 19 are native to Madagascar, including six endemics. In this study we conducted the most thorough phylogenetic analysis of Celastraceae yet completed with respect to both character and taxon sampling, and include representatives of five new endemic genera. Fifty-one new accessions, together with 328 previously used accessions of Celastrales, were sampled for morphological characters, two rDNA gene regions, and two plastid gene regions. The endemic Malagasy genera are resolved in two separate lineages-Xenodrys by itself and all other endemic genera in a clade that also includes four lineages inferred to have dispersed from Madagascar: Brexia madagascariensis (Mascarene Islands, coastal Africa), Elaeodendron (West Indies, Africa to New Caledonia), and Pleurostylia (Africa to New Caledonia). Of the 12 extant Malagasy Celastraceae lineages identified, eight are clearly of African origin. The origins of the remaining four lineages are less clear, but reasonable possibilities include America, Eurasia, Africa, southern India, Malesia, and Australia. Based on 95% credible age intervals from fossil-calibrated molecular dating, all 12 extant Malagasy Celastraceae lineages appear to have arisen following dispersal after the separation of Madagascar from other landmasses within the last 70 million years.

Revisiting the origin and diversification of vascular plants through a comprehensive Bayesian analysis of the fossil record.

Plants have a long evolutionary history, during which mass extinction events dramatically affected Earth's ecosystems and its biodiversity. The fossil record can shed light on the diversification dynamics of plant life and reveal how changes in the origination-extinction balance have contributed to shaping the current flora. We use a novel Bayesian approach to estimate origination and extinction rates in plants throughout their history. We focus on the effect of the 'Big Five' mass extinctions and on estimating the timing of origin of vascular plants, seed plants and angiosperms. Our analyses show that plant diversification is characterized by several shifts in origination and extinction rates, often matching the most important geological boundaries. The estimated origin of major plant clades predates the oldest macrofossils when considering the uncertainties associated with the fossil record and the preservation process. Our findings show that the commonly recognized mass extinctions have affected each plant group differently and that phases of high extinction often coincided with major floral turnovers. For instance, after the Cretaceous-Paleogene boundary we infer negligible shifts in diversification of nonflowering seed plants, but find significantly decreased extinction in spore-bearing plants and increased origination rates in angiosperms, contributing to their current ecological and evolutionary dominance.

Miocene dispersal drives island radiations in the palm tribe Trachycarpeae (Arecaceae).

The study of three island groups of the palm tribe Trachycarpeae (Arecaceae/Palmae) permits both the analysis of each independent radiation and comparisons across the tribe to address general processes that drive island diversification. Phylogenetic relationships of Trachycarpeae were inferred from three plastid and three low-copy nuclear genes. The incongruent topological position of Brahea in CISP5 was hypothesized to be caused by a gene duplication event and was addressed using uninode coding. The resulting phylogenetic trees were well-resolved and the genera were all highly supported except for Johannesteijsmannia and Serenoa. Divergence time analysis estimated the stem of the tribe to be approximately 86 Ma and the crown to be 38 Ma, indicating that significant extinction may have occurred along this branch. Historical biogeographic analysis suggested that Trachycarpeae are of southern North American, Central American, or Caribbean origin and supports previous hypotheses of a Laurasian origin. The biogeography and disjunctions within the tribe were interpreted with respect to divergence times, the fossil record, and geological factors such as the formation of the Greater Antilles--Aves Ridge, the Bering and the North Atlantic land bridges, tectonic movement in Southeast Asia, climatic shifts between the Eocene and Pliocene, and volcanism in the Pacific basin. In considering the three major island radiations within Trachycarpeae, Miocene dispersal appears to have been the driving force in allopatric speciation and is highlighted here as an emerging pattern across the tree of life.

Evaluating multiple criteria for species delimitation: an empirical example using Hawaiian palms (Arecaceae: Pritchardia).

BACKGROUND: Robust species delimitations are fundamental for conservation, evolutionary, and systematic studies, but they can be difficult to estimate, particularly in rapid and recent radiations. The consensus that species concepts aim to identify evolutionarily distinct lineages is clear, but the criteria used to distinguish evolutionary lineages differ based on the perceived importance of the various characteristics of evolving populations. We examined three different species-delimitation criteria (monophyly, absence of genetic intermediates, and diagnosability) to determine whether currently recognized species of Hawaiian Pritchardia are distinct lineages. RESULTS: Data from plastid and nuclear genes, microsatellite loci, and morphological characters resulted in various levels of lineage subdivision that were likely caused by differing evolutionary rates between data sources. Additionally, taxonomic entities may be confounded because of the effects of incomplete lineage sorting and/or gene flow. A coalescent species tree was largely congruent with the simultaneous analysis, consistent with the idea that incomplete lineage sorting did not mislead our results. Furthermore, gene flow among populations of sympatric lineages likely explains the admixture and lack of resolution between those groups. CONCLUSIONS: Delimiting Hawaiian Pritchardia species remains difficult but the ability to understand the influence of the evolutionary processes of incomplete lineage sorting and hybridization allow for mechanisms driving species diversity to be inferred. These processes likely extend to speciation in other Hawaiian angiosperm groups and the biota in general and must be explicitly accounted for in species delimitation.