Resolving an ancient, rapid radiation in Saxifragales

Despite the prior use of approximately 9000 bp, deep-level relationships within the angiosperm clade, Saxifragales remain enigmatic, due to an ancient, rapid radiation (89.5 to 110 Ma based on the fossil record). To resolve these deep relationships, we constructed several new data sets: (1) 16 genes representing the three genomic compartments within plant cells (2 nuclear, 10 plastid, 4 mitochondrial; aligned, analyzed length = 21,460 bp) for 28 taxa; (2) the entire plastid inverted repeat (IR; 26,625 bp) for 17 taxa; (3) "total evidence" (50,845 bp) for both 17 and 28 taxa (the latter missing the IR). Bayesian and ML methods yielded identical topologies across partitions with most clades receiving high posterior probability (pp = 1.0) and bootstrap (95% to 100%) values, suggesting that with sufficient data, rapid radiations can be resolved. In contrast, parsimony analyses of different partitions yielded conflicting topologies, particularly with respect to the placement of Paeoniaceae, a clade characterized by a long branch. In agreement with published simulations, the addition of characters increased bootstrap support for the putatively erroneous placement of Paeoniaceae. Although having far fewer parsimony-informative sites, slowly evolving plastid genes provided higher resolution and support for deep-level relationships than rapidly evolving plastid genes, yielding a topology close to the Bayesian and ML total evidence tree. The plastid IR region may be an ideal source of slowly evolving genes for resolution of deep-level angiosperm divergences that date to 90 My or more. Rapidly evolving genes provided support for tip relationships not recovered with slowly evolving genes, indicating some complementarity. Age estimates using penalized likelihood with and without age constraints for the 28-taxon, total evidence data set are comparable to fossil dates, whereas estimates based on the 17-taxon data are much older than implied by the fossil record. Hence, sufficient taxon density, and not simply numerous base pairs, is important in reliably estimating ages. Age estimates indicate that the early diversification of Saxifragales occurred rapidly, over a time span as short as 6 million years. Between 25,000 and 50,000 bp were needed to resolve this radiation with high support values. Extrapolating from Saxifragales, a similar number of base pairs may be needed to resolve the many other deep-level radiations of comparable age in angiosperms.     

Sequencing of whole plastid genomes and nuclear ribosomal DNA of Diospyros species (Ebenaceae) endemic to New Caledonia: many species,little divergence

Background and Aims Some plant groups, especially on islands, have been shaped by strong ancestral bottlenecks and rapid, recent radiation of phenotypic characters. Single molecular markers are often not informative enough for phylogenetic reconstruction in such plant groups. Whole plastid genomes and nuclear ribosomal DNA (nrDNA) are viewed by many researchers as sources of information for phylogenetic reconstruction of groups in which expected levels of divergence in standard markers are low. Here we evaluate the usefulness of these data types to resolve phylogenetic relationships among closely related Diospyros species.Methods Twenty-two closely related Diospyros species from New Caledonia were investigated using whole plastid genomes and nrDNA data from low-coverage next-generation sequencing (NGS). Phylogenetic trees were inferred using maximum parsimony, maximum likelihood and Bayesian inference on separate plastid and nrDNA and combined matrices.Key Results The plastid and nrDNA sequences were, singly and together, unable to provide well supported phylogenetic relationships among the closely related New Caledonian Diospyros species. In the nrDNA, a 6-fold greater percentage of parsimony-informative characters compared with plastid DNA was found, but the total number of informative sites was greater for the much larger plastid DNA genomes. Combining the plastid and nuclear data improved resolution. Plastid results showed a trend towards geographical clustering of accessions rather than following taxonomic species.Conclusions In plant groups in which multiple plastid markers are not sufficiently informative, an investigation at the level of the entire plastid genome may also not be sufficient for detailed phylogenetic reconstruction. Sequencing of complete plastid genomes and nrDNA repeats seems to clarify some relationships among the New Caledonian Diospyros species, but the higher percentage of parsimony-informative characters in nrDNA compared with plastid DNA did not help to resolve the phylogenetic tree because the total number of variable sites was much lower than in the entire plastid genome. The geographical clustering of the individuals against a background of overall low sequence divergence could indicate transfer of plastid genomes due to hybridization and introgression following secondary contact.     

