Phylogenomic analysis suggests Coreidae and Alydidae (Hemiptera: Heteroptera) are not monophyletic

Next‐generation sequencing technologies (NGS) allow systematists to amass a wealth of genomic data from non‐model species for phylogenetic resolution at various temporal scales. However, phylogenetic inference for many lineages dominated by non‐model species has not yet benefited from NGS, which can complement Sanger sequencing studies. One such lineage, whose phylogenetic relationships remain uncertain, is the diverse, agriculturally important and charismatic Coreoidea (Hemiptera: Heteroptera). Given the lack of consensus on higher‐level relationships and the importance of a robust phylogeny for evolutionary hypothesis testing, we use a large data set comprised of hundreds of ultraconserved element (UCE) loci to infer the phylogeny of Coreoidea (excluding Stenocephalidae and Hyocephalidae), with emphasis on the families Coreidae and Alydidae. We generated three data sets by including alignments that contained loci sampled for at least 50%, 60%, or 70% of the total taxa, and inferred phylogeny using maximum likelihood and summary coalescent methods. Twenty‐six external morphological features used in relatively comprehensive phylogenetic analyses of coreoids were also re‐evaluated within our molecular phylogenetic framework. We recovered 439–970 loci per species (16%–36% of loci targeted) and combined this with previously generated UCE data for 12 taxa. All data sets, regardless of analytical approach, yielded topologically similar and strongly supported trees, with the exception of outgroup relationships and the position of Hydarinae. We recovered a monophyletic Coreoidea, with Rhopalidae highly supported as the sister group to Alydidae + Coreidae. Neither Alydidae nor Coreidae were monophyletic; the coreid subfamilies Hydarinae and Pseudophloeinae were recovered as more closely related to Alydidae than to other coreid subfamilies. Coreinae were paraphyletic with respect to Meropachyinae. Most morphological traits were homoplastic with several clades defined by few, if any, synapomorphies. Our results demonstrate the utility of phylogenomic approaches in generating robust hypotheses for taxa with long‐standing phylogenetic problems and highlight that novel insights may come from such approaches.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Assessing support for Blaberoidea phylogeny suggests optimal locus quality

Phylogenomics seeks to use next‐generation data to robustly infer an organism's evolutionary history. Yet, the practical caveats of phylogenomics motivate investigation of improved efficiency, particularly when quality of phylogenies are questionable. To achieve improvements, one goal is to maintain or enhance the quality of phylogenetic inference while severely reducing dataset size. We approach this by assessing which kinds of loci in phylogenomic alignments provide the majority of support for a phylogenetic inference of cockroaches in Blaberoidea. We examine locus substitution rate, saturation, evolutionary divergence, rate heterogeneity, stabilizing selection, and a priori information content as traits that may determine optimality. Our controlled experimental design is based on 265 loci for 102 blaberoidean taxa and 22 outgroup species. Loci with high substitution rate, low saturation, low sequence distance, low rate heterogeneity, and strong stabilizing selection derive more support for phylogenetic relationships. We found that some phylogenetic information content estimators may not be meaningful for assessing information content a priori. We use these findings to design concatenated datasets with an optimized subsample of 100 loci. The tree inferred from the optimized subsample alignment was largely identical to that inferred from all 265 loci but with less evidence of long branch attraction, improved statistical support, and potential 4‐6x improvements to computation time. Supported by phylogenetic and morphological evidence, we erect three newly named clades (Anallactinae Evangelista & Wipfler subfam. nov., Orkrasomeria tax. nov. Evangelista, Wipfler, & Béthoux and Hemithyrsocerini Evangelista tribe nov.) and propose other taxonomic modifications. The diagnosis of Pseudophyllodromiidae Grandcolas, 1996 is modified to accommodate Anallactinae and Pseudophyllodromiinae Vickery & Kevan, 1983. The diagnosis of Ectobiidae Brunner von Wattenwyl, 1865 is modified to add novel morphological characters.     

Chloroplast phylogenomics of unicellular and colonial Volvocales provides perspectives on the evolution of morphological characters

