Phylogenetic analyses of a combined data set suggest that the Attheya lineage is the closest living relative of the pennate diatoms (Bacillariophyceae)

A Bayesian analysis of a seven gene data set was conducted to reconstruct phylogenetic relationships among a sample of centric and pennate diatoms and to test alternative hypotheses about the closest living relative of Bacillariophyceae. A lineage, composed of two Attheya species, was inferred to share the most recent common ancestor with Bacillariophyceae--a relationship that was also corroborated by the combined parsimony analysis. All competing hypotheses about the closest living relative of Bacillariophyceae were rejected because 100% of the trees in the post-burn-in sample in the Bayesian analysis supported the Attheya-Bacillariophyceae clade. According to a partitioned Bremer support analysis, the majority of the genes in the combined data matrix supported the Attheya--Bacillariophyceae relationship. The global topology of the phylogenetic tree indicated that a monophyletic group consisting of Thalassiosirales and Toxarium undulatum formed the deepest branch followed by a node uniting a clade composed of Bacillariophyceae/Attheya species and a lineage made up of Eucampia zoodiacus, Chaetocerotales, Lithodesmiales, Triceratiales, Biddulphiales and Cymatosirales. Except for the phylogenetic positions of Lithodesmiales, Thalassiosira sp and Skeletonema costatum, the optimal tree obtained from the combined parsimony analysis showed the same branching order of taxa as those seen in the consensus tree inferred from three independent Markov chain Monte Carlo analyses. Noteworthy findings are that Toxarium undulatum shares a strongly supported node with Thalassiosirales and that the genus Attheya is not a member of the Chaetocerotales lineage.     

Mitochondrial DNA,morphology, and the phylogenetic relationships of Antarctic icefishes (Notothenioidei: Channichthyidae)

The Channichthyidae is a lineage of 16 species in the Notothenioidei, a clade of fishes that dominate Antarctic near-shore marine ecosystems with respect to both diversity and biomass. Among four published studies investigating channichthyid phylogeny, no two have produced the same tree topology, and no published study has investigated the degree of phylogenetic incongruence between existing molecular and morphological datasets. In this investigation we present an analysis of channichthyid phylogeny using complete gene sequences from two mitochondrial genes (ND2 and 16S) sampled from all recognized species in the clade. In addition, we have scored all 58 unique morphological characters used in three previous analyses of channichthyid phylogenetic relationships. Data partitions were analyzed separately to assess the amount of phylogenetic resolution provided by each dataset, and phylogenetic incongruence among data partitions was investigated using incongruence length difference (ILD) tests. We utilized a parsimony-based version of the Shimodaira-Hasegawa test to determine if alternative tree topologies are significantly different from trees resulting from maximum parsimony analysis of the combined partition dataset. Our results demonstrate that the greatest phylogenetic resolution is achieved when all molecular and morphological data partitions are combined into a single maximum parsimony analysis. Also, marginal to insignificant incongruence was detected among data partitions using the ILD. Maximum parsimony analysis of all data partitions combined results in a single tree, and is a unique hypothesis of phylogenetic relationships in the Channichthyidae. In particular, this hypothesis resolves the phylogenetic relationships of at least two species (Channichthys rhinoceratus and Chaenocephalus aceratus), for which there was no consensus among the previous phylogenetic hypotheses. The combined data partition dataset provides substantial statistical power to discriminate among alternative hypotheses of channichthyid relationships. These findings suggest the optimal strategy for investigating the phylogenetic relationships of channichthyids is one that uses all available phylogenetic data in analyses of combined data partitions.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

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.     

Molecular phylogeny of songbirds (Aves: Passeriformes) and the relative utility of common nuclear marker loci

While the monophyly of the largest avian order Passeriformes as well as its suborders suboscines (Tyranni) and oscines (Passeri) is well established, lower phylogenetic relationships of this fast radiated taxon have been a continuous matter of debate, especially within the suborder oscines. Many studies analyzing phylogenetic relationships of the Passeriformes using molecular markers have been published, which led to a better resolved phylogeny. Conflicting hypotheses and still remaining uncertainties, especially within the Passerida, have repeatedly stimulated further research with additional new markers. In the present study we used a combination of established molecular markers (RAG‐1, RAG‐2, c‐myc) and the recently introduced ZENK. We accomplished phylogenetic analyses using maximum parsimony, maximum likelihood and Bayesian inference, both separately for all genes and simultaneously. To assess the phylogenetic utility of the different genes in avian systematics we analyzed the influence of each data partition on the phylogenetic tree yielded by the combined approach using partitioned Bremer support. Compared with the other single gene analyses, the ZENK trees exhibited by far the best resolution and showed the lowest amount of homoplasy. Our data indicate that this gene is—at least in passerines—suitable for inference of even old taxonomic splits. Our combined analysis yields well‐supported phylogenetic hypotheses for passerine phylogeny and apart from corroborating recently proposed hypotheses on phylogenetic relationships in the Passeriformes we provide evidence for some new hypotheses. The subdivision of the Passerida into three superfamilies, Sylvioidea, Passeroidea and Muscicapoidea, the first as sister to the two latter groups is strongly supported. We found evidence for a split between Paridae and the remaining Sylvioidea. © The Willi Hennig Society 2007.     

