Evolution of the mimetic African swallowtail butterfly Papilio dardanus: molecular data confirm relationships with P. phorcas and P. constantinus

Sister group relationships of the African Mocker Swallowtail, Papilio dardanus, were determined using mitochondrial 16S rRNA and Cytochrome B sequences, and nuclear Elongation Factor EF-la and the ITS-1 region of the ribosomal RNA locus. All four data sets placed P. dardanus as the sister species of P. phorcas, with P. constantinus as the next closest relative, and quite distinct from other African Papilios such as P. nobilis. These data support earlier morphological studies indicating that the only two African Papilio species with multiple female colour forms have a common ancestor not shared with any other living species. This information is important for conclusions about the evolution of female-limited mimicry in swallowtail butterflies, and successful use of the ITS-1 gene opens up new possibilities for studying this phenomenon within P. dardanus.     

Phylogenetic analysis of restriction-site variation in wild and cultivated Amaranthus species (Amaranthaceae)

Amaranthus includes approximately 60 species, of which three are cultivated as a grain source. Many wild Amaranthus species possess agriculturally desirable traits such as drought and salt tolerance, and pathogen resistance. We examined relationships among wild and cultivated Amaranthus species based upon restriction-site variation in two chloroplast DNA regions and in a nuclear DNA region. The chloroplast regions consisted of (1) an intergenic spacer in transfer RNA genes and (2) the ribulose-1,5-bisphosphate carboxylase gene with a flanking open reading frame. The nuclear region was the internal transcribed spacers ITS-1 and ITS-2 flanking the 5.8S gene in the ribosomal DNA. These regions were amplified by the polymerase chain reaction and digested with a total of 38 restriction endonucleases. We detected 11 potentially informative restriction-site mutations and seven length-polymorphisms among the 28 Amaranthus species. Parsimony analysis was used to find the shortest tree for each separate data set (chloroplast, nuclear, and length) and for two combined matrices (chloroplast/nuclear and all data sets). Overall, there was a low level of variation which generated poorly resolved trees among the 28 species. Congruence analyses revealed that the chloroplast and nuclear data sets were congruent with each other but not to the length data set. The congruence of the chloroplast and nuclear data sets suggested that cytoplasmic gene flow may not be a confounding factor in our analyses. The phylogeny also suggested that drought tolerance evolved independently several times. The molecular phylogeny provides a basis for selection of species pairs for crop development.     

Phylogenetic study of heliconiine butterflies based on morphology and restriction analysis of ribosomal RNA genes

Relationships of ten heliconiine butterflies (genera Dryas and Heliconius , family Nymphalidae) were elucidated by phylogenetic analysis of characters based on ribosomal DNA restriction site variation and morphology. Agraulis vanillae , also a heliconiine, was used as the outgroup species. Although neither the morphological nor the molecular data unambiguously resolve relationships among the heliconiines, a combined analysis of both data sets results in a tree that is similar to traditional systematic arrangements and previous views of radiation in the group. Both pupal-mating and nonpupal-mating species group as clades in the combined analysis. However, the restriction site data alone do not support the monophyly of the pupal-mating clade, and the morphological data alone do not support the monophyly of the non-pupal-mating clade. Furthermore, relationships of H. melpomene, H. cydno and the silvaniform species depart from traditional arrangements based on morphology and reproductive compatibility experiments. All trees support the independent evolution of similar wing patterns of species previously suggested to be members of mimicry complexes. Several mimicry complexes appear to have a member in each of the two major monophyletic groups (pupal-mating and non-pupal-mating clades).     

Systematics and evolutionary relationships of the mountain lizard Liolaemus monticola (Liolaemini): how morphological and molecular evidence contributes to reveal hidden species diversity

The delimitation of species is a major issue in systematic biology and has been a re-emerging discipline in the last decade. A number of studies have shown that the use of multiple data sets is critical for the identification of cryptic species, particularly in groups with complex evolutionary histories. Liolaemus monticola is a montane lizard species distributed in central Chile (32°–42°S), with four described subspecies in a latitudinal gradient from north to south: L. m. monticola , L. m. chillanensis , L. monticola ssp. and L. m. villaricensis . In order to test the systematic status and phylogenetic relationships of the taxa included in the L. monticola group, we analysed morphological (morphometric and meristic) and molecular (allozyme and mitochondrial DNA) data sets. The results of the morphological analyses showed that meristic variables correctly assigned individuals with higher accuracy than did morphometric characters. The results of the analyses of allozyme data revealed eight diagnostic loci that are evidence for significant differences among the four L. monticola subspecies. Phylogenetic analyses with mitochondrial DNA data, including additional species, showed that the L. monticola group is polyphyletic. We postulate that the four current subspecies represent independent evolutionary lineages and must be raised to the specific level as L. monticola , L. chillanensis and L. villaricensis . The taxonomic status of the unnamed L. monticola ssp. remains unresolved, although we provide a preliminary proposal.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 635–650.     