Volvocales forms a species-rich clade with wide morphological variety and is regarded as an ideal model for tracing the evolutionary transitions in multicellularity. The phylogenetic relationships among the colonial volvocine algae and its relatives are important for investigating the origin of multicellularity in the clade Reinhardtinia. Therefore, a robust phylogenetic framework of the unicellular and colonial volvocine algae with broad taxon and gene sampling is essential for illuminating the evolution of multicellularity. Recent chloroplast phylogenomic studies have uncovered five major orders in the Chlorophyceae, but the family-level relationships within Sphaeropleales and Volvocales remain elusive due to the uncertain positions of some incertae sedis taxa. In this study, we contributed six newly sequenced chloroplast genomes in the Volvocales and analyzed a dataset with 91 chlorophycean taxa and 58 protein-coding genes. Conflicting phylogenetic signals were detected among chloroplast genes that resulted in discordant tree topologies among different analyses. We compared the phylogenetic trees inferred from original nucleotide, RY-coding, codon-degenerate, and amino acid datasets, and improved the robustness of phylogenetic inference in the Chlorophyceae by reducing base compositional bias. Our analyses indicate that the unicellular Chlamydomonas and Vitreochlamys are close to or nested within the colonial taxa, and all the incertae sedis taxa are nested within the monophyletic Sphaeropleales s.l. We propose that the colonial taxa in the Reinhardtinia are paraphyletic and multicellularity evolved once in the volvocine green algae and might be lost in Chlamydomonas and Vitreochlamys.     

Gross brain morphology of Hypoptopomatinae and Neoplecostominae (Siluriformes: Loricariidae): Comparative anatomy and phylogenetic implications

Hypoptopomatinae and Neoplecostominae include about 250 valid species, a substantial portion of loricariid catfishes. Although the relationships among the members of these subfamilies have been inferred by many authors, the most recent hypotheses based on morphological and molecular data differ widely. Herein, we provide new data on the morphology of the central nervous system, and evaluate the usefulness of these characters in phylogenetic inference. To accomplish this, we characterized the gross brain morphology of those catfishes, and analyzed 54 neuroanatomical characters in a total of 40 terminal taxa representing Hypoptopomatinae and Neoplecostominae, and also members of Delturinae and Hypostominae as outgroups. Hypoptopomatinae and Neoplecostominae are recovered as separate subfamilies, and most of our results are compatible with morphology‐based analyses. We conclude that neuroanatomy provides an informative source of new characters with strong phylogenetic signal at all recovered taxonomic levels.     

A molecular phylogeny of historical and contemporary specimens of an under‐studied micro‐invertebrate group

Resolution of relationships at lower taxonomic levels is crucial for answering many evolutionary questions, and as such, sufficiently varied species representation is vital. This latter goal is not always achievable with relatively fresh samples. To alleviate the difficulties in procuring rarer taxa, we have seen increasing utilization of historical specimens in building molecular phylogenies using high throughput sequencing. This effort, however, has mainly focused on large‐bodied or well‐studied groups, with small‐bodied and under‐studied taxa under‐prioritized. Here, we utilize both historical and contemporary specimens, to increase the resolution of phylogenetic relationships among a group of under‐studied and small‐bodied metazoans, namely, cheilostome bryozoans. In this study, we pioneer the sequencing of air‐dried cheilostomes, utilizing a recently developed library preparation method for low DNA input. We evaluate a de novo mitogenome assembly and two iterative methods, using the sequenced target specimen as a reference for mapping, for our sequences. In doing so, we present mitochondrial and ribosomal RNA sequences of 43 cheilostomes representing 37 species, including 14 from historical samples ranging from 50 to 149 years old. The inferred phylogenetic relationships of these samples, analyzed together with publicly available sequence data, are shown in a statistically well‐supported 65 taxa and 17 genes cheilostome tree, which is also the most broadly sampled and largest to date. The robust phylogenetic placement of historical samples whose contemporary conspecifics and/or congenerics have been sequenced verifies the appropriateness of our workflow and gives confidence in the phylogenetic placement of those historical samples for which there are no close relatives sequenced. The success of our workflow is highlighted by the circularization of a total of 27 mitogenomes, seven from historical cheilostome samples. Our study highlights the potential of utilizing DNA from micro‐invertebrate specimens stored in natural history collections for resolving phylogenetic relationships among species.     

Is RAD‐seq suitable for phylogenetic inference? An in silico assessment and optimization