Attributes and congruence of three molecular data sets: inferring phylogenies among Septoria-related species from woody perennial plants

To improve our understanding of phylogenetic relationships within the anamorphic genus Septoria, three molecular data sets representing 2,417 bp of nuclear and mitochondrial genes were evaluated. Separate gene analyses and combined analyses were performed using first, the maximum parsimony criterion and second, a Bayesian framework. The homogeneity of data partitions was evaluated via a combination of homogeneity partition tests and tree topology incongruence tests before conducting combined analyses. A last incongruence re-evaluation using partitioned Bremer support was performed on the combined tree, which corroborated the previous estimates. After each separate data set attributes were examined, simple explanations were advocated as the causes of the significant incongruences detected. The analysis of multiple gene partitions showed unprecedented phylogenetic resolution within the genus Septoria that supported the results from previously published single gene phylogenies. Specifically, we have delimited distinct but closely related species representing monophyletic groups that frequently correlated with their respective host families. Conversely, the occurrence of well-supported groups including closely related but distinct molecular taxa sampled on unrelated host-plants allowed us to reject, in these particular cases, the co-evolutionary concept expected between a parasite and its host and to discuss alternative evolutionary models recently proposed for these pathogens.     

Phylogeny of sipunculan worms: A combined analysis of four gene regions and morphology

The intra-phyletic relationships of sipunculan worms were analyzed based on DNA sequence data from four gene regions and 58 morphological characters. Initially we analyzed the data under direct optimization using parsimony as optimality criterion. An implied alignment resulting from the direct optimization analysis was subsequently utilized to perform a Bayesian analysis with mixed models for the different data partitions. For this we applied a doublet model for the stem regions of the 18S rRNA. Both analyses support monophyly of Sipuncula and most of the same clades within the phylum. The analyses differ with respect to the relationships among the major groups but whereas the deep nodes in the direct optimization analysis generally show low jackknife support, they are supported by 100% posterior probability in the Bayesian analysis. Direct optimization has been useful for handling sequences of unequal length and generating conservative phylogenetic hypotheses whereas the Bayesian analysis under mixed models provided high resolution in the basal nodes of the tree.     

Effects of data incompleteness on the relative performance of parsimony and Bayesian approaches in a supermatrix phylogenetic reconstruction of Mustelidae and Procyonidae (Carnivora)

Missing data are commonly thought to impede a resolved or accurate reconstruction of phylogenetic relationships, and probabilistic analysis techniques are increasingly viewed as less vulnerable to the negative effects of data incompleteness than parsimony analyses. We test both assumptions empirically by conducting parsimony and Bayesian analyses on an approximately 1.5 × 10⁶‐cell (27 965 characters × 52 species) mustelid–procyonid molecular supermatrix with 62.7% missing entries. Contrary to the first assumption, phylogenetic relationships inferred from our analyses are fully (Bayesian) or almost fully (parsimony) resolved topologically with mostly strong support and also largely in accord with prior molecular estimations of mustelid and procyonid phylogeny derived with parsimony, Bayesian, and other probabilistic analysis techniques from smaller but complete or nearly complete data sets. Contrary to the second assumption, we found no compelling evidence in support of a relationship between the inferior performance of parsimony and taxon incompleteness (i.e. the proportion of missing character data for a taxon), although we found evidence for a connection between the inferior performance of parsimony and character incompleteness (i.e. no overlap in character data between some taxa). The relatively good performance of our analyses may be related to the large number of sampled characters, so that most taxa (even highly incomplete ones) are represented by a sufficient number of characters allowing both approaches to resolve their relationships. © The Willi Hennig Society 2009.     

The phylogeny of Iberian Aphodiini species (Coleoptera, Scarabaeoidea, Scarabaeidae, Aphodiinae) based on morphology

