Molecular Evolution, Functional Variation, and Proposed Nomenclature of the Gene Family That Includes Sphingomyelinase D in Sicariid Spider Venoms

Abstract

The venom enzyme sphingomyelinase D (SMase D) in the spider family Sicariidae (brown or fiddleback spiders [Loxosceles] and six-eyed sand spiders [Sicarius]) causes dermonecrosis in mammals. SMase D is in a gene family with multiple venom-expressed members that vary in functional specificity. We analyze molecular evolution of this family and variation in SMase D activity among crude venoms using a data set that represents the phylogenetic breadth of Loxosceles and Sicarius. We isolated a total of 190 nonredundant nucleotide sequences encoding 168 nonredundant amino acid sequences of SMase D homologs from 21 species. Bayesian phylogenies support two major clades that we name alpha and beta, within which we define seven and three subclades, respectively. Sequences in the alpha clade are exclusively from New World Loxosceles and Loxosceles rufescens and include published genes for which expression products have SMase D and dermonecrotic activity. The beta clade includes paralogs from New World Loxosceles that have no, or reduced, SMase D and no dermonecrotic activity and also paralogs from Sicarius and African Loxosceles of unknown activity. Gene duplications are frequent, consistent with a birth-and-death model, and there is evidence of purifying selection with episodic positive directional selection. Despite having venom-expressed SMase D homologs, venoms from New World Sicarius have reduced, or no, detectable SMase D activity, and Loxosceles in the Southern African spinulosa group have low SMase D activity. Sequence conservation mapping shows >98% conservation of proposed catalytic residues of the active site and around a plug motif at the opposite end of the TIM barrel, but alpha and beta clades differ in conservation of key residues surrounding the apparent substrate binding pocket. Based on these combined results, we propose an inclusive nomenclature for the gene family, renaming it SicTox, and discuss emerging patterns of functional diversification.