Inferring phylogenetic relationships between closely related taxa can be hindered by three factors: (1) the lack of informative molecular variation at short evolutionary timescale; (2) the lack of established markers in poorly studied taxa; and (3) the potential phylogenetic conflicts among different genomic regions due to incomplete lineage sorting or introgression. In this context, Restriction site Associated DNA sequencing (RAD‐seq) seems promising as this technique can generate sequence data from numerous DNA fragments scattered throughout the genome, from a large number of samples, and without preliminary knowledge on the taxa under study. However, divergence beyond the within‐species level will necessarily reduce the number of conserved and non‐duplicated restriction sites, and therefore the number of loci usable for phylogenetic inference. Here, we assess the suitability of RAD‐seq for phylogeny using a simulated experiment on the 12 Drosophila genomes, with divergence times ranging from 5 to 63 million years. These simulations show that RAD‐seq allows the recovery of the known Drosophila phylogeny with strong statistical support, even for relatively ancient nodes. Notably, this conclusion is robust to the potentially confounding effects of sequencing errors, heterozygosity, and low coverage. We further show that clustering RAD‐seq data using the BLASTN and SiLiX programs significantly improves the recovery of orthologous RAD loci compared with previously proposed approaches, especially for distantly related species. This study therefore validates the view that RAD sequencing is a powerful tool for phylogenetic inference.     

Phylogenetic relationships and taxonomic ranking of pipizine flower flies (Diptera: Syrphidae) with implications for the evolution of aphidophagy

The taxonomic rank and phylogenetic relationships of the pipizine flower flies (Diptera: Syrphidae: Pipizini) were estimated based on DNA sequence data from three gene regions (COI, 28S and 18S) and 111 adult morphological characters. Pipizini has been treated as a member of the subfamily Eristalinae based on diagnostic adult morphological characteristics, while the larval feeding mode and morphology is shared with members of the subfamily Syrphinae. We analysed each dataset, both separately and combined, in a total evidence approach under maximum parsimony and maximum likelihood. To evaluate the influence of different alignment strategies of rDNA 28S and 18S genes on the resulting topologies, we compared the topologies inferred from a multiple alignment using fast Fourier transform (MAFFT) program with those topologies resulting from aligning the secondary structure of these rDNA genes. Total evidence analyses resolved pipizines as a sister group of the subfamily Syrphinae. Although the structural alignment and the MAFFT alignment differed in the inferred relationships of some clades and taxa, there was congruence in the placement of pipizines. The homogeneous morphology of the Pipizini clade in combination with their unique combination of characters among the Syrphidae suggest a change of rank to subfamily. Thus, we propose to divide Syrphidae into four subfamilies, including the subfamily Pipizinae stat. rev.     

Cytochrome b and Bayesian inference of whale phylogeny

In the mid 1990s cytochrome b and other mitochondrial DNA data reinvigorated cetacean phylogenetics by proposing many novel and provocative hypotheses of cetacean relationships. These results sparked a revision and reanalysis of morphological datasets, and the collection of new nuclear DNA data from numerous loci. Some of the most controversial mitochondrial hypotheses have now become benchmark clades, corroborated with nuclear DNA and morphological data; others have been resolved in favor of more traditional views. That major conflicts in cetacean phylogeny are disappearing is encouraging. However, most recent papers aim specifically to resolve higher-level conflicts by adding characters, at the cost of densely sampling taxa to resolve lower-level relationships. No molecular study to date has included more than 33 cetaceans. More detailed molecular phylogenies will provide better tools for evolutionary studies. Until more genes are available for a high number of taxa, can we rely on readily available single gene mitochondrial data? Here, we estimate the phylogeny of 66 cetacean taxa and 24 outgroups based on Cytb sequences. We judge the reliability of our phylogeny based on the recovery of several deep-level benchmark clades. A Bayesian phylogenetic analysis recovered all benchmark clades and for the first time supported Odontoceti monophyly based exclusively on analysis of a single mitochondrial gene. The results recover the monophyly of all but one family level taxa within Cetacea, and most recently proposed super- and subfamilies. In contrast, parsimony never recovered all benchmark clades and was sensitive to a priori weighting decisions. These results provide the most detailed phylogeny of Cetacea to date and highlight the utility of both Bayesian methodology in general, and of Cytb in cetacean phylogenetics. They furthermore suggest that dense taxon sampling, like dense character sampling, can overcome problems in phylogenetic reconstruction.     