Abstract.  The phylogeny of Iberian Aphodiini species was reconstructed based on morphology. Wing venation, mouthparts, male and female genitalia, and external morphology provided ninety-four characters scored for ninety-three Aphodiini species. Phylogenetic analyses were based on maximum parsimony and Bayesian inference criteria. Maximum parsimony consensus trees recovered Acrossus species as a sister group of the remaining Aphodiini, followed by two other branches, one including Neagolius , Plagiogonus , Ahermodontus and Ammoecius species, and the other including Oxyomus , Nimbus , Heptaulacus and Euheptaulacus species. The remaining studied taxa clustered in an unresolved group. Bayesian inference trees recovered Acrossus as the sister group of the remaining Iberian Aphodiini, followed by Colobopterus erraticus and the rest of the Iberian Aphodiini, but this latter branch was unresolved. The general lack of statistical support for the inferred phylogenetic relationships at terminal nodes using both maximum parsimony and Bayesian inference suggests that variation in morphological characters useful for phylogenetic inference in the present study is small, perhaps as a consequence of a radiation event occurring at the origin of the tribe. A probable evolutionary pattern for Aphodiini is proposed which infers six groups, namely Acrossian, Ammoecian, Oxyomian, Aphodian s.str., Colobopteran and Aphodian s.l. clades.     

Phylogenetic relationships of Steinernema Travassos, 1927 (Nematoda: Cephalobina: Steinernematidae) based on nuclear, mitochondrial and morphological data

Entomopathogenic nematodes of the genus Steinernema are lethal parasites of insects that are used as biological control agents of several lepidopteran, dipteran and coleopteran pests. Phylogenetic relationships among 25 Steinernema species were estimated using nucleotide sequences from three genes and 22 morphological characters. Parsimony analysis of 28S (LSU) sequences yielded a well-resolved phylogenetic hypothesis with reliable bootstrap support for 13 clades. Parsimony analysis of mitochondrial DNA sequences (12S rDNA and cox 1 genes) yielded phylogenetic trees with a lower consistency index than for LSU sequences, and with fewer reliably supported clades. Combined phylogenetic analysis of the 3-gene dataset by parsimony and Bayesian methods yielded well-resolved and highly similar trees. Bayesian posterior probabilities were high for most clades; bootstrap (parsimony) support was reliable for approximately half of the internal nodes. Parsimony analysis of the morphological dataset yielded a poorly resolved tree, whereas total evidence analysis (molecular plus morphological data) yielded a phylogenetic hypothesis consistent with, but less resolved than trees inferred from combined molecular data. Parsimony mapping of morphological characters on the 3-gene trees showed that most structural features of steinernematids are highly homoplastic. The distribution of nematode foraging strategies on these trees predicts that S. hermaphroditum, S. diaprepesi and S. longicaudum (US isolate) have cruise forager behaviours.     

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.     

Molecular systematics of the Eastern Fence Lizard (Sceloporus undulatus): a comparison of Parsimony, Likelihood, and Bayesian approaches

Phylogenetic analysis of large datasets using complex nucleotide substitution models under a maximum likelihood framework can be computationally infeasible, especially when attempting to infer confidence values by way of nonparametric bootstrapping. Recent developments in phylogenetics suggest the computational burden can be reduced by using Bayesian methods of phylogenetic inference. However, few empirical phylogenetic studies exist that explore the efficiency of Bayesian analysis of large datasets. To this end, we conducted an extensive phylogenetic analysis of the wide-ranging and geographically variable Eastern Fence Lizard (Sceloporus undulatus). Maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses were performed on a combined mitochondrial DNA dataset (12S and 16S rRNA, ND1 protein-coding gene, and associated tRNA; 3,688 bp total) for 56 populations of S. undulatus (78 total terminals including other S. undulatus group species and outgroups). Maximum parsimony analysis resulted in numerous equally parsimonious trees (82,646 from equally weighted parsimony and 335 from weighted parsimony). The majority rule consensus tree derived from the Bayesian analysis was topologically identical to the single best phylogeny inferred from the maximum likelihood analysis, but required approximately 80% less computational time. The mtDNA data provide strong support for the monophyly of the S. undulatus group and the paraphyly of "S. undulatus" with respect to S. belli, S. cautus, and S. woodi. Parallel evolution of ecomorphs within "S. undulatus" has masked the actual number of species within this group. This evidence, along with convincing patterns of phylogeographic differentiation suggests "S. undulatus" represents at least four lineages that should be recognized as evolutionary species.     

Bayesian mixed models and the phylogeny of pitvipers (Viperidae: Serpentes)