Phylogenetic signal and noise: predicting the power of a data set to resolve phylogeny

A principal objective for phylogenetic experimental design is to predict the power of a data set to resolve nodes in a phylogenetic tree. However, proactively assessing the potential for phylogenetic noise compared with signal in a candidate data set has been a formidable challenge. Understanding the impact of collection of additional sequence data to resolve recalcitrant internodes at diverse historical times will facilitate increasingly accurate and cost-effective phylogenetic research. Here, we derive theory based on the fundamental unit of the phylogenetic tree, the quartet, that applies estimates of the state space and the rates of evolution of characters in a data set to predict phylogenetic signal and phylogenetic noise and therefore to predict the power to resolve internodes. We develop and implement a Monte Carlo approach to estimating power to resolve as well as deriving a nearly equivalent faster deterministic calculation. These approaches are applied to describe the distribution of potential signal, polytomy, or noise for two example data sets, one recent (cytochrome c oxidase I and 28S ribosomal rRNA sequences from Diplazontinae parasitoid wasps) and one deep (eight nuclear genes and a phylogenomic sequence for diverse microbial eukaryotes including Stramenopiles, Alveolata, and Rhizaria). The predicted power of resolution for the loci analyzed is consistent with the historic use of the genes in phylogenetics.     

Phylogeny of Saxifragales (angiosperms,eudicots): analysis of a rapid,ancient radiation

Rapid, ancient radiations pose one of the most difficult challenges for phylogenetic estimation. We used DNA sequence data of 9,006 aligned base pairs from five genes (chloroplast atpB, matK, rbcL, and 18S and 26S nrDNA) to elucidate relationships among major lineages of Saxifragales (angiosperms, eudicots). These relationships were poorly supported in previous studies, apparently because the lineages originated in rapid succession. Using an array of methods that explicitly incorporate assumptions about evolutionary process (weighted maximum parsimony, maximum likelihood, LogDet/paralinear transformed distances), we show that the initial diversification of Saxifragales was indeed rapid. We suggest that the poor resolution of our best phylogenetic estimate is not due to violations of assumptions or to combining data partitions having conflicting histories or processes. We show that estimated branch lengths during the initial diversification are exceedingly short, and we estimate that acquiring sufficient sequence data to resolve these relationships would require an extraordinary effort (approximately 10(7) bp), assuming a linear increase in branch support with branch length. However, our simulation of much larger data sets containing a distribution of phylogenetic signal similar to that of the five sampled gene sequences suggests a limit to achievable branch support. Using statistical tests of differences in the likelihoods of topologies, we evaluated whether the initial radiation of Saxifragales involved the simultaneous origin of major lineages. Our results are consistent with predictions that resolving the branching order of rapid, ancient radiations requires sampling characters that evolved rapidly at the time of the radiation but have since experienced a slower evolutionary rate.     

Phylogenomics of the hyperdiverse daisy tribes: Anthemideae,Astereae, Calenduleae,Gnaphalieae, and Senecioneae

Asteraceae account for 10% of all flowering plant species, and 35%–40% of these are in five closely related tribes that total over 10 000 species. These tribes include Anthemideae, Astereae, Calenduleae, Gnaphalieae, and Senecioneae, which form one of two enormous clades within Subfamily Asteroideae. We took a phylogenomics approach to resolve evolutionary relationships among these five tribes. We sampled the nuclear and plastid genomes via HybSeq target enrichment and genome skimming, and recovered 74 plastid genes and nearly 1000 nuclear loci, known as Conserved Orthologous Sequences. We tested for conflicting support in both data sets and used network analyses to assess patterns of reticulation to explain the early evolutionary history of this lineage, which has experienced whole‐genome duplications and rapid radiations. We found concordance and conflicting support in both data sets and documented four ancient hybridization events. Due to the timing of the early radiation of this five‐tribe lineage, shortly before the Eocene–Oligocene extinction event (34 MYA), early lineages were likely lost, obscuring some details of their early evolutionary history.     

