Multiple genes elucidate the evolution of venomous snail-hunting Conus species

The species-rich Cone snails (Conus sp.) are predatory, marine gastropods known for small venom peptides that are valuable for pharmacological research applications. Phylogenetic analyses with mitochondrial rRNA sequences have facilitated peptide discovery. However, these relatively conserved genes leave unresolved the closer relationships among many species. We sequenced 26 internal transcribed spacer 2 (ITS2) sequences from genomic ribosomal DNA to elucidate the evolutionary relationships among molluscivorous species and to piscivorous and vermivorous species. We show that ITS2 sequences are well conserved within species but are sufficiently variable among species to resolve recent divergences. Using Bayesian, maximum likelihood and log-determinant methods, we use the ITS sequences to resolve portions of the tree that could not be resolved using the more conventional mt rRNA sequences. When the ITS2 sequences are added to existing COI and to the more conserved rRNA sequences and then properly modeled, support throughout the tree is increased. This enables us to show finer relationships among the molluscivorous species that reveal three well-supported clades (Conus, Cylinder, and Darioconus) and renders the ITS2 sequences an essential component in advancing the discovery and pharmacological characterization of novel peptides from the venoms of these molluscs.