The quality of the fossil record and the accuracy of phylogenetic inferences about sampling and diversity

Because phylogenies can be estimated without stratigraphic data and because estimated phylogenies also infer gaps in sampling, some workers have used phylogeny estimates as templates for evaluating sampling from the fossil record and for "correcting" historical diversity patterns. However, it is not known how sampling intensity (the probability of sampling taxa per unit time) and completeness (the proportion of taxa sampled) affect the accuracy of phylogenetic inferences, nor how phylogenetically inferred estimates of sampling and diversity respond to inaccurate estimates of phylogeny. Both issues are addressed with a series of simulations using simple models of character evolution, varying speciation patterns, and various rates of speciation, extinction, character change, and preservation. Parsimony estimates of simulated phylogenies become less accurate as sampling decreases, and inaccurate trees chronically underestimate sampling. Biotic factors such as rates of morphologic change and extinction both affect the accuracy of phylogenetic estimates and thus affect estimated gaps in sampling, indicating that differences in implied sampling need not reflect actual differences in sampling. Errors in inferred diversity are concentrated early in the history of a clade. This, coupled with failure to account for true extinction times (i.e., the Signor-Lipps effect), inflates relative diversity levels early in clade histories. Because factors other than differences in sampling predict differences in the numbers of gaps implied by phylogeny estimates, inferred phylogenies can be misleading templates for evaluating sampling or historical diversity patterns.     

Molecular systematics of teioid lizards (Teioidea/Gymnophthalmoidea: Squamata) based on the analysis of 48 loci under tree‐alignment and similarity‐alignment

We infer phylogenetic relationships within Teioidea, a superfamily of Nearctic and Neotropical lizards, using nucleotide sequences. Phylogenetic analyses relied on parsimony under tree‐alignment and similarity‐alignment, with length variation (i.e. gaps) treated as evidence and as absence of evidence, and maximum‐likelihood under similarity‐alignment with gaps as absence of evidence. All analyses produced almost completely resolved trees despite 86% of missing data. Tree‐alignment produced the shortest trees, the strict consensus of which is more similar to the maximum‐likelihood tree than to any of the other parsimony trees, in terms of both number of clades shared, parsimony cost and likelihood scores. Comparisons of tree costs suggest that the pattern of indels inferred by similarity‐alignment drove parsimony analyses on similarity‐aligned sequences away from more optimal solutions. All analyses agree in a majority of clades, although they differ from each other in unique ways, suggesting that neither the criterion of optimality, alignment nor treatment of indels alone can explain all differences. Parsimony rejects the monophyly of Gymnophthalmidae due to the position of Alopoglossinae relative to Teiidae, whereas support of Gymnophthalmidae by maximum‐likelihood was low. We address various nomenclatural issues, including Gymnophthalmidae Fitzinger, 1826 being an older name than Teiidae Gray, 1827. We recognize three families in the arrangement Alopoglossidae + (Teiidae + Gymnophthalmidae). Within Gymnophthalmidae we recognize Cercosaurinae, Gymnophthalminae, Rhachisaurinae and Riolaminae in the relationship Cercosaurinae + (Rhachisaurinae + (Riolaminae + Gymnophthalminae)). Cercosaurinae is composed of three tribes—Bachiini, Cercosaurini and Ecpleopodini—and Gymnophthalminae is composed of three—Gymnophthalmini, Heterodactylini and Iphisini. Within Teiidae we retain the currently recognized three subfamilies in the arrangement: Callopistinae + (Tupinambinae + Teiinae). We also propose several genus‐level changes to restore the monophyly of taxa.     

Inferring phylogeny at low taxonomic levels: utility of rapidly evolving cpDNA and nuclear ITS loci