Phylogenetic analysis of Hepatozoon species (Apicomplexa: Adeleorina) infecting frogs of Nova Scotia, Canada, determined by ITS-1 sequences

Species of Hepatozoon are apicomplexan parasites infecting tetrapod vertebrates and hematophagous arthropods. Two species, Hepatozoon catesbianae and Hepatozoon clamatae, have been described inhabiting the erythrocytes of bullfrogs and green frogs. A number of characteristics typically used to distinguish between members of this genus are shared between these 2 species, prompting speculation as to whether or not these organisms are in fact distinct species. To test the species distinction, bullfrogs and green frogs were captured at various sites across Nova Scotia, blood samples were collected, and DNA was extracted from samples containing parasites. The internal transcribed spacer 1 (ITS-1) from geographically diverse samples of both species was amplified by PCR, sequenced, and analyzed. ITS-1 sequences from the 2 species revealed single-nucleotide polymorphisms at 6 sites. Phylogenetic analysis of these molecular data and cytopathological features place isolates of each species in separate monophyletic groups. Comparison of the ITS-1 sequences between isolates from Nova Scotia and Ontario revealed that ITS-1 sequences of H. catesbianae from a previous study were mischaracterized as being those of H. clamatae. Phylogenetic data based on molecular variation and cytopathological features from this study provide the strongest evidence to date supporting the distinction between these 2 species.     

Untangling the phylogeny of neotropical lianas (Bignonieae, Bignoniaceae)

The tribe Bignonieae (Bignoniaceae) is a large and morphologically diverse clade of neotropical lianas. Despite being a conspicuous component of the neotropical flora, the systematics of the tribe has remained uncertain due to confusing patterns of morphological variation within the group. Chloroplast (ndhF) and nuclear (PepC) DNA sequences were used here to reconstruct the phylogeny of Bignonieae. Individual analyses of ndhF and PepC were highly similar to one another, yet localized differences in the placement of six species suggests some conflict between data sets. Combined analyses result in trees that are consistent with those from the individual analyses and provide greater support for the suggested relationships. This phylogeny provides important new insights into the systematics of the tribe. It identifies 21 strongly supported species groups, eight of which broadly correspond to currently recognized genera. In addition, each of these 21 species groups is supported by morphological synapomorphies. The consistency between morphological and molecular data suggests that the current phylogeny provides a solid framework for a formal revision of the generic-level classification and for addressing other aspects of the biology of Bignonieae.     

Combining DNA sequences and morphology in systematics: testing the validity of the dragonfly species Cordulegaster bilineata

Morphological and molecular techniques are rarely combined when answering questions of taxonomic validity. In this study, we combine morphological techniques with DNA sequences to determine the validity of the dragonfly species Cordulegaster bilineata. The two dragonfly species C. bilineata and C. diastatops are very similar in size, body color, and morphological characters, and due to these similarities, the status of C. bilineata as a valid species is in question. In this study we compare morphological measurements of males and internal transcribed spacer 1 (ITS-1) sequences of rDNA between the two taxa. The hamule measurements (where copulation occurs) of males show little difference between the taxa in question, but the anal appendage measurements (where the male first contacts the female) show marked divergence between the two taxa. Cluster analysis with these anal appendage measurements correctly assigns almost all individuals measured into their respective taxon. PCR amplification products of ITS-1 display a approximately 50 bp size difference between C. bilineata (n = 4) and C. diastatops (n = 5) regardless of collection site. Sequence data for these amplifications show 51 bp missing in one locus in the ITS-1 of C. bilineatarelative to C. diastatops. A lone population of C. diastatops from Wisconsin has three individuals with ITS-1 products that match the size of both C. bilineata and C. diastatops. One individual from this population appears to yield two ITS-1 amplification products that match both C. bilineata and C. diastatops. Although this population may be evidence for hybridization between the two taxa, such hybridization is not necessarily sufficient to disqualify the validity of a separate species designation for C. bilineata. Morphology and ITS-1 sequences depict a high degree of divergence that is consistent with species-level differences.     