Molecular phylogeny of Malagasy reed frogs, Heterixalus, and the relative performance of bioacoustics and color-patterns for resolving their systematics

The members of the genus Heterixalus constitute one of the endemic frog radiations in Madagascar. Here we present a complete species-level phylogeny based on DNA sequences (4876 base pairs) of three nuclear and four mitochondrial markers to clarify the phylogenetic relationships among and within all known species of this genus, as well as the phylogenetic position of the monospecific Seychellean Tachycnemis seychellensis. Although the performance to resolve supported clades of Heterixalus species differed among the investigated gene fragments when analyzed separately, we could identify five well-supported species groups within Heterixalus in the combined analysis of all gene fragments. Our data strongly support a Heterixalus-Tachycnemis clade, and indicate the probable monophyly of Heterixalus placed sister to Tachycnemis. However, the diversification of these lineages may have happened in a short interval of time, leading to an unstable placement of Tachycnemis in the single-gene fragment phylogenies. Referring to the hitherto existing classification of Heterixalus, which is predominantly based on chromatic and bioacoustic characters, we examined the relative performance of these data sets relative to our molecular phylogeny. A Bayesian tree reconstructed with a bioacoustic data set yielded a higher resemblance to the molecular phylogeny than a tree constructed using a chromatic data set, which supports the importance of bioacoustic characters for systematic analyses of these anurans.     

Conflicting phylogenetic signals in genomic data of the coffee family (Rubiaceae)

Reconstructions of phylogenetic relationships in the flowering plant family Rubiaceae have up until now relied heavily on single‐ or multi‐gene data, primarily from the plastid compartment. With the availability of cost‐ and time‐efficient techniques for generating complete genome sequences, the opportunity arises to resolve some of the relationships that, up until now, have proven problematic. Here, we contribute new data from complete 58 plastid genome sequences, representing 55 of the currently 65 recognized tribes of the Rubiaceae. Also contributed are new data from the nuclear rDNA cistrons for corresponding taxa. Phylogenetic analyses are conducted on two plastid data sets, one including data from the protein coding genes only, and a second where protein coding data are combined with non‐coding regions, and on a nuclear rDNA data set. Our results clearly show that simply adopting a “more characters” approach does not resolve the relationships in the Rubiaceae. More importantly, we identify conflicting phylogenetic signals in the data. Analyses of the same plastid data, treated as nucleotides or as codon‐degenerated data, resolve and support conflicting topologies in the subfamily Cinchonoideae. As these analyses use the same data, we interpret the conflict to result from erroneous assumptions in the models used to reconstruct our phylogenies. Conflicting signals are also identified in the analyses of the plastid versus the nuclear rDNA data sets. These analyses use data from different genomic compartments, with different inheritance patterns, and we interpret the conflicts as representing “real” conflicts, reflecting biological processes of the past.     

Phylogeny of the Liliales (Monocotyledons) with special emphasis on data partition congruence and RNA editing

A phylogenetic analysis of the monocot order Liliales was performed using sequence data from three mitochondrial (atp1, cob, nad5) and two plastid genes (rbcL, ndhF). The complete data matrix includes 46 terminals representing all 10 families currently included in Liliales. The two major partitions, mitochondrial and plastid data, were congruent, and parsimony analysis resulted in 50 equally parsimonious trees and a well resolved consensus tree confirming monophyly of all families. Mitochondrial genes are known to include RNA edited sites, and in some cases unprocessed genes are replaced by retro‐processed gene copies, that is processed paralogs. To test the effects on phylogeny reconstruction of predicted edited sites and potentially unintentionally sampled processed paralogs, a number of analyses were performed using subsets of the complete data matrix. In general, predicted edited sites were more homoplasious than the other characters and increased incongruence among most data partitions. The predicted edited sites have a non‐random phylogenetic signal in conflict with the signal of the non‐edited sites. The potentially misleading signal was caused partially by the apparent presence of processed paralogs in Galanthus (Amaryllidaceae), part of the outgroup, but also by a deviating evolutionary pattern of predicted edited sites in Liliaceae compared with the remainder of the Liliales. Despite the problems that processed paralogs may cause, we argue that they should not a priori be excluded from phylogenetic analysis.     

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.     