Plant molecular systematic studies of closely related taxa have relied heavily on sequence data from nuclear ITS and cpDNA. Positive attributes of using ITS sequence data include the rapid rate of evolution compared to most plastid loci and availability of universal primers for amplification and sequencing. On the other hand, ITS sequence data may not adequately track organismal phylogeny if concerted evolution and high rDNA array copy number do not permit identification of orthologous copies. Shaw et al. (American Journal of Botany 92: 142-166) recently identified nine plastid regions that appear to provide more potentially informative characters than many other plastid loci. In the present study, sequences of these loci and ITS were obtained for six taxonomic groups in which phylogenetic relationships have been difficult to establish using other data. The relative utility of these regions was compared by assessing the number of parsimony informative characters, character congruence, resolution of inferred trees, clade support, and accuracy. No single locus emerged as the best in all lineages for any of these measures of utility. Results further indicated that in preliminary studies, sampling strategy should include at least four exemplar taxa. The importance of sampling data from independent distributions is also discussed.     

A phylogeny of robber flies (Diptera: Asilidae) at the subfamilial level: molecular evidence

We present the first formal analysis of phylogenetic relationships among the Asilidae, based on four genes: 16S rDNA, 18S rDNA, 28S rDNA, and cytochrome oxidase II. Twenty-six ingroup taxa representing 11 of the 12 described subfamilies were selected to produce a phylogenetic estimate of asilid subfamilial relationships via optimization alignment, parsimony, and maximum likelihood techniques. Phylogenetic analyses support the monophyly of Asilidae with Leptogastrinae as the most basal robber fly lineage. Apocleinae+(Asilinae+Ommatiinae) is supported as monophyletic. The laphriinae-group (Laphriinae+Laphystiinae) and the dasypogoninae-group (Dasypogoninae+Stenopogoninae+Stichopogoninae+ Trigonomiminae) are paraphyletic. These results suggest that current subfamilial classification only partially reflects robber fly phylogeny, indicating the need for further phylogenetic investigation of this group.     

Phylogenetic relationships of African microhylid frogs inferred from DNA sequences of mitochondrial 12S and 16S rRNA genes

The phylogenetic relationships of microhylid frogs are poorly understood. The first molecular phylogeny for continental African microhylids is presented, including representatives of all subfamilies, six of the eight genera, and the enigmatic hemisotid Hemisus. Mitochondrial 12S and 16S rRNA sequence data were analysed using parsimony, likelihood and Bayesian methods. Analyses of the data are consistent with the monophyly of all sampled subfamilies and genera. Hemisus does not nest within either brevicipitines or non-brevicipitines. It is possibly the sister group to brevicipitines, in which case brevicipitines might not be microhylids. Phrynomantis and Hoplophryne potentially group with non-African, non-brevicipitine microhylids, in agreement with recent morphological and molecular data. Within brevicipitines, Breviceps is recovered as the sister group to a clade of Callulina+Spelaeophryne+Probreviceps. The relationships among the genera within this latter clade are unclear, being sensitive to the method of analysis. Optimal trees suggest the Probreviceps macrodactylus subspecies complex might be paraphyletic with respect to P. uluguruensis, corroborating preliminary morphological studies indicating that P. m. rungwensis may be a distinct species. P. m. loveridgei may be paraphyletic with respect to P. m. macrodactylus, though this is not strongly supported. Some biogeographic hypotheses are examined in light of these findings.     

Pushing short DNA fragments to the limit: Phylogenetic relationships of ‘hydrobioid’ gastropods (Caenogastropoda: Rissooidea)

Although phylogenetic studies are increasingly utilizing multi-locus datasets, a review of GenBank data for the Gastropoda indicates a strong bias towards a few short gene fragments (most commonly COI, LSU rRNA, and SSU rRNA). This is particularly the case for the Rissooidea, one of the largest and most taxonomically difficult gastropod superfamilies. Here we analyze fragments of these three genes from 90 species to determine whether they can well resolve higher relationships within this superfamily, whether structurally aligned sequence datasets increase phylogenetic signal, and whether the inclusion of highly variable regions introduces noise. We also used the resulting phylogenetic data in combination with morphological/anatomical evidence to re-evaluate the taxonomic status of ‘hydrobioid’ family-level groups.Our results indicate that all three of the alignment strategies that were used resulted in phylogenies having similar signal levels. However, there was a slight advantage to using structural alignment for inferring family-level relationships. Moreover, the set of ‘standard’ gastropod genes supported recognition of many previously recognized families and provides new insight into the systematics of several problematic groups. However, some family-group taxa were unresolved and the relationships among families were also poorly supported, suggesting a need for more extensive sampling and inclusion of additional genes.     