COMBINED AND SEPARATE ANALYSES OF MORPHOLOGICAL AND MOLECULAR DATA IN THE PLANT FAMILY RUBIACEAE

Abstract— The Rubiaceae are one of the largest of the families of angiosperms, with over 10000 species. The tribal and subfamilial classification is provisional due to the lack of phylogenetic hypotheses. The present study of the Rubiaceae is based on 33 genera and three data sets, one morphological and two molecular from chloroplast DNA, restriction sites andrbcL sequences. There is much congruence between the morphological and the molecular data sets, but also conflict. For parsimony reasons, the best phylogenetic hypothesis is a tree based on an analysis of the combined data sets. The so-called “total evidence” criterion for the combined analysis is simply a reiteration of the principle of parsimony. In this particular study, the classification would be almost the same even if based on the separate analyses instead of the combined. Despite the inapplicability of consensus trees or trees from separate analyses for phylogenetic hypotheses and classification, separate analyses may provide important information. It is the best way to reveal conflicts between different data sets. Knowledge of the conflicts can promote further detailed investigation in order to improve understanding of characters and phylogenetic hypotheses. In this study, the tribe Vanguerieae provides such an example; morphological data support a position in the subfamily Cinchonoideae, but DNA and a tree based on the combined data support a position in subfamily Ixoroideae. The tribe's position in the morphological tree is probably due to missing information concerning the correct pollen presentation system. Bootstrap fractions and K. Bremer's branch support values are used to evaluate the stability of particular nodes in the trees. Interestingly these values are not always correlated, e.g. in the morphological tree, the node with the highest branch support value has very low bootstrap fraction. The reasons for these differences are unclear, but large differences are presumably more likely to occur on short branches.     

Phylogenetic analyses of marine sponges within the order Verongida: a comparison of morphological and molecular data

Abstract. Because the taxonomy of marine sponges is based primarily on morphological characters that can display a high degree of phenotypic plasticity, current classifications may not always reflect evolutionary relationships. To assess phylogenetic relationships among sponges in the order Verongida, we examined 11 verongid species, representing six genera and four families. We compared the utility of morphological and molecular data in verongid sponge systematics by comparing a phylogeny constructed from a morphological character matrix with a phylogeny based on nuclear ribosomal DNA sequences. The morphological phylogeny was not well resolved below the ordinal level, likely hindered by the paucity of characters available for analysis, and the potential plasticity of these characters. The molecular phylogeny was well resolved and robust from the ordinal to the species level. We also examined the morphology of spongin fibers to assess their reliability in verongid sponge taxonomy. Fiber diameter and pith content were highly variable within and among species. Despite this variability, spongin fiber comparisons were useful at lower taxonomic levels (i.e., among congeneric species); however, these characters are potentially homoplasic at higher taxonomic levels (i.e., between families). Our molecular data provide good support for the current classification of verongid sponges, but suggest a re-examination and potential reclassification of the genera Aiolochroia and Pseudoceratina . The placements of these genera highlight two current issues in morphology-based sponge taxonomy: intermediate character states and undetermined character polarity.     

Phylogenetic inferences in Antennaria (Asteraceae: Gnaphalieae: Cassiniinae) based on sequences from Nuclear Ribosomal DNA internal transcribed spacers (ITS)

The phylogenetic relationships among sexually reproducing species of Antennaria (Asteraceae) are poorly understood. An earlier cladistic analysis based on morphology did not fully resolve the phylogeny of these taxa and therefore a different approach using molecular data was explored. The internal transcribed spacer regions (ITS-1 and ITS-2) of nuclear ribosomal DNA were sequenced for 30 species of Antennaria and one species from each of the outgroup genera Anaphalis, Ewartia, Leontopodium, and Pseudognaphalium. The ITS-1 sequence in Antennaria ranged from 253 to 260 base pairs (bp) in length, and the proportion of nucleotide differences between pairs of species of Antennaria ranged from 1 to 14%. For ITS-2, the divergence between pairs of species of Antennaria ranged from 0 to 8%. ITS-2 is shorter than ITS-1, ranging from 213 to 219 bp. Phylogenetic analysis indicates that, relative to the outgroups included, Antennaria is a well-supported monophyletic group. Based on the genera surveyed, Leontopodium appears to be the sister genus of Antennaria. The general topology of the molecular trees agrees with that based on previous morphological analyses and indicates that Antennaria is composed of six clades of equal rank, corresponding to the traditionally recognized informal groups, the Geyeriae, Argenteae, Arcuatae, Dimorphae, Pulcherrimae, and Catipes. Sequence and morphological data indicate that the Alpinae and Dioicae are unnatural, polyphyletic units that should be abandoned and redefined as the monophyletic Catipes group. Phylogenetic analysis of ITS sequences also suggests the dissociation of A. stenophylla from the Dimorphae, where it is traditionally placed, and its affiliation with the Argenteae, as well as the placement of A. arcuata in its own group.     