Resolving phylogenetic signal from noise when divergence is rapid: a new look at the old problem of echinoderm class relationships

Resolving evolutionary relationships in groups that underwent fast radiation in deep time is a problem for molecular phylogeny, as the scant phylogenetic signal that characterises short internal branches is generally swamped by more recent substitutions. We implement an approach, that maps how the support for rival phylogenies changes when analysing subsets of sites with either faster and more heterogeneous rates or slower and more homogeneous rates, to address a long-standing problem in deuterostome phylogeny - the interrelationships of the eleutherozoan echinoderm classes. We show that miRNA genes are phylogenetically uninformative as to the relationships of asteroids, echinoids and ophiuroids, consistent with a rapid radiation of these groups as suggested by their fossil record. Using three nuclear rRNAs and seven nuclear housekeeping genes, we map the support for the three possible phylogenetic arrangements of asteroids, ophiuroids and echinoids when moving between subsets of the data with very similar or very different rates of evolution. Only one of the three possible topologies (asteroids (ophiuroids + echinoids)) strengthens when the most rate-homogeneous subset of data are analysed. The other two possible pairings become stronger in a less reliable data subset, which includes the fastest and thus homoplasy-rich data in our alignment. Thus, while superficial analysis of our concatenated alignment identifies asteroids and ophiuroids as sister taxa, more thorough analyses suggest that ophiuroids may be more closely related to echinoids. Divergence of these echinoderm groups, using a relaxed molecular clock, is estimated to have occurred within ∼5 million years. Our results illustrate that the analytic approach of phylogenetic signal dissection can be a powerful tool to investigate rapid radiations in deep geologic time.     

The Case for a Single‐Plastid Origin Revisited: Convergent Evolution in Organellar Gene Content

Three eukaryotic lineages generally are believed to have plastids that are primary in origin; that is, descended directly from a cyanobacterial endosymbiont. The recovery of these plastids as a monophyletic group in most molecular phylogenetic analyses, along with similarities in genome content and protein targeting mechanisms, have been cited as strong evidence in support of the hypothesis of a single endosymbiotic origin of all plastids. Although these data indeed are consistent with a single plastid origin, they also are consistent with the proposition of multiple endosymbiotic origins. Each hypothesis requires certain evolutionary assumptions in order to be reconciled with all existing data; at present, it is unclear which of these assumptions most likely reflect the historical process that gave rise to plastid diversity. Here we examine similarities in gene content among representatives of the three primary plastid lineages, using as a control the genome of a mitochondrion that almost certainly originated as an independent endosymbiotic association. To minimize metabolic constraints on gene retention we focus on two datasets, ribosomal protein and transfer RNA genes, neither of which is tied directly to specific organellar functions. Analyses of all possible pair‐wise comparisons among the three plastids and mitochondrion indicate that genomic similarities are most consistent with convergent evolution due to constraints on gene loss, rather than with hypothesized shared evolutionary histories. We find no evidence of phylogenetic signal in the pattern of gene loss overlying this convergence. In light of these results, we address other lines of evidence and arguments that have been raised in support of a single plastid origin.     

Profiling phylogenetic informativeness

