The phylogenetic relationships of the members of the phylum Sipuncula are investigated by means of DNA sequence data from three nuclear markers, two ribosomal genes (18S rRNA and the D3 expansion fragment of 28S rRNA), and one protein-coding gene, histone H3. Phylogenetic analysis via direct optimization of DNA sequence data using parsimony as optimality criterion is executed for 12 combinations of parameter sets accounting for different indel costs and transversion/transition cost ratios in a sensitivity analysis framework. Alternative outgroup analyses are also performed to test whether they affected rooting of the sipunculan topology. Nodal support is measured by parsimony jackknifing and Bremer support values. Results from the different partitions are highly congruent, and the combined analysis for the parameter set that minimizes overall incongruence supports monophyly of Sipuncula, but nonmonophyly of several higher taxa recognized for the phylum. Mostly responsible for this is the split of the family Sipunculidae in three main lineages, with the genus Sipunculus being the sister group to the remaining sipunculans, the genus Phascolopsis nesting within the Golfingiiformes, and the genus Siphonosoma being associated to the Phascolosomatidea. Other interesting results are the position of Phascolion within Golfingiidae and the position of Antillesoma within Aspidosiphonidae. These results are not affected by the loci selected or by the outgroup chosen. The position of Apionsoma is discussed, although more data would be needed to better ascertain its phylogenetic affinities. Monophyly of the genera with multiple representatives (Themiste, Aspidosiphon, and Phascolosoma) is well supported, but not the monophyly of the genera Nephasoma or Golfingia. Interesting phylogeographic questions arise from analysis of multiple representatives of a few species.