Multigene Barcoding and Phylogeny of Geographically Widespread Muricids (Gastropoda: Neogastropoda) Along the Coast of China

The identification and phylogeny of muricids have been in a state of confusion for a long time due to the morphological convergence and plasticity. DNA-based identification and phylogeny methods often offer an analytically powerful addition or even an alternative. In this study, we employ a DNA barcoding method to identify 17 known and easily confused muricid species (120 individuals) from the whole China coast based on mitochondrial cytochrome c oxidase subunit I (COI) and 16S rRNA sequences, and nuclear ITS-1 and 28S rRNA sequences. The phylogeny of muricid subfamilies is also analysed based on all mitochondrial and nuclear sequences. The universal COI and 16S rRNA primers did not work broadly across the study group, necessitating the redesign of muricid specific COI and 16S rRNA primers in this paper. Our study demonstrates that COI gene is a suitable marker for barcoding muricids, which can distinguish all muricid species studied. Phylogenetic analysis of 16S rRNA, ITS-1 and 28S rRNA data also provide good support for the species resolution observed in COI data. The relationships of muricid subfamilies are resolved based on the separate and combined gene data that showed the monophyly of each the subfamilies Ergalataxinae, Rapaninae, Ocenebrinae and Muricinae, especially that Ergalataxinae did not fall within Rapaninae.     

Ecological and morphological differentiation among cryptic evolutionary lineages in freshwater limpets of the nominal form-group Ancylus fluviatilis (O.F. Müller, 1774)

The phylogeny and potential mode of speciation of the river limpet Ancylus fluviatilis (Basommatophora) was examined using mitochondrial DNA sequences from 16S ribosomal RNA, cytochrome c oxidase subunit I (COI) and nuclear DNA from internal transcribed spacer (ITS-1) regions from 103 populations across Europe. Four highly divergent lineages were observed within Ancylus. Clade 1, representing the nominal taxon Ancylus fluviatilis (O.F. Müller, 1774), is mainly found in central and northern Europe, Clade 2 is present in a single Portuguese population, Clade 3 is distributed on the Canary islands, North Africa and the eastern Mediterranean region, whereas Clade 4 inhabits the Northern Mediterranean coasts. Phylogenetic analyses revealed an overall consistent topology of nuclear and mitochondrial gene trees. Based on a molecular clock, we estimated that the basic radiation occurred in the late Pliocene. Although clades differ significantly in size independent shell shape, morphological differentiation of lineages is not feasible without genetic data. Environmental data related to climate (precipitation, temperature, etc.) showed a significant differentiation of clades. Clade 1 dwells in relatively colder and more stable habitats than Clades 3 and 4, whose habitats in turn differ in a low or high amount of precipitation during spring and autumn, respectively. Based on the combined data sets on mitochondrial DNA, nuclear DNA, morphological and ecological differentiation, we conclude that Ancylus represents a cryptic species complex of reproductively and genetically isolated lineages. In addition, the joint analysis suggests that ecological speciation is probable to explain current patterns.     

Nuclear rDNA-based molecular clock of the evolution of triatominae (Hemiptera: reduviidae), vectors of Chagas disease

The evolutionary history and times of divergence of triatomine bug lineages are estimated from molecular clocks inferred from nucleotide sequences of the small subunit SSU (18S) and the second internal transcribed spacer (ITS-2) of the nuclear ribosomal DNA of these reduviids. The 18S rDNA molecular clock rate in Triatominae, and Prosorrhynchan Hemiptera in general, appears to be of 1.8% per 100 million years (my). The ITS-2 molecular clock rate in Triatominae is estimated to be around 0.4-1% per 1 my, indicating that ITS-2 evolves 23-55 times faster than 18S rDNA. Inferred chronological data about the evolution of Triatominae fit well with current hypotheses on their evolutionary histories, but suggest reconsideration of the current taxonomy of North American species complexes.     