The resolution of four controversial topics in phylogenetic experimental design hinges upon the informativeness of characters about the historical relationships among taxa. These controversies regard the power of different classes of phylogenetic character, the relative utility of increased taxonomic versus character sampling, the differentiation between lack of phylogenetic signal and a historical rapid radiation, and the design of taxonomically broad phylogenetic studies optimized by taxonomically sparse genome-scale data. Quantification of the informativeness of characters for resolution of phylogenetic hypotheses during specified historical epochs is key to the resolution of these controversies. Here, such a measure of phylogenetic informativeness is formulated. The optimal rate of evolution of a character to resolve a dated four-taxon polytomy is derived. By scaling the asymptotic informativeness of a character evolving at a nonoptimal rate by the derived asymptotic optimum, and by normalizing so that net phylogenetic informativeness is equivalent for all rates when integrated across all of history, an informativeness profile across history is derived. Calculation of the informativeness per base pair allows estimation of the cost-effectiveness of character sampling. Calculation of the informativeness per million years allows comparison across historical radiations of the utility of a gene for the inference of rapid adaptive radiation. The theory is applied to profile the phylogenetic informativeness of the genes BRCA1, RAG1, GHR, and c-myc from a muroid rodent sequence data set. Bounded integrations of the phylogenetic profile of these genes over four epochs comprising the diversifications of the muroid rodents, the mammals, the lobe-limbed vertebrates, and the early metazoans demonstrate the differential power of these genes to resolve the branching order among ancestral lineages. This measure of phylogenetic informativeness yields a new kind of information for evaluation of phylogenetic experiments. It conveys the utility of the addition of characters a phylogenetic study and it provides a basis for deciding whether appropriate phylogenetic power has been applied to a polytomy that is proposed to be a rapid radiation. Moreover, it provides a quantitative measure of the capacity of a gene to resolve soft polytomies.     

Examining monophyly in a large radiation of Madagascan butterflies (Lepidoptera: Satyrinae: Mycalesina) based on mitochondrial DNA data

The satyrine butterfly subtribe Mycalesina has undergone one of the more spectacular evolutionary radiations of butterflies in the Old World tropics. Perhaps the most phenotypically pronounced diversification of the group has occurred in the Malagasy region, where 68 currently recognized species are divided among five genera. Here, we report the results of phylogenetic analyses of sequence data from the cytochrome c oxidase II and cytochrome b mitochondrial genes, for a total of 54 mycalesine taxa, mostly from Madagascar. These molecular data complement an existing data set based on male morphological characters. The molecular results support the suggestion from morphology that three of the five Malagasy genera are paraphyletic and support the monophyly of at least three major morphological clades. Novel hypotheses of terminal taxon pairs are generated by the molecular data. Dense taxon sampling appears to be crucial for elucidating phylogenetic relationships within this large radiation. A potentially complex scenario for the origin of Malagasy mycalesines is proposed.     

The phylogenetic history of Selaginellaceae based on DNA sequences from the plastid and nucleus: extreme substitution rates and rate heterogeneity

Molecular phylogenetic research on Selaginellaceae has focused on the plastid gene rbcL, which in this family has unusually high substitution rates. Here we develop a molecular data set from the nuclear 26S ribosomal DNA gene with the aim of evaluating and extending the results of previous phylogenetic research. The 26S rDNA and the rbcL regions were sequenced for a sample of 23 species, which represent the main elements of species diversity in the family. The data were analysed independently and in combination using both maximum parsimony and Bayesian inference. Although several between genome differences were found, the general pattern of relationships uncovered by all analyses was very similar. Results corroborate the previous study supporting new groupings not previously recognised on morphological grounds. Substitution rates in the 26S rDNA were also found to be high (26% informative) for the region analysed, but lower than for rbcL (37% informative). These data indicate that high substitution rates might be widespread in all three genomes (i.e., plastid, mitochondrion, and nucleus).     

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.     

Phylogenetic signal and its decay in mitochondrial SSU and LSU rRNA gene fragments of Anisoptera