Molecular phylogeny and systematics of Dipodoidea: a test of morphology‐based hypotheses

Lebedev, V.S., Bannikova, A.A., Pagès, M., Pisano, J., Michaux, J.R. & Shenbrot, G.I. (2012). Molecular phylogeny and systematics of Dipodoidea: a test of morphology‐based hypotheses. —Zoologica Scripta, 42, 231–249. The superfamily Dipodoidea (Rodentia, Myomorpha) in its current interpretation contains a single family subdivided into six subfamilies. Four of them include morphologically specialized bipedal arid‐dwelling jerboas (Dipodinae – three‐toed jerboas, Allactaginae – five‐toed jerboas, Cardiocraniinae – pygmy jerboas and Euchoreutinae – long‐eared jerboas), the other two are represented by more generalized quadrupedal taxa (Zapodinae – jumping mice and Sminthinae – birch mice). Despite considerable effort from morphologists, the taxonomy as well as the phylogeny of the Dipodoidea remains controversial. Strikingly, molecular approach has never been envisaged to investigate these questions. In this study, the phylogenetic relationships among the main dipodoid lineages were reconstructed for the first time using DNA sequence data from four nuclear genes (IRBP, GHR, BRCA1, RAG1). No evidence of conflict among genes was revealed. The same robustly supported tree topology was inferred from the concatenated alignment whatever the phylogenetic methods used (maximum parsimony, maximum‐likelihood and Bayesian phylogenetic methods). Sminthinae branches basally within the dipodoids followed by Zapodinae. Monophyletic Cardiocraniinae is sister to all other jerboas. Within the latter, the monophyly of both Dipodinae and Allactaginae is highly supported. The relationships between Dipodinae, Allactaginae and Euchoreutinae should be regarded as unresolved trichotomy. Morphological hypotheses were confronted to findings based on the presented molecular data. As a result, previously proposed sister group relationships between Euchoreutes and Sicista, Paradipus and Cardiocraniinae as well as the monophyly of Cardiocaniinae + Dipodinae were rejected. However, the latter association is consistently supported by most morphological analyses. The basis of the obvious conflict between genes and morphology remains unclear. Suggested modifications to the taxonomy of Dipodoidea imply recognition of three families: Sminthidae, Zapodidae and Dipodidae, the latter including Cardiocraniinae, Euchoreutinae, Allactaginae and Dipodinae as subfamilies.     

Phylogenetic relationships of colubroid snakes based on mitochondrial DNA sequences

We used a 694 bp length of the mitochondrial ND4 gene from 40 genera to infer phylogenetic relationships among colubroid snakes. The goals of this study were to identify conserved subsets of ND4 sequence data that could be used to address (1) which nominal higher-level colubroid taxa are monophyletic, and (2) the relationships among the monophyletic lineages identified. Use of transversions only proved the most reliable and efficient means of retrieving colubroid relationships. Transversion parsimony and neighbour-joining analyses identify similar monophyletic higher-level taxa, but relationships among these lineages differ considerably between the two analyses. These differences were affected by the inclusion/exclusion of (1) transitions, (2) autapomorphies, and (3) the boid outgroups. Saturation effects among the transitions, uninformativeness of autapomorphies for clustering taxa, and long-branch and base-compositional problems among the boids lead us to regard the tree resulting from transversion parsimony analysis rooted with Acrochordus as the best current estimate of colubroid phylogenetic relationships. However, several aspects of this proposed phylogeny need further testing (e.g. the apparent diphyly of Natricinae is especially controversial). Relationships retrieved using all colubroid taxa are not obtained when sparsely or unevenly sampled experimental subsets of taxa are used instead, suggesting that long-branch problems can severely compromise elucidation of colubroid relationships if limited taxonomic sampling strategies are followed. We discuss the importance of this finding for previous molecular attempts to assess colubroid relationships. Our analyses confirm the historical validity of several nominal colubroid families and subfamilies, establish polyphyly of a few, but generally fail to resolve relationships among the monophyletic taxa we identify. More conservative character information will be required to confidently resolve the last issue.