The subfamily Crotalinae (pitvipers) contains over 190 species of venomous snakes distributed in both the Old and New World. We incorporated an extensive sampling of taxa (including 28 of 29 genera), and sequences of four mitochondrial gene fragments (2.3kb) per individual, to estimate the phylogeny of pitvipers based on maximum parsimony and Bayesian phylogenetic methods. Our Bayesian analyses incorporated complex mixed models of nucleotide evolution that allocated independent models to various partitions of the dataset within combined analyses. We compared results of unpartitioned versus partitioned Bayesian analyses to investigate how much unpartitioned (versus partitioned) models were forced to compromise estimates of model parameters, and whether complex models substantially alter phylogenetic conclusions to the extent that they appear to extract more phylogenetic signal than simple models. Our results indicate that complex models do extract more phylogenetic signal from the data. We also address how differences in phylogenetic results (e.g., bipartition posterior probabilities) obtained from simple versus complex models may be interpreted in terms of relative credibility. Our estimates of pitviper phylogeny suggest that nearly all recently proposed generic reallocations appear valid, although certain Old and New World genera (Ovophis, Trimeresurus, and Bothrops) remain poly- or paraphyletic and require further taxonomic revision. While a majority of nodes were resolved, we could not confidently estimate the basal relationships among New World genera and which lineage of Old World species is most closely related to this New World group.     

Evolutionary relationships of the cup-fungus genus Peziza and Pezizaceae inferred from multiple nuclear genes: RPB2, beta-tubulin, and LSU rDNA

To provide a robust phylogeny of Pezizaceae, partial sequences from two nuclear protein-coding genes, RPB2 (encoding the second largest subunit of RNA polymerase II) and beta-tubulin, were obtained from 69 and 72 specimens, respectively, to analyze with nuclear ribosomal large subunit RNA gene sequences (LSU). The three-gene data set includes 32 species of Peziza, and 27 species from nine additional epigeous and six hypogeous (truffle) pezizaceous genera. Analyses of the combined LSU, RPB2, and beta-tubulin data set using parsimony, maximum likelihood, and Bayesian approaches identify 14 fine-scale lineages within Pezizaceae. Species of Peziza occur in eight of the lineages, spread among other genera of the family, confirming the non-monophyly of the genus. Although parsimony analyses of the three-gene data set produced a nearly completely resolved strict consensus tree, with increased confidence, relationships between the lineages are still resolved with mostly weak bootstrap support. Bayesian analyses of the three-gene data, however, show support for several more inclusive clades, mostly congruent with Bayesian analyses of RPB2. No strongly supported incongruence was found among phylogenies derived from the separate LSU, RPB2, and beta-tubulin data sets. The RPB2 region appeared to be the most informative single gene region based on resolution and clade support, and accounts for the greatest number of potentially parsimony informative characters within the combined data set, followed by the LSU and the beta-tubulin region. The results indicate that third codon positions in beta-tubulin are saturated, especially for sites that provide information about the deeper relationships. Nevertheless, almost all phylogenetic signal in beta-tubulin is due to third positions changes, with almost no signal in first and second codons, and contribute phylogenetic information at the "fine-scale" level within the Pezizaceae. The Pezizaceae is supported as monophyletic in analyses of the three-gene data set, but its sister-group relationships is not resolved with support. The results advocate the use of RPB2 as a marker for ascomycete phylogenetics at the inter-generic level, whereas the beta-tubulin gene appears less useful.     

Plethodontid salamander mitochondrial genomics: A parsimony evaluation of character conflict and implications for historical biogeography

A new parsimony analysis of 27 complete mitochondrial genomic sequences is conducted to investigate the phylogenetic relationships of plethodontid salamanders. This analysis focuses on the amount of character conflict between phylogenetic trees recovered from newly conducted parsimony searches and the Bayesian and maximum likelihood topology reported by Mueller et al. (2004 ; PNAS, 101, 13820–13825). Strong support for Hemidactylium as the sister taxon to all other plethodontids is recovered from parsimony analyses. Plotting area relationships on the most parsimonious phylogenetic tree suggests that eastern North America is the origin of the family Plethodontidae supporting the “Out of Appalachia” hypothesis. A new taxonomy that recognizes clades recovered from phylogenetic analyses is proposed. © The Willi Hennig Society 2005.     

Molecular phylogenetics and evolution of the endemic Hawaiian genus Adenophorus (Grammitidaceae)