Comparative osteology of the Danio (Cyprinidae: Ostariophysi) axial skeleton with comments on Danio relationships based on molecules and morphology

Danio is a diverse group of small, colourful and easily bred freshwater fishes native to Southern Asia. Biological interest in danios has increased in recent years because the zebrafish, Danio rerio , has become an important model organism, particularly for studies of vertebrate developmental biology and genetics. Though several phylogenetic studies of Danio have been done on a subset of Danio species, the resulting phylogenies conflict in detail. To examine the utility of osteology for systematics of this group at the species level, we studied the axial skeleton for 11 species of Danio. We analyse our morphological data alone and combined with DNA sequence data for five gene sequences generated in earlier phylogenetic studies. The axial skeleton of Danio exhibits 14 characters that prove useful in phylogenetic analysis. Both molecular and morphological data support the monophyly of the danios included in our analysis and both data sets support the monophyly of two subclades: a deep-bodied group and a slender-bodied group. Separate analysis of molecular and morphological data sets show that the molecular data resolves relationships within the slender subclade whereas the topology of the deep-bodied subclade is determined by morphological data.  © 2002 The Linnean Society of London. Zoological Journal of the Linnean Society , 2002, 135 , 529–546.     

The nuclear rDNA region of Gyrodactylus arcuatus and G. branchicus (Monogenea: Gyrodactylidae)

The primary structure of the ribosomal DNA internal transcribed spacers (ITS-1 and ITS-2) and 5.8S rRNA gene were used to characterize and identify 2 monogenean species of Gyrodacrylus living externally on the threespine stickleback (Gasterosteus aculeatus). The ITS region was amplified by PCR from freshwater, brackish, and marine isolates of Gyrodactylus arcuatus and G. branchicus, and the ends of the coding regions were identified by comparative alignment. No intraspecific and very low interspecific variation were observed in the 5.8S rRNA gene; high inter- and low intraspecific variation were revealed in the ITS-1 and ITS-2 regions. The morphological species identification was in all cases confirmed by the molecular identification. Intraspecifically, samples from 2 locations in the North Sea could be differentiated, but the Baltic sample resembled North Sea genotypes. Our approach offers perspectives for a multimetric genetical, morphometrical, and ecological taxonomy of the genus Gyrodactylus.     

Evolutionary relationships of trichostrongyloid nematodes (Strongylida) inferred from ribosomal DNA sequence data

The evolutionary relationships of 21 species of trichostrongyloid nematodes were determined by use of sequence data of the second internal transcribed spacer of the ribosomal DNA aligned according to secondary structure information. Irrespective of the method of analysis used, the topologies of the phylogenetic trees derived from the molecular data differed with respect to all four hypotheses proposed previously for the evolutionary relationships of the different subfamilies within the Trichostrongylidae based on morphological data. Thus, the molecular data set did not resolve the conflict between the four previous proposals for the subfamilial relationships. Nonetheless, all trees derived from the molecular data showed strong support for the exclusion of the genera Filarinema and Amidostomum from the clade containing the species within the family Trichostrongylidae. This represents a major difference from the most recent proposal of the systematics of the Trichostrongyloidea in which these two genera were included within the Trichostrongylidae. Therefore, the molecular data support an earlier systematic framework in which Filarinema and Amidostomum were considered to be sister groups of the Trichostrongyloidea.     

Cryptic diversity,sympatry, and other integrative taxonomy scenarios in the Mexican <Emphasis Type="Italic">Ceratozamia miqueliana</Emphasis> complex (Zamiaceae)