The phylogeny of Anisoptera, dragonflies in the strict sense, has proven to be notoriously difficult to resolve. Based on morphological characters, several recent publications dealing with the phylogeny of dragonflies proposed contradicting inter- and intrafamily relationships. We explored phylogenetic information content of mitochondrial large-subunit (LSU) and small-subunit (SSU) ribosomal gene fragments for these systematic problems. Starting at published universal primers, we developed primer sets suitable for amplifying large parts of the LSU and SSU rRNA genes within dragonflies. These fragments turned out to harbor sufficient phylogenetic information to satisfyingly resolve intrafamily relationships, but they contain insufficient phylogenetic structure to permit reliable conclusions about several interfamily relationships. We demonstrate that decay of phylogenetic signal progresses from intrafamily to interfamily to outgroup relationships and is correlated with an increase of genetic distances. As expected, signal decay is most pronounced in fast-changing sites. Additionally, base composition among fast-changing sites significantly deviates from the expected homogeneity. Homogeneity of base composition among all included taxa was restored only after removing fast-changing sites from the data set. The molecular data tentatively support interfamily relationships proposed by the most recent publication based on morphological characters of fossil and extant dragonflies.     

Phylogenomic analysis supports the monophyly of cryptophytes and haptophytes and the association of rhizaria with chromalveolates

Here we use phylogenomics with expressed sequence tag (EST) data from the ecologically important coccolithophore-forming alga Emiliania huxleyi and the plastid-lacking cryptophyte Goniomonas cf. pacifica to establish their phylogenetic positions in the eukaryotic tree. Haptophytes and cryptophytes are members of the putative eukaryotic supergroup Chromalveolata (chromists [cryptophytes, haptophytes, stramenopiles] and alveolates [apicomplexans, ciliates, and dinoflagellates]). The chromalveolates are postulated to be monophyletic on the basis of plastid pigmentation in photosynthetic members, plastid gene and genome relationships, nuclear "host" phylogenies of some chromalveolate lineages, unique gene duplication and replacements shared by these taxa, and the evolutionary history of components of the plastid import and translocation systems. However the phylogenetic position of cryptophytes and haptophytes and the monophyly of chromalveolates as a whole remain to be substantiated. Here we assess chromalveolate monophyly using a multigene dataset of nuclear genes that includes members of all 6 eukaryotic supergroups. An automated phylogenomics pipeline followed by targeted database searches was used to assemble a 16-protein dataset (6,735 aa) from 46 taxa for tree inference. Maximum likelihood and Bayesian analyses of these data support the monophyly of haptophytes and cryptophytes. This relationship is consistent with a gene replacement via horizontal gene transfer of plastid-encoded rpl36 that is uniquely shared by these taxa. The haptophytes + cryptophytes are sister to a clade that includes all other chromalveolates and, surprisingly, two members of the Rhizaria, Reticulomyxa filosa and Bigelowiella natans. The association of the two Rhizaria with chromalveolates is supported by the approximately unbiased (AU)-test and when the fastest evolving amino acid sites are removed from the 16-protein alignment.     

The enigmatic monotypic crab plover Dromas ardeola is closely related to pratincoles and coursers (Aves, Charadriiformes, Glareolidae)

The phylogenetic placement of the monotypic crab plover Dromasardeola (Aves, Charadriiformes) remains controversial. Phylogenetic analysis of anatomical and behavioral traits using phenetic and cladistic methods of tree inference have resulted in conflicting tree topologies, suggesting a close association of Dromas to members of different suborders and lineages within Charadriiformes. Here, we revisited the issue by applying Bayesian and parsimony methods of tree inference to 2,012 anatomical and 5,183 molecular characters to a set of 22 shorebird genera (including Turnix). Our results suggest that Bayesian analysis of anatomical characters does not resolve the phylogenetic relationship of shorebirds with strong statistical support. In contrast, Bayesian and parsimony tree inference from molecular data provided much stronger support for the phylogenetic relationships within shorebirds, and support a sister relationship of Dromas to Glareolidae (pratincoles and coursers), in agreement with previously published DNA-DNA hybridization studies.