Recent studies of the phylogeny of several groups of native Hawaiian vascular plants have led to significant insights into the origin and evolution of important elements of the Hawaiian flora. No groups of Hawaiian pteridophytes have been subjected previously to rigorous phylogenetic analysis. We conducted a molecular phylogenetic analysis of the endemic Hawaiian fern genus Adenophorus employing DNA sequence variation from three cpDNA fragments: rbcL, atpbeta, and the trnL-trnF intergenic spacer (IGS). In the phylogenetic analyses we employed maximum parsimony and Bayesian inference. Bayesian phylogenetic inference often provided stronger support for hypothetical relationships than did nonparametric bootstrap analyses. Although phylogenetic analyses of individual DNA fragments resulted in different patterns of relationships among species and varying levels of support for various clades, a combined analysis of all three sets of sequences produced one, strongly supported phylogenetic hypothesis. The primary features of that hypothesis are: (1) Adenophorus is monophyletic; (2) subgenus Oligadenus is paraphyletic; (3) the enigmatic endemic Hawaiian species Grammitis tenella is strongly supported as the sister taxon to Adenophorus; (4) highly divided leaf blades are evolutionarily derived in the group and simple leaves are ancestral; and, (5) the biogeographical origin of the common ancestor of the Adenophorus-G. tenella clade remains unresolved, although a neotropical origin seems most likely.     

Molecular systematics of the Jacks (Perciformes: Carangidae) based on mitochondrial cytochrome b sequences using parsimony,likelihood, and Bayesian approaches

The Carangidae represent a diverse family of marine fishes that include both ecologically and economically important species. Currently, there are four recognized tribes within the family, but phylogenetic relationships among them based on morphology are not resolved. In addition, the tribe Carangini contains species with a variety of body forms and no study has tried to interpret the evolution of this diversity. We used DNA sequences from the mitochondrial cytochrome b gene to reconstruct the phylogenetic history of 50 species from each of the four tribes of Carangidae and four carangoid outgroup taxa. We found support for the monophyly of three tribes within the Carangidae (Carangini, Naucratini, and Trachinotini); however, monophyly of the fourth tribe (Scomberoidini) remains questionable. A sister group relationship between the Carangini and the Naucratini is well supported. This clade is apparently sister to the Trachinotini plus Scomberoidini but there is uncertain support for this relationship. Additionally, we examined the evolution of body form within the tribe Carangini and determined that each of the predominant clades has a distinct evolutionary trend in body form. We tested three methods of phylogenetic inference, parsimony, maximum-likelihood, and Bayesian inference. Whereas the three analyses produced largely congruent hypotheses, they differed in several important relationships. Maximum-likelihood and Bayesian methods produced hypotheses with higher support values for deep branches. The Bayesian analysis was computationally much faster and yet produced phylogenetic hypotheses that were very similar to those of the maximum-likelihood analysis.     

Complete mitochondrial genome of Shinisaurus crocodilurus (Squamata: Shinisaurus) and its genetic relationship with related species

The 16,585 base pairs mitochondrial genome of Shinisaurus crocodilurus was determined by using PCR amplification and DNA sequencing. To determine the phylogenetic position of S. crocodilurus with related species within Squamata, the phylogenetic tree was reconstructed with the concatenated nucleotide sequences of the 12 heavy-strand-encoded protein genes. Phylogenetic analyses based on maximum parsimony and Bayesian inference methods consistently support that the S. crocodilurus was closely related to the Helodermatidae within a monophyletic Anguimorpha group. And the result here contradicted the monophyly of Varanoidea (Varanidae + Helodermatidae). In addition, the Gekkonidae was found to possess a basal phylogenetic position within squamata and the traditional hypothesis of monophyletic lineages of Iguania and Scleroglossa was not supported in this study.     

Phylogeny of shorebirds, gulls, and alcids (Aves: Charadrii) from the cytochrome-b gene: parsimony, Bayesian inference, minimum evolution, and quartet puzzling

Charadrii (shorebirds, gulls, and alcids) have exceptional diversity in ecological, behavioral, and life-history traits. A phylogenetic framework is necessary to fully understand the relationships among these traits. Despite several attempts to resolve the phylogeny of the Charadrii, none have comprehensively utilized molecular sequence data. Complete and partial cytochrome-b gene sequences for 86 Charadrii and five Falconides species (as outgroup taxa) were obtained from GenBank and aligned. We analyzed the resulting matrices using parsimony, Bayesian inference, minimum evolution, and quartet puzzling methods. Posterior probabilities, decay indices, and bootstrapping provide strong support for four major lineages consisting of gulls, alcids, plovers, and sandpipers, respectively. The broad structure of the trees differ significantly from all previous hypotheses of Charadrii phylogeny in placing the plovers at the base of the tree below the sandpipers in a pectinate sequence towards a large clade of gulls and alcids. The parsimony, Bayesian, and minimum evolution models provide strong evidence for this phylogenetic hypothesis. This is further corroborated by non-tree based measures of support and conflict (Lento plots). The quartet puzzling trees are poorly resolved and inconclusive.