According to the integrative taxonomy (IT) framework, delimitation of taxonomic specimens into bona fide species requires the incorporation of multiple sources of biological evidence that jointly perform a role as systematic criteria. Here, we collected a combination of geographic, morphological (qualitative and quantitative), and molecular data sets from a group of specimens in the Ceratozamia miqueliana species complex. We analyzed them under methods suitable for each of these partitions and then synthesized results under the taxonomic circle—i.e., an operational tool for inference in IT. Taxonomic circle-based species delimitations in this Ceratozamia species complex have downstream taxonomic-nomenclatural consequences given that morphometric affinity between these specimens does not allow their species-level identification. However, joint patterns of qualitative morphological (i.e., diagnostics in standard sense) and molecular (i.e., character-based, DNA diagnostic or “bar coding”) characters between the species suggest the following scenarios: (1) C. miqueliana is a morphologically diagnosable species that lacks DNA diagnostics; (2) Ceratozamia subroseophylla is diagnosable both in terms of morphological and molecular evidence and sympatric with C. miqueliana; (3) Ceratozamia euryphyllidia is diagnosable both in morphological and molecular terms; (4) Ceratozamia zoquorum and Ceratozamia santillanii are cryptic species, i.e., both have exclusive DNA diagnostics but are not diagnosable morphologically; and (5) Ceratozamia becerrae is not diagnosable by any criterion and should therefore fall in synonymy. The present work contributes to establish a robust platform for subsequent systematic assessments of Mexican cycad diversity, which could include new hypotheses on the evolutionary processes involved in their current biogeographic distribution, ecological relationships, and other life-history aspects, in turn useful for conservation biology concerns.     

Phylogenetic analysis of nuclear and mitochondrial DNA reveals a complex of cryptic species in Crassicutis cichlasomae (Digenea: Apocreadiidae), a parasite of Middle-American cichlids

We obtained nuclear ITS-1 and mitochondrial cox1 sequences from 225 Crassicutis cichlasomae adults collected in 12 species of cichlids from 32 localities to prospect for the presence of cryptic species. This trematode is commonly found in species of cichlids over a wide geographic range in Middle-America. Population-level phylogenetic analyses of ITS-1 and cox1, assessments of genetic and haplotype diversity, and morphological observations revealed that C. cichlasomae represents a complex of seven cryptic species for which no morphological diagnostic characters have been discovered thus far. Bayesian and Maximum Likelihood analyses of concatenated datasets (906 bp) recovered eight lineages of C. cichlasomae, all with high posterior probabilities and bootstrap branch support. Values of genetic divergence between clades ranged from 1.0% to 5.2% for ITS-1, and from 7.2% to 30.0% for cox1. Morphological study of more than 300 individuals did not reveal structural diagnostic traits for the species defined using molecular evidence. These observations indicate that some traditional morphological characters (e.g., testes position) have substantial intra-specific variation, and should be used with caution when classifying C. cichlasomae and their sister taxa. Additionally, phylogenetic analyses did not reveal a strict correlation between these cryptic species and their host species or geographic distribution, however it appears that genetic distinctiveness of these cryptic species was influenced by the diversification and biogeographical history of Middle-American cichlids.     

Reconstructing the phylogeny of the Sipuncula

Sipunculans are marine spiralian worms with possible close affinities to the Mollusca or Annelida. Currently 147 species, 17 genera, 6 families, 4 orders and 2 classes are recognized. In this paper we review sipunculan morphology, anatomy, paleontological data and historical affiliations. We have conducted cladistic analyses for two data sets to elucidate the phylogenetic relationships among sipunculan species. We first analyzed the relationships among the 45 species of Phascolosomatidea with representatives of the Sipunculidea as outgroups, using 35 morphological characters. The resulting consensus tree has low resolution and branch support is low for most branches. The second analysis was based on DNA sequence data from two nuclear ribosomal genes (18S rRNA and 28S rRNA) and one nuclear protein-coding gene, histone H3. Outgroups were chosen among representative spiralians. In a third analysis, the molecular data were combined with the morphological data. Data were analyzed using parsimony as the optimality criterion and branch support evaluated with jackknifing and Bremer support values. Branch support for outgroup relationships is low but the monophyly of the Sipuncula is well supported. Within Sipuncula, the monophyly of the two major groups, Phascolosomatidea and Sipunculidea is not confirmed. Of the currently recognized families, only Themistidae appears monophyletic. The Aspidosiphonidae, Phascolosomatidae and Golfingiidae would be monophyletic with some adjustments in their definition. The Sipunculidae is clearly polyphyletic, with Sipunculus nudus as the sister group to the remaining Sipuncula, Siphonosoma cumanense the sister group to a clade containing Siphonosoma vastumand the Phascolosomatidea, and Phascolopsis gouldi grouping within the Golfingiiformes, as suggested previously by some authors. Of the genera with multiple representatives, only Phascolosoma and Themiste are monophyletic as currently defined. We are aiming to expand our current dataset with more species in our molecular database and more detailed morphological studies.